CN103370635B - 用于确定时钟校正的方法和系统 - Google Patents

Abstract

基于伪距代码测量结果、载波相位测量结果、和由接收器网络提供的广播卫星时钟误差确定每个导航卫星的卫星时钟误差(S402)。确定每个卫星的计算卫星时钟误差与广播时钟误差之间的差值(S406)。对于每个集群，从导航卫星中选出时钟基准卫星，其中所述时钟基准卫星具有卫星集群的时钟误差差值的中值(S408)。通过向一个或多个集群中的每个卫星的广播时钟误差应用所述基准卫星的时钟误差的函数，确定所述广播时钟误差的校正(S410)。

Description

用于确定时钟校正的方法和系统

本文依据35U.S.C.119(e)要求于2011年1月14日提交的名称为“用于确定时钟校正的方法和系统”的第61/432,646号美国临时申请的优先权。

技术领域

本发明涉及用于确定卫星导航系统的时钟校正的方法和系统。

背景技术

定位接收器，例如全球导航卫星系统(GNSS)接收器，可估算物体或车辆的位置和速度。由于一个或多个卫星的不准确的时钟数据，定位接收器可能得出不精确的位置和速度。发送GNSS误差的校正信息的GNSS增强(augmentation)系统可以向用户提供能够获得比其他可能的系统更精确导航的信息。这种增强系统通常包括对发送的卫星时钟和轨道信息的时钟和轨道校正。

GNSS增强系统通常包括位于已知位置处的接收器网络。这些接收器向处理中心发送信息，在处理中心中，所述信息被合并并且GNSS校正信息被计算。由于增强系统网络中的每个接收器通常具有其自己的与其他网络时钟不同步的本地时钟，因此需要具有唯一时间源，使来自网络接收器的所有信息都可以参考该时间源。

现有技术中，在一个或多个网络接收器处，有精确的时钟(通常是网络接收器之外的原子钟)用于向网络接收器提供定时信号(timingsignal)。从而网络接收器具有非常稳定的定时参考，并且所述网络接收器可用作供来自其他网络接收器的信息参照的唯一的时间源。由于该精确的时钟可能会失效或该网络接收器与处理中心之间的通信连接可能失效，因此通常采用多个精确的时钟以提高系统的稳定性。这些精确的时钟的性能通常被实时地评估和检测，以便于在这些精确的时钟之间切换。

实践中，这种现有方法需要与网络接收器一同定位在多个不同位置处的昂贵的原子钟。因此，当前需要一种可用于确定卫星导航系统的时钟校正的改进的、经济的方法和系统。特别地，需要一种用于确定时钟校正的方法和系统，其在任意网络接收器处都不需要精确的时钟。

发明内容

根据一个实施例，一种用于确定卫星时钟校正的方法和系统有利于一个或多个移动定位接收器做出精确的位置确定。所述方法和系统包括通过由导航卫星的一个或多个集群传送的卫星信号的固定接收器网络收集伪距代码测量结果、载波相位测量结果和导航信息数据(例如，广播时钟误差或广播时钟偏差)。在每个时间间隔(例如，历元(epoch))中，基于所述伪距代码测量结果、载波相位测量结果和导航信息数据计算每个导航卫星的时钟误差。在每个导航卫星的计算出的时钟误差和广播时钟误差之间形成多个时钟误差差值。对于每个集群，从该集群中的导航卫星中选出时钟基准卫星，该时钟基准卫星具有所述卫星集群的多个时钟误差差值的中值。随后通过在所述集群中的每个卫星的广播时钟误差上增加所述时钟基准卫星的时钟误差的函数(例如，在每个集群中以逐个卫星的方式进行)，确定对该集群中的每个卫星的广播时钟误差的校正。每一历元或其他时间间隔(例如，多个历元)内的所述校正的值基于控制参数极限被限定(例如，其中，所述控制参数极限小于或近似等于所述时钟基准卫星的时钟误差，或者其中，在卡尔曼滤波器获得所述广播时钟误差的解之后，所述控制参数极限包括用距离单位表示的固定参数。

附图说明

图1是用于确定时钟校正的系统的第一实施例的方块图。

图2是用于确定时钟校正的系统的第二实施例的方块图。

图3是依据图1或2的系统的移动定位接收器的方块图。

图4是依据所述系统的任一实施例的用于确定时钟误差校正的方法的流程图的第一例。

图5是依据所述系统的任一实施例的用于确定时钟误差校正的方法的流程图的第二例。

图6是用于确定接收时钟误差校正的接收器的位置估算或速度估算的方法的流程图。

具体实施方式

依据用于确定卫星导航系统(例如，GNSS系统)的时钟校正的系统11的一个实施例，系统11有利于对广播卫星时钟误差数据的精确和平稳的卫星时钟校正进行实时的确定，以向可实时使用时钟校正数据44的移动定位接收器42的提供增强的可靠性和精确度。所述方法和系统生成时钟校正数据44(例如，时钟误差校正数据)，该数据非常适于在实时确定移动定位接收器42或相关物体的位置和速度时提供增强的精度和分辨率。

在一个实施例中，系统11包括能够通过一个或多个通信链路13与中央电子数据处理系统16进行通讯的接收器网络15。中央电子数据处理系统16被连接到基于地面的卫星通信装置38，该通信装置38可通过位于地球表面上方的轨道中的远程卫星通信装置40(例如，人造卫星)向一个或多个移动定位接收器42提供时钟校正数据44或其他校正数据(例如，位置校正数据)。移动定位接收器42与校正通信装置41(例如，卫星校正接收器)关联，所述校正通信装置41能够从远程卫星通信装置40接收卫星信号(例如，或时钟校正数据44)。尽管图1描绘了基于地面的卫星通信装置38、远程卫星通信装置40和用于通信或发送校正数据的校正通信装置41，但可通过电子通信、电信、电磁信号或无线通信中的任一方式发送校正数据或校正信息。

在一个实施例中，接收器网络15包括第一基准站10至第N基准站14，其中N等于大于2的任意正整数或整数。如图1所示，例如，接收器网络15包括第一基准站10、第二基准站12和第三基准站，其中N等于3。在一个实施例中，接收器网络15可包括全球地分布在不同的已知位置(例如，基准站的地理坐标)处的基准站(10，12和14)，尽管基准站(10，12和14)的区域的或本地的分布也是可行的。

每个基准站(10，12和14)包括位于相应已知或固定位置处的基准接收器51、基准电子数据处理系统53、和传送器或收发器55。每个基准接收器51接收或检测以下数据中的一个或多个：(1)伪距代码测量结果或数据，(2)载波相位测量结果或数据，(3)被播送的卫星时钟误差数据，(4)星历数据，和(5)其他基准导航数据。基准接收器51可从一个或多个全球定位系统(GPS)卫星、从一个或多个全球导航卫星系统(GLONASS)卫星、或同时从以上二者接收卫星信号、载波信号和编码信号。

接收器数据处理系统53可存储、检索和处理与接收到的卫星信号、载波信号和编码信号相关联的基准导航数据。在一个实施例中，每个基准站(10，12或14)的传送器或收发器55通过一个或多个通信链路13向中央电子数据处理系统16或计算机网络中心发送下述数据中的一个或多个：(1)伪距代码测量结果或数据，(2)载波相位测量结果或数据，(3)被播送的卫星时钟误差数据，(4)星历(ephemeris)数据，(5)基准导航数据，和(6)通过分析或处理任意前述数据项产生的衍生数据。例如，上述被传送的数据可涉及第一导航卫星集群(例如，全球定位系统(GPS)卫星)、第二卫星集群(例如，GLONASS卫星)，或同时涉及以上二者。

基准导航数据可包括下述数据项中的一个或多个：(a)每个相应基准站(10，12或14)的相应的基准站位置或地理坐标，(b)每个基准接收器51的相应的天线地理坐标，(c)每个基准站(10，12或14)的相应基准站标识符(identifier)，和(d)提供相应的广播时钟误差数据和相应的星历数据的相应卫星的卫星标识符。

通信链路13支持在基准站(10，12或14)和中央电子数据处理系统16之间的实时数据传送。通信链路13可包括以下一项或多项：数据包网络、数据网络、虚拟专用网络、互联网数据通道、电信设备、光学通信设备、无线射频通信设备、微波通信设备、光纤链路、点对点微波通信链路、专用线路连接、多导体电缆、同轴电缆传输线路、光传输线路、通信线路、或用于实时传输数据的其他适当的通信链路。

伪距代码测量结果或数据可以被基准接收器51从用于调制或编码由每个基准接收器51接收的卫星信号的信息中解码出来。

载波相位测量数据可通过基准接收器51中的相位测量装置获得，其中所述基准接收器例如通过接收器前端接收一个或多个载波信号。在一个实施例中，相位测量装置可提供GPS载波信号、GLONASS载波信号或二者的载波相位测量结果。

星历数据包括卫星轨道信息或按照计划时间参照地球表面上或附近的一个或多个基准位置的卫星位置的数据收集。星历数据可以被储存在数据库或一个或多个文件中。卫星可在能够被基准接收器51接收的导航信息中发送该星历数据。

在可替换的实施例中，基准接收器51检测与一个或多个载波信号相关的载波相位数据或载波相位测量结果；基准站(10，12或14)中的基准数据处理系统53确定用于基准接收器51接收的多个卫星信号的初始载波相位模糊度的解(solution)或载波相位模糊度的解集。基准数据处理系统53接收载波相位数据，并且初始模糊度解通过使用基准接收器51的已知的或固定的位置校正它，以确定增强的模糊度解。在一个例子中，校正后的载波相位数据包括估算的载波相位和增强的模糊度解，或者由它们衍生出的其他数据。依次地，校正后的载波相位数据或其他校正数据经由无线信号或电磁信号被传送至中央数据处理系统16。

在一个实施例中，中央电子数据处理系统16包括由一个或多个计算机构成的中央网络。中央电子数据处理系统16包括数据处理器18、数据存储装置22、和连接到数据总线20的通信界面36。数据处理器18可包括微处理器、控制器、可编程逻辑阵列、专用集成电路、算术运算逻辑单元、或适于处理数据、执行数据的数学运算、或执行数据的布尔逻辑或算术运算逻辑的其他装置。数据存储装置22可包括电子存储器、光学存储装置、磁性存储装置、或能够存储和找回数据的其他装置。如图1所示，数据存储装置22存储、支持或与下述模块中的一个或多个交互：数据收集模块24、参数估算器26、主预测性滤波器(predictive filter)27、时钟差分模块28、基准卫星选择器30、时钟校正模块32、和通信模块34。数据存储装置22中的任意上述模块(24，26，27，28，30，32和34)可包括软件模块、电子模块、软件指令、或用于执行等同的软件指令的硬件或电子模块。如果任意上述模块包括软件指令或程序，那么那些软件指令或软件程序可以被储存在数据存储装置22、光学存储介质、磁性存储介质或其他非暂时性或充分持久的存储介质中。

在一个实施例中，数据收集模块24提供软件指令用于或能够收集、组织、和管理由接收器网络15接收到的伪距代码测量结果或数据、载波相位测量结果或数据、和导航信息数据(例如，包括星历和广播卫星时钟误差数据)的处理，其中所述接收器网络15在多个频率上接收一个或多个导航卫星集群发送的信号。导航卫星集群可以例如是GPS集群，GLONASS集群，或其他集群。在一个实施例中，数据收集模块24追踪那些能够通过使用基准站的基准站标识符或地理坐标对基准数据进行识别、组织或分类的数据库、图表、文件(例如，反向文件)或其他数据结构中的上述信息。

参数估算器26提供软件指令用于或能够基于伪距测量结果或数据、载波相位测量结果或数据、基准站(或基准接收器)的相应的已知位置、和广播卫星时钟误差数据以逐历元的方式为每个卫星计算或确定计算卫星时钟误差数据(例如，时钟状态数据)。在一种构造中，参数估算器26可基于除伪距测量结果或数据、载波相位测量结果或数据、基准站(或基准接收器)的相应的已知位置之外的导航信息数据(例如，其包括广播卫星时钟误差数据)计算或确定计算卫星时钟误差。在一个实施例中，参数估算器26可包括、使用或使用以下器件中的一个或多个：中心轨道估算器、中心时钟估算器、预测性滤波器、降误(error-reducing)滤波器、和卡尔曼滤波器。

一个历元可由下述之一或多个定义：基准接收器51获取输入的射频或微波测距信号并产生伪距/载波相位测量结果时的起始时间、终止时间、持续或间隔、或离散的时刻。当前历元指的是产生最新测量结果(例如，GPS测量结果)的最近的时刻。

时钟差分模块28提供软件指令或能够形成或确定一个或多个集群中的每个卫星的成对的计算时钟误差数据和播送的时钟数据之间的时钟误差差值。在一个实施例中，对于由基准接收器51接收的每个卫星而言，可基于相位测量结果、伪距代码测量结果、基准站(10，12和14)的相应的已知固定位置、和广播卫星时钟误差确定或估算计算卫星时钟误差确定计算卫星时钟误差。时钟差分模块28可计算出与GPS集群有关的第一组计算卫星时钟误差数据，和用于GLONASS集群的第二组计算卫星时钟误差数据。

通常，每个卫星和定位接收器的绝对时钟在每个用于轨道和时钟确定的全球导航卫星系统(GNSS)(例如，GPS和GLONASS)中是不可观测的。需要至少一个额外的约束条件，例如每个全球导航卫星系统(GNSS)的伪距代码测量结果或载波相位测量结果，以精确地确定卫星时钟误差和接收器时钟误差。

在一个实施例中，广播卫星时钟误差数据可表现为相对于相应集群系统时钟(例如，GPS系统时钟或GLONASS系统时钟)的时钟偏差、时钟偏置或差分时钟数据。时钟差分模块28或数据处理系统16按照逐个卫星的方式确定每个集群中的每个卫星或仅限于每个集群中的特定的卫星子集内的每个卫星的计算卫星时钟误差数据与广播时钟数据之间的时钟误差差值或误差数据。时钟误差差值或误差数据可以被存储在数据存储装置22中的一个或多个阵列或矩阵中，在那里每个阵列对应着特定的集群。尽管每个阵列或矩阵可覆盖集群中的所有卫星，但也可能矩阵被限定为每个集群的子集，以改进在中央电子数据处理系统16中确定和汇集计算卫星时钟误差数据或时钟误差差值以实现向移动定位接收器42的终端用户提供校正数据(例如，时钟误差校正)的实时响应性所需的(例如，数据处理器18的)数据处理吞吐时间或计算时间。为此，特定的核心卫星，例如AMC-1卫星或相对于集群系统时间(例如，GPS系统时间)具有小于10纳秒精度的卫星时钟数据的历史精确度的其他卫星可以被包含在每个集群的卫星子集中，其中所述子集的差分或误差数据将被计算。

基准卫星选择器30提供软件指令用于或能够(例如，动态的以逐历元的方式)从每个集群中选出一个特定的卫星作为每个集群的时钟基准卫星，其中，对于特定的相应的历元而言，时钟基准卫星具有在每个卫星集群(例如，GPS和GLONASS)中按照上述方式形成的差值的中值。例如，该中值(而非平均值)可以用于排除任何错误的影响或来自外部卫星时钟数据对所生成的误差数据或差值的完全不同的冲击。如果使用平均值，则所生成的误差数据可能会被来自非正常工作的卫星时钟的外部卫星时钟数据不适当地歪曲。

时钟校正模块32提供软件指令用于或能够通过逐个卫星时钟地增加或调整广播卫星时钟数据计算或确定广播卫星时钟误差数据的校正信息(例如，时钟校正数据44或维修时钟校正数据)。例如，时钟校正模块32适于通过在集群中的其他卫星的每个播报卫星时钟误差上增加与选定的基准卫星有关的时钟差值或误差数据的函数(例如，数学函数，线性代数函数，或矩阵表达式)来确定时钟误差校正信息。在一个实施例中，时钟校正模块32包括下述一个或多个：预测性滤波器、降误滤波器和卡尔曼滤波器。在另一实施例中，时钟校正模块32和参数估算器26可利用、使用共享或包括主预测性滤波器27。主预测性滤波器27包括预测性滤波器、降误滤波器、或卡尔曼滤波器。

通信模块34提供软件指令用于或能够支持被计算出的时钟校正数据44向用户接收器的发送(例如，通过如图1所示的卫星通信系统，或者通过无线通信系统)。

在一个实施例中，轨道中的每个卫星上的一个或多个精准的卫星时钟(例如原子钟)被用作广播卫星时钟，用作潜在的主卫星时钟、或用作卫星时钟数据的来源。一个或多个卫星时钟可被选作用于校正数据的或移动定位接收器42位置分辨的时钟数据的源。例如，GPS卫星时钟、GLONASS卫星时钟、或GPS卫星时钟与GLONASS卫星时钟的混合可被选作时钟数据的源。在第一例子中，用在时钟校正数据44中的卫星时钟数据包括整个GPS集群或服务于广域的集群的具体部分(例如，主要的或统计学上重要的样本尺寸)中的GPS卫星时钟的中值。在第二例子中，用在时钟校正数据44中的卫星时钟数据包括整个GLONASS集群或服务于广域的集群的具体部分(例如，主要的或统计学上重要的样本尺寸)中的GLONASS卫星时钟的中值。在第三例子中，用在时钟校正数据中的卫星时钟数据包括整个GPS集群中的GPS卫星时钟的中值和整个GLONASS集群中的GLONASS卫星时钟的中值。

根据一个实施例，用于播报主时钟的系统和方法非常适于消除外精度时钟频率输入(例如，原子接收器时钟或原子基准振荡器)和降低全球或广域差分校正系统的日常操作复杂性和开发成本。如本文所述，主时钟应当指的是GPS集群和GLONASS集群的时钟基准卫星时钟，其对于任一历元都具有时钟误差差值的中值，使得主时钟在不同的历元中可以改变。由于在当前的方法和系统中大量卫星被用作可能的时钟数据的源，因此，时钟数据的可靠性可以被进一步提高。该系统和方法的另一优点为有助于具有零平均统计值(例如，接近零平均统计值)的时钟校正或时钟校正数据和可能受预建立或预设定的约束条件限制的校正漂移率(shift ratio)。

除了图2的系统111使用无线基准站138取代了接地卫星通信装置38和远程卫星通信装置40之外，图2的系统111与图1的系统11类似。无线基准站138可代表例如无线通信网络的一部分。此外，移动定位接收器42与校正通信装置141相关联，所述校正通信装置141包括用于与无线基准站138通信的收发器或无线通信装置(例如，蜂窝网络、全球移动通信系统(GSM)设备、码分多址(CDMA)设备或其他无线通信设备)。图1和图2中的类似的附图标记指代相同的元件。

在图2中，时钟校正数据44(例如，时钟误差校正数据)被传送给与无线基准站138的覆盖区域或无线基准站138的网络内的相应的校正通信装置41相关联的一个或多个移动定位接收器42，并使之可获得所述数据。

图3描绘了用在图1或图2的实施例中的移动定位接收器42的例子。移动定位接收器42包括被连接到模数转换器306的输入端的前端302。模数转换器306的输出端被连接到接收器电子数据处理系统308。

接收器电子数据处理系统308包括接收器数据处理器310、接收器数据存储装置312、接收器数据总线324、数据端口322、解码器320、和相位测量装置304。接收器数据处理器310、接收器数据存储装置312、数据端口322、相位测量装置304和解码器320被连接到接收器数据总线324，以支持在接收器电子数据处理系统308的前述部件之间的通信。

接收器数据处理器310可包括微处理器、控制器、可编程逻辑阵列、专用集成电路、算术运算逻辑单元、逻辑装置、数据处理电子装置，或用于执行软件、逻辑、算法、或程序指令的其他装置。

接收器数据存储装置312可包括电子存储器、光学存储装置、磁性存储装置、或其他存储装置。接收器数据存储装置312包括实时动态引擎314、预测性滤波器316(例如，卡尔曼滤波器)和位置估算器318。

定位接收器42与校正通信装置(41或141)相关联。校正通信装置(41或141)可以被集成到定位接收器42中或者可通过数据端口322与定位装置进行通信。校正通信装置(41或141)通过一个或多个中继的无线或卫星通信装置(例如，图1中的38、40，或者图2中的138)接收来自中央电子数据处理系统16的校正数据(例如，时钟误差校正数据(即，44)和/或基准载波相位校正数据)。

接收器前端302可包括任意适当的电路以接收由一个或多个卫星(例如，导航卫星)传送的卫星信号。例如，接收器前端302能够同时接收GPS和GLONASS卫星信号。接收器前端302可包括能够接收由卫星集群内的一个或多个卫星传送的多个载波信号的扩展频谱接收器或码分多址(CDMA)接收器。例如，接收器前端302可包括用于放大卫星信号的前置放大器或放大器，混合器和基准振荡器，其中放大器输入端被连接到天线，放大器输出端被连接到一个混合器输入端，基准振荡器被连接到另一个混合器输入端，以及混合器输出端被连接到接收器数据处理系统308。

在一个实施例中，模数转换器306在接收器前端302和接收器电子数据处理系统308之间提供了接口。模数转换器306将模拟相位测量结果转化为能够被接收器电子数据处理系统308处理或操作的数字载波相位测量数据。

解码器320确定伪距代码测量结果并向接收器电子处理系统308提供所述伪距代码测量结果。在一个实施例中，解码器320包括解调器或用于解调伪随机噪声代码(例如，粗略获取代码(coarse acquisition code)(C/A)或其他更精确的民用或军用编码)的其他装置，其中所述伪随机噪声代码用于调制一个或多个载波信号。例如，解码器320可包括一组相关器351(例如，每个GPS信道一个相关器，以及每个GLONASS信道一个相关器)，其中每个相关器351被连接到伪随机噪声代码生成器390以提供解调的同相(I)和90度移相(Q)的信号成分(例如，振幅和相位或向量)。对于GPS而言，L1载波信号被过程获取代码(C/A)和加密的精确代码P(Y)调制，而L2信号被加密的P(Y)代码调制。在一个实施例中，解码器可包括被连接到输入延迟模块的代码生成器，其中延迟模块的输出被连接到相关器，用于测量基准伪随机噪声代码(其可以被延迟模块用已知的增量延迟)与通过接收器前端302接收到的伪随机代码之间的相关性。解码器320还可有利于解码调制载波信号的导航信息，例如星历数据。

相位测量装置304包括用于测量载波信号的相位的任意装置、集成电路、电子模块或数据处理器。在一种构造中，相位测量装置304包括信号发生器，和连接到第一相关器和第二相关器的相位延迟模块。相位测量装置测量或估算由接收器前端302提供的一个或多个载波信号的被观测的相位。被测量的相位可以被表示为载波信号的整数波长、载波信号的分数波长、和/或载波信号的角度。

相位测量装置304可确定下述一项或多项：(1)第一载波信号、第二载波信号或两者的分数阶波长的第一测得相位成分；(2)第一载波信号、第二载波信号或两者的整个波长的第二测得相位成分。后面的第二测得相位成分可由计数器(例如，过零计数器)确定，其中，计数器对被接收、重构或处理后的载波信号在时域中以基准量级(如，0电压)与X轴相交的次数计数，其中，X代表时间，Y轴代表载波信号的量级。然而，相位测量装置304依赖于在接收器数据处理系统308中的进一步处理才能确定或分辨可能导致第二测得相位成分错误或被偏移整数个波长周期(例如，用于估算相应的卫星与定位接收器之间的距离或路程)的整个周期整周模糊度。相位测量装置304确定并向接收器电子处理系统308提供载波相位测量数据。

实时动态引擎314、预测性滤波器316、和位置估算器308可彼此通信。实时动态引擎314、预测性滤波器316、和位置估算器308可通过逻辑数据通路、物理数据通路或同时通过上述二者在彼此之间通信。逻辑数据通路指的是例如软件模块之间的或者一个或多个软件程序之间的虚拟的数据通路或数据通道。物理数据通路指的是传输线或支持数据、逻辑信号、电信号或电磁信号等通信的一个或多个接收器数据总线324。

在一个实施例中，实时动态引擎314包括搜索引擎或用于搜索或确定一个或多个接收自多个卫星的载波信号的相位的整周模糊度解集的其他软件指令。动态引擎314例如可检索由相位测量装置304提供的载波相位数据。整周模糊度解集指的是例如接收到的由一个或多个卫星发送的载波信号(例如，用于GPS的1.57542GHz的L1信号、1.22760GHz的L2信号或类似信号)的载波相位中的整数周期相位模糊度。动态引擎314或其搜索引擎可使用最小二乘方或卡尔曼滤波器技术以减少搜索空间或获得用于由卫星发送的载波信号的整数周期相位模糊度的一个或多个模糊度解集(ambiguity set solutions)。

在一种构造中，预测性滤波器316(例如，卡尔曼滤波器)能够促进搜索引擎的搜索更有效率、更快活或更精确地获得由卫星发送的载波信号的整数周期相位模糊度的模糊度解集。在可替换实施例中，替换性的不基于检索的技术(例如，宽莱茵相位载波组合和数据处理)可以用于解析载波相位整周模糊度解集。

估算器318包括接收用于估算与定位接收器42相关联的物体或车辆的位置、速度、姿态、或任一前述属性的接收器数据处理器310或其他数据处理装置。估算器318被连接到实时动态引擎314和接收器数据处理器310或与它们通信。一旦整周模糊度解集被确定，位置估算器318或定位接收器42可使用载波相位测量数据以基于已知的传播速度(即，光速)提供对于每个卫星和定位接收器之间的距离或路程的精确估算。然后，三个或更多个卫星与定位接收器之间的路程可用于估算接收器42的位置、速度或姿态。

预测性滤波器316包括降误滤波器、卡尔曼滤波器或其变形，用于减少或减轻误差，例如测量误差。卡尔曼滤波器可包括预测性滤波器装置或电路，其使用信号的叠加、延迟、和反馈来处理数据并补偿噪声的影响和测得数据中的不特定性等。重置或重新初始化所指的都是对降误或卡尔曼滤波器的状态进行重新初始化。

位置估算器318包括用于估算与定位接收器相关联的物体或车辆的位置的电子模块、软件模块、或以上二者。位置估算器可使用下述数据源中的一种或多种以确定定位接收器的天线或与之相关的物体或车辆的估算位置或姿态：解码的伪随机噪声代码、载波相位测量数据、加密的精确代码(例如，P(Y)代码)、粗略获取代码、导航信息、和整周相位模糊度数据、以及基准站载波相位数据，其中基准站载波相位数据可以被集成在整周相位模糊度数据中。

预测性滤波器316可包括接收来自实时动态引擎314的输入数据的降误滤波器，其中，所述输入数据包括用于相应载波相位测量数据的模糊度解数据(例如，整周模糊度集)。在一个实施例中，除了基准接收器51不连接到校正通信装置(41或141)之外，图1和图2中的基准接收器51可包括与图3的移动定位接收器42基本相同的元件。

图4是用于确定时钟校正数据44的方法的流程图的第一例。图4的方法开始于步骤S400。

在步骤S400中，数据收集模块24或中央电子数据处理系统16收集由基准站(例如，10，12和14)的接收器网络15接收到的伪距代码测量结果、载波相位测量结果、和播送的时钟误差。例如，基准站的接收器网络15包括定位在各个已知位置或已知地理坐标处的固定的基准接收器51，其中基准接收器51接收由一个或多个导航卫星集群发送的多个频率(例如，多个载波频率)的信号。所述一个或多个集群可包括GPS集群、GLONASS集群、或其他导航卫星集群。对于GPS卫星集群而言，多个频率例如可包括L1载波频率和L2载波频率。对于GLONASS卫星集群而言，多个载波频率被分配给不同的卫星。

在基准接收器51中，相位测量装置或接收器数据处理系统53测量被基准接收器51接收的相应载波信号的载波相位。载波相位测量结果例如可包括GPS载波相位测量数据和GLONASS载波相位测量数据。基准接收器51中的解码器可解码所传送的星历数据、广播卫星时钟误差、和被传输的或用于调制一个或多个卫星信号的至少一个载波信号的其他导航数据。在一种构造中，通过一个或多个通信链路13将载波相位测量结果、伪距代码测量结果、和导航信息数据从接收器15的网络传送给中央数据处理系统16。

每个卫星通常播送星历数据和卫星时钟误差数据(例如，对于特定卫星而言的相对于GPS系统时间或GLONASS系统时间的差分时钟数据)。广播时钟数据是供接收器使用以校正定位接收器(例如，基准接收器51)的伪距和载波相位测量结果的二阶多项式拟合因数的集合。每个导航卫星发送关于其自己的轨道和卫星时钟的星历数据。广播卫星时钟误差的精确度，例如相对于GPS时间而言，通常在2-3米水平之内。

在步骤S402中，在每个时间间隔(例如，历元)，参数估算器26或中央电子数据处理系统16基于以下收集到的数据中的一项或多项估算每个卫星的计算卫星时钟误差数据：每个卫星的接收器网络15的伪距代码测量结果、每个卫星的接收器网络15的载波相位测量结果、接收器网络15之内的固定的基准接收器51的相应已知的位置、和广播卫星时钟误差。

可根据各种技术执行步骤S402，这些技术可以被可替换地或累积地应用。根据第一技术，在每个时间间隔(例如，历元)，参数估算器26或中央电子数据处理系统16基于伪距代码测量结果、载波相位测量结果、固定的基准接收器51的相应已知的位置、和广播卫星时钟误差或导航信息数据(例如，其可包括广播卫星时钟误差)估算在一个或多个集群(例如，GPS、GLONASS或二者)中活动的或适当地起作用的每个卫星的卫星时钟误差数据。

根据第二技术，如果卫星在特定的时间周期内提供变坏的或不可靠的卫星信号，或变坏的或不可靠的数据，则参数估算器26或中央数据处理系统16可在一段时间周期内暂停处理该有缺陷的或可疑的卫星的卫星时钟误差数据，直到接收到的卫星信号具有所接收的数据的足够的信号质量或可靠性水平。

根据第三技术，参数估算器26或中央电子数据处理系统16使用单差分程序或其他程序确定一个或多个集群(例如，GPS或GLONASS集群)中的每个相应卫星的计算卫星时钟误差数据。例如，在一个实施例中，通过在基准接收器51的载波相位测量结果之间重复地应用单差分，每个可能的不同的接收器对、卫星组合或可用的交换(例如，在任一对基准接收器51处接收足够的信号强度和信号质量)被用于估算特定集群中的每个卫星的时钟误差或偏差(例如，计算卫星时钟误差)。计算卫星时钟误差数据可以被表述为相对于GPS系统时间或GLONASS系统时间的差分时间或偏差。

根据第四技术，时钟差分模块28或中央电子数据处理系统16使用二重差分程序或另一程序去确定计算卫星时钟误差数据。例如，二重差分程序指的是来自两个基准接收器51(例如，位于已知的固定位置或地理坐标处)的载波相位测量结果之间的差值，其中所述两个基准接收器51基本同时地接收来自同一对的两个卫星的载波信号，以确定计算卫星时钟误差数据具有降低的时钟误差或偏差。例如，在一个实施例中，依照基准接收器51的载波相位测量结果之间的二重差分的重复应用，每个可能的不同的接收器对-卫星对组合或可用的交换(例如，在任一对基准接收器51处接收足够的信号强度和信号质量)被用于估算特定集群中的每个卫星的时钟误差或偏差(例如，计算卫星时钟误差)。实践中，可通过获得两个单差分确定的第二个差分确定二重差分。计算卫星时钟误差数据可以被表述为相对于GPS系统时间或GLONASS系统时间的差分时间或偏差。

在步骤S406中，时钟差分模块28或中央电子数据处理系统16形成或确定每个卫星的相应的成对的计算卫星时钟误差数据与相应的广播时钟数据之间的时钟误差差值(例如，差分数据或差分时钟误差数据)。例如，计算卫星时钟误差数据是与特定基准站(例如，10，12或14)的基准接收器51的已知的固定为位置以及特定卫星的特定基准站(例如，10，12或14)的载波相位测量结果相一致的卫星时钟数据。

广播时钟数据由每个卫星上的精确的时钟(例如，铷、铯或原子钟)提供或由此得出。在一个实施例中，播送的时钟数据包括相对于集群系统时钟的差分的时钟数据或偏差时钟数据。在另一实施例中，播送的时钟数据代表被传送的卫星时钟数据(以及与之相关的星历数据)。

可根据能够独立地或累积地应用的各种技术执行步骤S406。

根据第一技术，对于一个或多个集群以及每个连续的历元，时钟差分模块28、数据处理系统16、中央电子数据处理系统16、或数据处理系统16中的主预测性滤波器27(例如，卡尔曼滤波器)使用一个或多个下述项逐个卫星地解析时钟误差差值：(1)同一集群中的一个或多个卫星的相应的计算卫星时钟误差数据，(2)同一集群中的一个或多个卫星的相应的广播卫星时钟误差数据，和(3)基于特定卫星的附带的星历数据经调节或校正后的同一集群中的相应的广播卫星时钟误差数据。所述同一集群指的是同一卫星集群，例如其可以包括GPS集群或GLONASS集群。

根据第二技术，时钟差分模块28或数据处理系统16逐个卫星地确定每个集群中的每个卫星、或者如以其它方式限于每个集群内的卫星的特定子集中的每个卫星的计算卫星时钟误差数据和广播时钟误差数据之间的差分或误差数据。差分或误差数据可以被存储在数据存储装置22中的一个或多个阵列或矩阵中，其中每个阵列对应着特定的集群。

在步骤S408中，基准卫星选择器30或电子数据处理系统16选择每个集群内的其中一个卫星作为该集群的时钟基准卫星，其中，时钟基准卫星具有该集群的时钟误差差值的中值(例如，参见步骤S406)。在执行步骤S408的第一例子中，基准卫星选择器30或电子数据处理系统16选择GPS集群中的其中一个卫星作为GPS集群的时钟基准卫星，其中，时钟基准卫星具有GPS集群的差分数据或时钟误差数据的中值。在执行步骤S408的第二例子中，基准卫星选择器30或电子数据处理系统16选择GLONASS集群中的其中一个卫星作为GPS集群的时钟基准卫星，其中，时钟基准卫星具有GLONASS集群的差分数据或时钟误差数据的中值。

在步骤S410中，时钟校正模块32或电子数据处理系统16通过调节一个或多个集群中的每个卫星的广播时钟误差数据确定或计算广播时钟数据的校正或校正数据。例如，时钟校正模块32或数据处理系统26通过向一个或多个集群中的每个卫星的广播时钟误差增加时钟基准卫星的时钟误差的函数，从而调节每个卫星的广播时钟误差数据，其中，每个历元或其他时间间隔(例如，多个相邻的历元)的校正的值可基于控制参数极限而被限定。

可依照那些可独立地或合并地应用的各种程序执行步骤S410。

在第一程序中，控制参数极限与时钟基准卫星的时钟误差成比例。在第二程序中，控制参数极限小于或近似等于时钟基准卫星的时钟误差。在第三程序中，在主预测性滤波器27(例如，卡尔曼滤波器)找到广播卫星时钟误差的校正的解之后，控制参数极限包括用距离单位表示的固定参数。

在第四程序中，对于一个或多个集群中的第一集群中的每个导航卫星而言，被计算出的校正在每个历元中被限定为第一集群的时钟基准卫星的第一时钟误差；并且，对于一个或多个集群中的第二集群中的每个导航卫星而言，被计算出的校正在每个历元中被限定为第二集群的时钟基准卫星的第二时钟误差，其中第一时钟误差与第二时钟误差是独立的。例如，第一集群可包括全球定位系统，以及其中第二集群可包括GLONASS集群。

在第五程序中，在找到广播时钟误差的校正的解之后，被计算出的校正在每个历元中被限定为不超出一阈值。在第六程序中，所述阈值大约为1毫米距离单位，或等同于0.00333纳秒时间单位，其中距离单位与时间单位之间的转换因数为光速(例如，3×108米/秒)。在第七程序中，时钟基准卫星的时钟误差被用作主预测性滤波器27(例如，卡尔曼滤波器)的额外的约束条件。

在第八程序中，时钟校正模块32或电子数据处理系统16通过将广播时钟误差的时钟误差校正限定为每个历元或其他时间间隔的最大时钟校正量来确定或计算广播时钟误差的校正数据，其中，所述其他时间间隔小于或等于下面的一项或多项：(a)相应的最大阈值(例如，近似等于或小于1毫米/历元或其他时间间隔，或者近似等于或小于0.00333纳秒)，(b)均方根误差(RMS)，或(c)距离单位或等同的时间单位中表达的相应的最大阈值，其中，光速是转换因数。例如，时钟误差校正数据中每单位时间内的变化可基于导航信息数据(例如，星历数据和广播时钟误差数据)加以限制，从而时钟误差校正数据或定位接收器42的临时校正数据偏移小于每单位时间的最大偏移量，其中所述最大偏移量可能导致与终端用户移动定位接收器42确定的位置估算有关的不连续性、矛盾、或不精确性。

在第九程序中，在规则的或周期性的间隔(例如，历元)中，时钟校正模块32或电子数据处理系统16通过动态地实时调整每单位时间内的时钟误差数据的变化，确定广播时钟误差数据的下一个时钟误差校正数据，其中所述变化小于或等于相应的最大阈值(例如，小于约1毫米/历元或其他时间间隔，或者小于约0.00333纳秒/历元或其他时间间隔)。

在第十程序中，可通过基于导航卫星的特定集群内的成对的计算卫星时钟误差数据与广播时钟误差数据(或者星历时钟数据)之间的中值差值定义时钟误差校正数据的每单位时间的变化的控制参数极限，执行步骤S140。例如，通过限定由中央数据处理系统16或时钟校正模块32确定的用于传送到定位接收器42的时钟误差校正(或最大校正偏移)，可如下执行步骤S410：

在方程1中，由参数估算器26计算出的卫星时钟误差数据与从接收器网络15接收到的星历时钟数据或广播卫星时钟误差数据之间的中值GPS时钟差值可根据下面的方式确定：

针对GPS，dClk_GPS＝中值(clockⁱ(hub)-clockⁱ(eph))，i＝1∶32   (1)

其中，dClk_GPS是GPS集群的中值GPS时钟差值，clockⁱ(hub)是第i个卫星的计算时钟误差数据，clockⁱ(eph)是第i个卫星的星历时钟数据或广播时钟误差数据，并且i是32个GPS卫星中的每个卫星的标识符。

在一个实施例中，在计算上述GPS集群的中值GPS时钟差值之前，计算卫星时钟误差数据和星历时钟数据或广播卫星时钟误差数据被正态化成当量单位以及相同的或基本相同的历元。如之后结合方程3和4的描述，根据中值GPS时钟差值的上述定义，dClk_GPS可被数据处理系统16或者其预测性滤波器(例如，卡尔曼滤波器)使用，以确定时钟误差校正数据(例如，时钟误差校正)的极限，或者每单位时间内的时钟误差校正数据中的最大变化，或者其在移动定位接收器42中的应用。

在方程2中，由参数估算器26计算出的卫星时钟误差数据与从接收器网络15接收到的星历时钟数据或广播卫星时钟误差数据之间的中值GLONASS时钟差值可根据下面的方式确定：

针对GLONASS，dClk_GLONASS＝中值(clock^j(hub)-clock^j(eph))，j＝1∶24  (2)

其中，dClk_GLONASS是GLONASS集群的中值GLONASS时钟差值，clock^j(hub)是第j个卫星的计算时钟误差数据，clock^j(eph)是第j个卫星的星历时钟数据或广播时钟误差数据，并且j是24个GLONASS卫星中的每个卫星的标识符。

在一个实施例中，在计算上述GPS集群的中值GLONASS时钟差之前，计算卫星时钟误差数据和星历时钟数据或广播卫星时钟误差数据被正态化成当量单位以及相同的或基本相同的历元。如之后结合方程3和4的描述，根据中值GLONASS时钟差值的上述定义，dClk_GLONASS可被数据处理系统16或者其预测性滤波器27(例如，卡尔曼滤波器)使用，以确定时钟误差校正数据(例如，时钟误差校正)的极限，或者每单位时间内的时钟误差校正数据中的最大变化，或者其在移动定位接收器42中的应用。

在一个实施例中，数据处理系统16、时钟校正模块32、或预测性滤波器(例如，卡尔曼滤波器)维持或确定时钟误差校正是平稳或连续的(例如，尽可能或实际可行的)以避免时钟误差校正过程中的任何可能的校正不连续或跳跃；并从而基于移动定位接收器42的时钟误差校正估算与速度相关的位置。参数估算器26、校正模块32、或预测性滤波器27在中央数据处理系统16处的约束条件更新(例如，预测性滤波器27的修正数据更新或状态更新)的每单位时间(例如，每个历元或其他时间间隔)内的最大时钟变化应当被限定为方程3：

$\mathrm{dX}=\frac{H^{T}P}{H^T\mathrm{PH}+R}\left(\mathrm{dClk}\right)<=\mathrm{Threshold}---\left(3\right)$

其中，dX是时钟校正估算的状态更新；H是包含灵敏度系数的设计矩阵；P是方差-协方差矩阵；dClk是时钟基准卫星的时钟误差，Threshold是以距离为单位(例如，1mm)或以等同的时间为单位的偏移的最大校正变化的阈值；R是所应用的校正数据的方差。

基于方程3，GPS和GLONASS卫星的播送主时钟的方差R或自适应方差(adaptive variance)可由主预测性滤波器27(例如，卡尔曼滤波器)或参数估算器26反向地计算，随后自适应方差(例如，用于GPS集群或GLONASS集群的自适应主时钟或基准卫星时钟的自适应校正数据)可被应用到主预测性滤波器27或参数估算器26中。

阈值方差水平T_R可依据方程4加以定义：

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

其中，T_R是方差R的阈值方差水平，并且其中，其他参数按照上文所述加以定义。

在一个实施例中，如果计算出的方差R大于阈值方差水平T_R，则时钟校正模块32或主预测性滤波器27使用计算出的方差R去限制或更新校正状态或时钟校正数据(例如，时钟误差校正)。在另一实施例中，如果计算出的方差R小于阈值方差水平T_R，则时钟校正模块32或主预测性滤波器27使用阈值方差水平T_R去限制或更新校正数据(例如，时钟误差校正)。

在可替换的实施例中，阈值方差水平T_R可包括预定的方差水平或用户可基于主时钟性能的实证研究、或基于移动定位接收器42的位置估算的历史精确度进行定义的恒定参数。例如，阈值方差水平可以被预设或预定为0.01m²。

参数估算器26(例如，中心轨道/时钟估算器)和由接收器网络15接收到的广播时钟误差(例如，星历时钟)之间的时钟误差差值的中值(例如，GPS中值或GLONASS中值)改进了可靠性，因为某些卫星可能(至少偶尔)会具有不正常的、无关的或有偏差的卫星时钟。如果输入主预测性滤波器27(例如，卡尔曼滤波器)或参数估算器26的约束过紧或过于受限，则它可能在产生的校正数据中导致显著的校正数据跳跃。如果输入主预测性滤波器27(例如，卡尔曼滤波器)或参数估算器26的约束过松，则由于卫星时钟与接收器(例如，基准接收器51)时钟之间的不可观测性，它可能在所产生的校正数据中导致慢漂移和对于时钟估算而言较大的方差。为了克服前述技术问题，根据一个实施例提出了一种自适应性的播送主时钟，其中，时钟校正中的变化被限制在特定的预定或预设的范围之内。

在步骤S412中，通信模块34、通信界面36或数据处理系统16将被计算出的时钟校正数据44(例如，时钟误差校正)通过卫星、无线或其他通信系统分配和传送给移动定位接收器42(例如，用户定位接收器)。在与图1所示的系统相一致的一个例子中，数据处理系统16通过基于地面的卫星通信装置38分配和传送时钟校正数据44，其中所述通信装置38提供了通向环绕地球运行的卫星上的远程卫星通信装置40的卫星上行通信信道。远程卫星通信装置40提供了与校正通信装置41通信的下行信道，其中所述校正通信装置41用于接收校正数据以输入移动定位接收器42。

在与图2所示的系统相一致的另一例子中，数据处理系统16通过无线基准站138向校正通信装置141或与移动定位接收器42相关联的传送器发送所述时钟校正数据44。

根据所述系统和方法，所有卫星时钟和定位接收器时钟的估算参数可以相对于约束条件而被定义，称之为特定集群的主时钟数据或时钟基准卫星时钟。对于实时全球差分校正系统或广域校正系统或本发明公开的系统和方法的任意实施例而言，电子数据处理系统都可以很好地适于产生平稳的时钟校正信息(例如，无瞬间跳跃)和有良好界限范围的时钟校正数据，以保持差分校正系统的各部分之间的通信链路的带宽并以最小的处理延时实时地预测时钟校正。

图5是用于确定时钟校正数据44的方法的流程图的第二例。图4和图5中类似的附图标记指代类似的过程或方法。图5的方法开始于步骤S401。

在步骤S401中，数据收集模块24或中央电子数据处理系统16收集由基准站(例如，10，12和14)的接收器网络15接收到的伪距代码测量结果、载波相位测量结果、和导航信息数据(例如，星历数据和时钟数据)。例如，基准站的接收器网络15包括位于相应的已知的位置或已知的地理坐标上的固定的基准接收器51，其中，基准接收器51接收由一个或多个导航卫星集群传送的多个频率(例如，多载波频率)的信号。一个或多个集群可包括GPS集群、GLONASS集群、或导航卫星的其他机器。对于GPS卫星集群而言，多频率例如可包括L1载波频率和L2载波频率。对于GLONASS卫星集群而言，多载波频率可以被分配给不同的卫星。

在基准接收器51中，相位测量装置或接收器数据处理系统53测量基准接收器51接收到的相应载波信号的载波相位。载波相位测量例如可包括GPS载波相位测量数据和GLONASS载波相位测量数据。基准接收器51中的解码器可解码所传送的星历数据、卫星时钟数据、和被传送或调制一个或多个卫星信号中的至少一个载波信号的其他导航数据。载波相位测量结构、伪距代码测量结果、和导航信息数据被从接收器网络15通过一个或多个通信链路13传送给中央数据处理系统16。

每个卫星通常播送广播星历数据和被播送的卫星时钟误差数据(例如，与特定卫星的GPS系统时间或GLONASS系统时间相对的差分时钟数据)。广播星历数据是由接收器使用用于预测定位接收器(例如，基准接收器51)的接收器时钟行为的一组二阶多项式拟合系数。每个导航卫星传送绕其自身轨道的星历数据和卫星时钟。例如，相对于GPS时间，播送星历时钟的精确度通常在2-3米之内。

在步骤S403中，在每一个时间间隔(例如，历元)内，参数估算器26或中央电子数据处理系统16基于以下被收集的数据中的一项或多项计算每个卫星的卫星时钟误差数据：每个卫星的接收器网络15的伪距代码测量结果、每个卫星的接收器网络15的载波相位测量结果、接收器网络15内的固定的基准接收器51的相应的已知位置、和导航信息数据(例如，星历数据和广播卫星时钟误差数据)。

可根据各种技术执行步骤S403，这些技术可以被可替换地或合并地应用。根据第一技术，在每个时间间隔(例如，历元)中，参数估算器26或中央电子数据处理系统16基于伪距代码测量结果、载波相位测量结果、固定的基准接收器51的相应已知的位置、和导航信息数据估算在一个或多个集群(例如，GPS、GLONASS或二者)中活动的或适当地起作用的每个卫星的卫星时钟误差数据。

在第一程序中，控制参数极限与时钟基准卫星的时钟误差成比例。在第二程序中，控制参数极限小于或近似等于时钟基准卫星的时钟误差。在第三程序中，在主预测性滤波器27(例如，卡尔曼滤波器)找到广播卫星时钟误差的校正的解之后，控制参数极限包括用距离单位表示的固定参数，。

根据第三技术，参数估算器26或中央电子数据处理系统16使用单差分程序或其他程序确定一个或多个集群(例如，GPS或GLONASS集群)中的每个相应卫星的计算卫星时钟误差数据。例如，单差分过程指的是两个基准接收器51(例如，处于已知的固定位置或地理坐标处)接收到的载波相位测量结果之间的差值，所述两个基准接收器51基本上同时接收来自相同的卫星的载波信号，以确定具有减少的时钟误差或偏差的计算卫星时钟误差数据。计算卫星时钟误差数据可以被表述为相对于GPS系统时间或GLONASS系统时间的差分时间或偏差。

根据第四技术，时钟差分模块28或中央电子数据处理系统16使用二重差分程序或另一程序去确定计算卫星时钟误差数据。例如，二重差分程序指的是来自两个(例如，位于已知的固定位置或地理坐标出的)基准接收器51的载波相位测量结果之间的差值，其中所述两个基准接收器51基本同时地接收来自同一对的两个卫星的载波信号以确定具有减少的时钟误差或偏差的计算卫星时钟误差数据。实践中，可通过获得两个单差分确定的第二个差分确定二重差分，其中，单差分指的是两个基准接收器51(例如，处于已知的固定位置或地理坐标处)接收到的载波相位测量结果之间的差值，所述两个基准接收器51基本上同时接收来自相同的卫星的载波信号，以确定被具有减少的时钟误差或偏差的计算卫星时钟误差数据。计算卫星时钟误差数据可以被表述为相对于GPS系统时间或GLONASS系统时间的差分时间或偏差。

图5中的其他步骤或过程，包括步骤S406，S408，S410和S412与图4的方法中的那些相同或基本类似。

图6是用于在接收时钟误差校正的接收器处确定位置估算或速度估算的方法的流程图。图6的方法开始于步骤S600。图6的方法可使用图3中所描绘的定位接收器的实施例或者能够接收和使用时钟误差校正以提供比其他可能方式更精确的位置估算的定位接收器的其他构造。

在步骤S600中，解码器320或接收器电子数据处理系统308对由一个或多个导航卫星发送的、相应载波信号上的、由定位接收器42接收到的伪距代码测量结果进行解码。

在步骤S602中，相位测量装置304或接收器电子数据处理系统308确定相应的载波信号的载波相位测量结果。

在步骤S604中，解码器320或接收器电子数据处理系统308接收来自一个或多个导航卫星的载波信号上的广播卫星时钟误差。

在步骤S606中，校正通信装置(41或141)或接收器电子数据处理系统308的数据端口322接收对于每个导航卫星的广播时钟误差的校正，所述校正包括被添加到一个或多个集群中的每个导航卫星的广播时钟误差上的时钟基准卫星的时钟误差的函数，其中，基于控制参数极限限制每历元的校正的值。

在步骤S608中，数据处理系统308或位置估算器308基于被解码的伪距代码测量结果、被确定的载波相位测量结果、接收到的广播卫星时钟误差和接收到的校正确定定位接收器的估算位置。例如，数据处理系统308或位置估算器308可分辨所述被确定的载波相位测量结果中的相位模糊度或整数相位模糊度，以获得与所接收到的广播时钟误差相一致的位置估算。

可依据各种技术执行步骤S608，这些技术可以被单独地或者合并地适用。

根据第一技术，数据处理系统308或位置估算器318可基于被解码的伪距代码测量结果、被确定的载波相位测量结果、接收到的广播卫星时钟误差和接收到的时钟误差校正(例如，以逐历元的方式)确定定位接收器的估算位置和估算速度。

在第二技术中，控制参数极限与时钟基准卫星的时钟误差成比例。在第三技术中，控制参数极限小于或近似等于时钟基准卫星的时钟误差。在第四技术中，每个历元或其他时间间隔(例如，多个相邻的历元)之内的校正或时钟校正中的变化由下面的方程限定：

$\mathrm{dX}=\frac{H^{T}P}{H^T\mathrm{PH}+R}\left(\mathrm{dClk}\right)<=\mathrm{Threshold}$

其中，dX是时钟校正估算的状态更新；H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；R是所应用的校正数据的方差；dClk是时钟基准卫星的时钟误差；以及Threshold是如在权利要求中定义的每历元所容许的最大校正变化。

在第五技术中，如果所应用的校正数据的计算出的方差R大于下文指定的值T_R，则计算出的方差R被用于校正状态的更新。

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

在第六技术中，如果计算出的方差R小于值T_R，则T_R被用于替代R，用于依据第四和第五技术中提供的方程的校正状态的更新。

已经描述了优选的实施例，但将变得明显的是，在不背离本发明的由所附权利要求界定的范围的前提下，可以做出各种修改。

Claims (25)

1.一种确定卫星时钟校正以允许至少一个定位接收器提供位置估算或速度估算或时间估算中的至少一项的方法，所述方法包括下述步骤：

从由导航卫星的一个或多个集群传送的卫星信号的固定接收器网络收集伪距代码测量结果、载波相位测量结果和广播卫星时钟误差；

在每个历元处，基于所述伪距代码测量结果、载波相位测量结果和广播卫星时钟误差计算每个导航卫星的计算卫星时钟误差；

在每个导航卫星的相应成对的计算卫星时钟误差和广播卫星时钟误差之间形成多个时钟误差差值；

对于每个集群，从该集群中的导航卫星中选出时钟基准卫星，该时钟基准卫星具有所述卫星集群的所述多个时钟误差差值的中值；以及

通过在所述一个或多个集群中的每个导航卫星的广播时钟误差上增加所述时钟基准卫星的时钟误差的数学函数，计算针对每个导航卫星的广播时钟误差的校正，其中，基于控制参数极限限制每一历元或其他时间间隔内的所述校正的值。

2.根据权利要求1所述的方法，其中，所述控制参数极限与所述时钟基准卫星的时钟误差成比例。

3.根据权利要求1所述的方法，其中，所述控制参数极限小于或等于所述时钟基准卫星的时钟误差。

4.根据权利要求1所述的方法，其中，在预测性滤波器获得所述广播卫星时钟误差的校正的解之后，所述控制参数极限包括用距离或时间单位表示的固定参数。

5.根据权利要求1所述的方法，其中，对于所述一个或多个集群中的第一集群中的每个相应的导航卫星，计算出的校正在每个历元中被限定为所述第一集群的时钟基准卫星的第一时钟误差；并且，其中，对于所述一个或多个集群中的第二集群中的每个相应的导航卫星，计算出的校正在每个历元中被限定为所述第二集群的时钟基准卫星的第二时钟误差，所述第一时钟误差独立于所述第二时钟误差。

6.根据权利要求5所述的方法，其中，所述第一集群包括全球定位系统集群，并且其中所述第二集群包括GLONASS集群。

7.根据权利要求1所述的方法，其中，在获得所述广播时钟误差的校正的解之后，计算出的校正在每一历元内被限定为不超出一阈值。

8.根据权利要求7所述的方法，其中，所述阈值为一毫米或等同的时间阈值。

9.根据权利要求1所述的方法，其中，应用所述时钟基准卫星的时钟误差作为所述控制参数极限被执行为预测性滤波器的额外的约束条件。

10.根据权利要求1所述方法，其中，每一历元中的所述校正或时钟校正中的变化由下述方程限定：

$\mathrm{dX}=\frac{H^{T}P}{H^T\mathrm{PH}+R}\left(\mathrm{dClk}\right)<=\mathrm{Threshold}$

其中，dX是时钟校正估算的状态更新；H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；R是所应用的校正数据的方差；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

11.根据权利要求1所述的方法，其中，在所应用的校正数据的计算出的校正数据的方差R大于依据下述方程定义的阈值方差水平T_R时，将计算出的方差R用于校正状态的更新：

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

其中，H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

12.根据权利要求1所述的方法，其中，在计算出的方差R小于依据下述方程定义的阈值方差水平T_R时，T_R替代R被用于校正状态的更新：

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

其中，H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

13.一种确定卫星时钟校正以允许至少一个定位接收器提供位置估算或速度估算或时间估算中的至少一项的系统，所述系统包括：

数据收集模块，该数据收集模块用于通过由导航卫星的一个或多个集群传送的卫星信号的固定接收器网络收集伪距代码测量结果、载波相位测量结果和广播卫星时钟误差；

参数估算器，该参数估算器基于所述伪距代码测量结果、载波相位测量结果和广播卫星时钟误差计算每个导航卫星的计算卫星时钟误差；

时钟差分模块，该时钟差分模块在每个导航卫星的相应成对的计算卫星时钟误差和广播卫星时钟误差之间形成多个时钟误差差值；

基准卫星选择器，该基准卫星选择器用于从导航卫星中选出时钟基准卫星，其中所述时钟基准卫星具有所述卫星集群的多个时钟误差差值的中值；以及

时钟校正模块，该时钟校正模块通过在所述一个或多个集群中的每个导航卫星的广播时钟误差上增加所述时钟基准卫星的时钟误差的数学函数，计算针对每个导航卫星的广播时钟误差的校正，其中，每一历元或其他时间间隔内的所述校正的值基于控制参数极限被限定。

14.根据权利要求13所述的系统，其中，所述控制参数极限与所述时钟基准卫星的时钟误差成比例。

15.根据权利要求13所述的系统，其中，所述控制参数极限小于或等于所述时钟基准卫星的时钟误差。

16.根据权利要求13所述的系统，其中，在预测性滤波器获得所述广播卫星时钟误差的校正的解之后，所述控制参数极限包括用距离或时间单位表示的固定参数。

17.根据权利要求13所述的系统，其中，对于所述一个或多个集群中的第一集群中的每个相应的导航卫星，计算出的校正在每个历元中被限定为所述第一集群的时钟基准卫星的第一时钟误差；并且，其中，对于所述一个或多个集群中的第二集群中的每个相应的导航卫星，计算出的校正在每个历元中被限定为所述第二集群的时钟基准卫星的第二时钟误差，所述第一时钟误差独立于所述第二时钟误差。

18.根据权利要求17所述的系统，其中，所述第一集群包括全球定位系统集群，并且其中所述第二集群包括GLONASS集群。

19.根据权利要求13所述的系统，其中，在获得所述广播时钟误差的校正的解之后，计算出的校正在每一历元内被限定为不超出一阈值。

20.根据权利要求19所述的系统，其中，所述阈值为一毫米或等同的时间阈值。

21.根据权利要求13所述的系统，其中，所述校正模块包括预测性滤波器，该预测性滤波器适于接收所述时钟基准卫星的时钟误差作为额外的约束条件。

22.根据权利要求13所述的系统，其中，每一历元中的所述校正或时钟校正中的变化由下述方程限定：

$\mathrm{dX}=\frac{H^{T}P}{H^T\mathrm{PH}+R}\left(\mathrm{dClk}\right)<=\mathrm{Threshold}$

其中，dX是时钟校正估算的状态更新；H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；R是所应用的校正数据的方差；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

23.根据权利要求13所述的系统，其中，在时钟校正模块的主预测性滤波器中，在所应用的校正数据的计算出的校正数据的方差R大于依据下述方程定义的阈值方差水平T_R时，将计算出的方差R用于校正状态的更新：

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

其中，H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

24.根据权利要求13所述的系统，其中，在时钟校正模块的主预测性滤波器中，在计算出的方差R小于依据下述方程定义的阈值方差水平T_R时，T_R替代R用于校正状态的更新：

$T_R=\frac{H^{T}P}{\mathrm{Threshold}}\left(\mathrm{dClk}\right)-H^T\mathrm{PH}$

其中，H是具有灵敏度系数的设计矩阵；P是方差-协方差矩阵；dClk是时钟基准卫星的时钟误差；以及Threshold是每一历元中的最大校正变化或所述控制参数极限。

25.根据权利要求13所述的系统，其中，所述校正模块包括预测性滤波器，其中，所述预测性滤波器一旦初始化，则用于被播送的时钟数据的校正数据、计算出的校正数据被限定为每个时间间隔的不超过时钟基准卫星的时钟误差的最大值。