专利名称：用以估计mw偏差的gnss信号处理的制作方法
技术领域：
本发明涉及GNSS信号处理，并且具体涉及包含精确卫星数据的GNSS信号处理。
发明内容
描述了用于处理从多个接收机处对多个历元上的多个卫星的GNSS信号的码观测结果和载波-相位观测结果得到的一组GNSS信号数据的方法和设备，所述GNSS 信号具有至少两个载波频率，包括形成在每个历元处针对每个接收机-卫星对的丽(Melbourne-WUbbena)组合，以获得针对每个历元的丽数据组，以及从针对每个历元的MW数据组估计针对每个卫星的MW偏差和一组WL (宽巷)�：龋鯩W偏差可从一个历元到另一个历元地变化，每个WL�：榷杂τ诮邮栈�-卫星链路和卫星-接收机-卫星链路中的一个，其中，针对每个卫星的所述MW偏差被建模成以下之一 (i)单个估计参数和 ( )所估计的偏移加上具有所估计的振幅的谐波变化。在一些实施例中，根据1的方法，校正被施加到所述GNSS信号数据。在一些实施例中，在估计针对每个卫星的MW偏差之前，GNSS信号数据的至少一个线性组合被平滑。在一些实施例中，施加至少一个MW偏差约束。在一些实施例中，施加至少一个整数WL模糊度约束。在一些实施例中，在观测图和滤波图之一上使用生成树(ST)以约束所述莫糊度。 在一些实施例中，在观测图和滤波图之一上使用最小生成树(MST)以约束所述莫糊度。在一些实施例中，所述最小生成树基于边权重，每个边权重是从接收机-卫星几何得到的。在一些实施例中，所述边权重是关于以下之一来限定的(i)从接收机到卫星的几何距离，(ii)卫星仰角角度，以及(iii)从-卫星到接收机到卫星的几何距离，以及(iv) 仰角的组合，在所述仰角下，在站处可见单差组合中的两个卫星。在一些实施例中，所述最小生成树基于边权重，每个边权重都基于WL�：刃畔ⅲ⑶沂枪赜谝韵轮焕聪薅ǖ� (i)m4莫糊度到整数的差，(ii)m4莫糊度方差，以及(iii)⑴和(ii)的组合。在一些实施例中，所述莫糊度中的至少一个被固定为整数值。在一些实施例中，确定WL整数�：戎档亩喔龊蜓∽椋纬伤龆喔龊蜓∽榈募尤ㄗ楹希约敖鯳L �：戎械闹辽僖桓龉潭ㄎ∽运黾尤ㄗ楹系闹�。在一些实施例中，所述估计步骤包括， 引入所固定的WL�：龋佣兰朴胨潭ǖ腤L�：认嗳莸呐锲�。在一些实施例中，所述估计步骤包括，将迭代滤波器施加到针对每个历元的所述MW数据，并且其中，引入所固定的莫糊度包括以下之一 (i)将所固定的莫糊度作为观测结果置入所述滤波器，(ii)将所固定的莫糊度作为观测结果置入在多个观测更新中的每个更新之后生成的滤波器副本中，以及(iii)将所述MW组合减去所固定的莫糊度，并且将得到的减少的 MW组合置入没有�：茸刺牡诙瞬ㄆ鳎怨兰普攵悦扛鑫佬堑闹辽僖桓鯩W偏差。在一些实施例中，用整数个数的WL周期对至少一个MW偏差进行移位。在一些实施例中，至少一个MW偏差及其相应的WL�：缺灰莆徽鍪腤L周期。在一些实施例中，所述导航消息包括轨道信息。一些实施例提供了用于执行一个或多个上述方法的设备。一些实施例提供了包括被配置为当在计算机处理单元上执行时实施一个或多个上述方法的指令的计算机程序。一些实施例提供了包括此类计算机程序的有形计算机可读介质。
根据下面参考附图所描述的实施例，本发明的所述和其他方面和特征将变得更容易理解，其中图1示出了根据本发明一些实施例的系统的高层级的视图；图2示出了根据本发明一些实施例的系统数据和系统的高层级的视图；图3是根据本发明一些实施例的网络处理器架构的示意图；图4是根据本发明一些实施例的数据校正的示意图；图5是根据本发明一些实施例的观测结果的线性组合的示意图6是一般卡尔曼(Kalman)滤波处理的示意图；图7是根据本发明一些实施例的码分级(code-leveled)的时钟处理器的示意图；图8、图9和图10被删除；图11是根据本发明一些实施例的Melbourne-WUbbena偏差处理流的示意图；图12是根据本发明一些实施例的Melbourne-WUbbena偏差处理流的示意图；图13A示出了根据本发明一些实施例的非差Melbourne-WiAbena偏差处理器的滤波状态；图1 示出了根据本发明一些实施例的单差Melbourne-WUbbena偏差处理器的滤波状态；图14是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图15是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图16是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图17是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图18是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图19A是GNSS站和卫星的观测图；图19B是示出作为顶点的站和卫星以及作为边缘的站-卫星观测结果的抽象图；图19C示出了图19B的最小生成树；图19D示出了具有受约束的边缘的最小生成树；图19E是GNSS站和卫星的非差观测图；图19F是与图19E的观测图对应的滤波图；图19G是GSNN站和卫星的单差观测图；图19H是与图19G的观测图对应的滤波图；图191是比较非差和单差处理中的约束的一组观测结果图；图20是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图21A示出了非差观测图上的生成树；图21B示出了非差观测图上的最小生成树；图21C示出了单差观测图上的生成树；图21D示出了单差观测图上的最小生成树；图22是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图23A是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图2 是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图24A是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图24B是根据本发明一些实施例的Melbourne-WUbbena滤波处理的示意图；图MC是根据本发明一些实施例的Melbourne-WUbbena滤波处理的示意图；图24D是根据本发明一些实施例的Melbourne-WUbbena滤波处理的示意图；图25A是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图；图25B示出了根据本发明一些实施例的移位偏差效应；图25C是根据本发明一些实施例的Melbourne-WUbbena偏差处理器的示意图26A是根据本发明一些实施例的轨道处理器的启动的示意图；图26B是根据本发明一些实施例的轨道处理器的示意图；图26C是根据本发明一些实施例的轨道处理器的轨道映射器的示意图；图26D是根据本发明一些实施例的轨道处理器的轨道映射器的示意图；图27是根据本发明一些实施例的码分级的时钟处理的时序图；图28A是根据本发明一些实施例的高速码分级的卫星时钟处理器的示意图；图28B是根据本发明一些实施例的高速码分级的卫星时钟处理器的示意图；图28C是根据本发明一些实施例的高速码分级的卫星时钟处理器的示意图；图四是根据本发明一些实施例的高速相位分级(phase-leveled)的卫星时钟处理器的示意图；图30A是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图30B是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图30C是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图31是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图32是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图33是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图34是根据本发明一些实施例的高速相位分级的卫星时钟处理器的示意图；图35被删除；图36是根据本发明一些实施例的网络处理器计算机系统的示意图；图37是根据本发明一些实施例的积分GNSS接收机系统的简化示意图；图38是根据本发明一些实施例的使用合成基站数据的GNSS漫游器处理的示意图；图39示出了根据本发明一些实施例的观测时钟预测；图40是根据本发明一些实施例的用于生成合成基站数据的处理的示意图；图41被删除；图42是根据本发明一些实施例的使用合成基站数据的替代性GNSS漫游器处理的示意图；图43是根据本发明一些实施例的使用合成基站数据的替代性GNSS漫游器处理的示意图；图44是根据本发明一些实施例的使用合成基站数据的低时延GNSS漫游器处理的时序图；图45是根据本发明一些实施例的使用合成基站数据的高准确性GNSS漫游器处理的时序图；图46是根据本发明一些实施例的使用合成基站数据的替代性GNSS漫游器处理的示意图；图47示出了与不使用模糊度固定的GNSS漫游器处理相比，根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的性能；图48是根据本发明一些实施例的使用模糊度固定的GNSS漫游器处理的示意图；图49是根据本发明一些实施例的使用模糊度固定的GNSS漫游器处理的示意图50是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图51是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图52是根据本发明一些实施例的使用模糊度固定的GNSS漫游器处理的示意图；图53是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图M是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图55是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图56是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图57是根据本发明一些实施例的使用模糊度固定的GNSS漫游器处理的示意图；图58是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图59是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图60是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图；图61是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图； 禾口图62是根据本发明一些实施例的使用�：裙潭ǖ腉NSS漫游器处理的示意图。
具体实施例方式部分1 系统概览全球导航卫星系统(GNSQ包括GPS、feilileo、Glonass、Compass和其他类似的定位系统。尽管这里给出的例子是针对GPS处理的，然而其原理适用于任何这种定位系统。对实时的定义在本文件中提到了几次术语“实时”。在由下列实施例所涵盖的本发明的范围中，该术语是指只要一个动作所需要的信息可用就有该动作(例如数据被处理、计算结果)。因此，存在特定的时延，并且这依赖于取决于系统元件的不同方面。如下文所述，对于本文件中所涵盖的应用，所需要的信息通常是GNSS数据，和/或GNSS校正。实时运行的网络处理器能够在以下动作之后提供关于来自监控接收机网络的数据的一个历元(epoch)的结果(la)数据被每个监控接收机收集(通常小于1兆秒)；(Ib) 数据被从每个接收机发送给处理中心(通常小于2秒)；(Ic)数据被处理器处理。网络处理器对结果的计算通常占0. 5至5秒之间，这取决于处理器类型和要使用的数据量。通常，在传输延迟上不遵循特定限制(例如3秒)的数据被拒绝或缓存，并且因而没有被立即用于当前的历元更新。这避免了在一个或多个站正在以不可接受的延迟量发送数据的情况下增大系统的时延。实时运行的漫游器接收机能够在接收机收集了数据(通常小于1兆秒)以及以下动作之后，提供关于数据的一个历元的结果Oa)处理中心生成校正数据(参见la、lb、 lc) ； (2b)接收到来自处理中心的校正数据(如果需要的话)(通常小于5秒)数据被处理(通常小于1兆秒)。为了避免或最小化由Qa)和Qb)引发的数据延迟效应，可以使用一种三角相位 (delta phase)方法，使得更新的接收机位置能够在数据被收集之后就被立即与校正数据流进行计算(通常小于1兆秒)。该三角相位方法例如描述在2009年8月18日授权的 U. Vollath的美国专利7 576 690中。图1和图2示出了根据本发明一些实施例的系统100的高层级的视图。地球上分布了全球跟踪网络的参考站，例如参考站105、110、. . . 115。非常精确地知道每个参考站的位置，例如，在小于2厘米以内。每个参考站都配备有天线，并跟踪由该站看到的卫星发送的GNSS信号，所述卫星例如是GNS卫星120、125、. . . 130。GNSS信号具有在两个或多个载波频率的每一个上调制的码。每个参考站获得GNSS数据205，该数据代表了针对在每个历元处看到的每个卫星，至少两个载波的载波-相位(载波)观测结果(observations) 210， 以及在至少两个载波上调制的各自的码的伪距(码)观测结果215。参考站还从卫星信号获得卫星的历书和星历表220。历书包括GNSS的所有卫星的粗略位置，而所谓的广播星历表提供了在特定时间间隔中卫星的时钟误差(大约1.5米)以及卫星位置的更精确的预测 (大约1米)。在参考站收集的GNSS数据经由通信信道135被发送给网络处理器140。如下文详述的，网络处理器140使用来自参考站的GNSS数据以及其他信息来生成包含精确的卫星位置和时钟数据的校正消息。该校正消息被发送给任何数量的GNSS漫游器接收机使用。如图1所示，该校正消息被经由上行链路150和通信卫星155传送以便用于大范围广播；可使用任何其他合适的传输介质，包括但不限于无线电广播或移动电话链路。漫游器160是 GNSS漫游器接收机的实例，其具有用于接收并跟踪在其位置可看到的GNSS卫星的信号的 GNSS天线165，并且可选地具有通信天线170。根据该校正消息的传输频带，其可由漫游器 160经由GNSS天线165或通信天线170来接收。部分2:网络架构图3是示出了根据本发明一些实施例的网络处理器140的处理流300的主要元件的示意图。将来自参考站310的全球网络的GNSS数据作为GNSS数据305不带校正地、或者在通过可选的数据校正器310校正后作为已校正GNSS数据315，提供给四个处理器标准时钟处理器320、Melbourne-Wubbena(MW)偏差处理器325、轨道处理器330以及相位时钟处理器335。数据校正器310可选地分析来自每个参考站的原始GNSS数据305，以检查所接收的观测结果的质量，并且在可能的情况下校正关于周跳(cycle slip)的数据，所述周跳在发生例如每次接收机丢失时钟时，在载波相位观测结果中跳动。商业可用的参考站通常检测周跳并且相应地标记数据。周跳检测和校正技术在例如G. Seeber所著的 SATELLITEGEODESY，第2版(2003年)，第277-281页有所概述。数据校正器310可选地应用其他校正。尽管对于所有处理器而言并非需要所有校正，然而如果被应用到数据，校正确实没有坏处。例如，如下文所述，一些处理器使用码和载波观测结果的线性组合，其中一些未校正误差在形成组合时被抵消。观测结果在每个参考站被逐个历元地获�。⒈淮奔浔昵┑卮透绱砥� 140。对于一些站来说，观测结果延迟到达。这个延迟散布在毫秒到分钟之间。因此，可选的同步器318对在预定时间跨度内的已校正参考站数据的数据进行收集，并将关于每个历元的观测结果作为集合传递给处理器。这使得以合理延迟到达的数据能够被包含在数据的历元中。丽偏差处理器325将未校正GNSS数据305或已校正GNSS数据315作为输入，因为它使用Melbourne-WUbbena线性组合，该组合除了相位和码观测结果的偏差和模糊度之外抵消了所有内容。因此，对于宽巷处理器325，只有接收机和卫星天线校正是重要的。基于该线性组合，计算针对每个卫星的一个MW偏差和针对每个接收机-卫星对的一个宽巷�：�。偏差是平滑的(无噪声)，并且仅表现出一些日以内(sub-daily)的低速变化。宽巷模糊度是恒量，并且只要在相应的卫星-接收机链路上的观测结果中没有发生周跳，就可被使用。因此，时间对于偏差估计而言并不十分关键，并且偏差估计能够例如以15分钟的更新速度来运行。这是有利的，因为计算时间是以站和卫星的数目的三次方来增长的。作为例子，对于具有80个站的合理网络而言，计算时间可以大约是15秒。固定宽巷�：�340 和/或宽巷偏差345的值可选地被用于轨道处理器330和/或相位时钟处理器335，和/或被提供给调度器355。MW偏差处理器325在下文的部分7中进行了详细描述。轨道处理器330的一些实施例基于预测-校正策略。使用精确强制模型(precise force model)并以对卫星参数的未知值的初始推测(初始位置、初始速度和动态强制模型参数)作为开始，每个卫星的轨道通过对卫星的非线性动态系统的积分来预测。包括对未知参数的当前位置的偏导数的敏感度矩阵被同时计算。初始卫星状态的敏感度在针对整个预测的计算的同时被计算。也就是说，起始于对未知参数的预测，求解微分方程组，这将轨道积分到当前时间或将来时间。该预测可被线性化到未知参数的方向中。因此，如果未知参数改变，则偏导数(敏感度)用作对当前卫星状态的改变大小的测量，反之亦然。在一些实施例中，这些偏导数被用于卡尔曼(Kalman)滤波器，以便通过将GNSS观测结果投射到卫星的未知参数来改善初始推测。该精确的初始状态估计用于再次对卫星的动态系统进行积分并确定精确的轨道。不时地执行初始卫星状态到当前历元的时间更新。 在一些实施例中，无电离层�：纫彩强ǘ瞬ㄆ鞯淖刺�。固定的宽巷�：戎�340被用来固定轨道处理器330的无电离层模糊度，以便增强所估计的轨道的准确性。卫星轨道非常平滑，并且能够对若干分钟和小时进行预测。精确的轨道预测结果350可选地被转发给标准时钟处理器320和相位时钟处理器335以及调度器355。超速轨道360，例如国际GNSS服务(IGS)提供的I⑶轨道，可用作对精确轨道预测结果355的替代。I⑶轨道一天更新四次，并且以三小时时延地可用。标准时钟处理器320使用GNSS数据305或已校正GNSS数据315，并使用精确轨道预测结果355或超速轨道365，来计算码分级的卫星时钟360 (也被称为标准卫星时钟)。 码分级的表明时钟与无电离层码观测结果一起使用时是足够的，但是与载波-相位观测结果一起使用时则不是足够的，这是因为码分级的时钟不保持�：鹊恼匦�。由标准时钟处理器320计算的码分级的时钟360表征卫星之间的时钟误差差异。标准时钟处理器 320将广播星历表的时钟误差用作伪观测结果，并将所估计的时钟调整到GPS时间，使得它们可用于计算例如，卫星信号的传输的确切时间。时钟误差变得很快，但是对于与相当吵杂的码测量一起的使用来说，厘米级的准确度是足够的。因此，30秒至60秒的“低速”更新速度是适当的。这是有利的，因为计算时间以站和卫星数目的三次方来增长。标准时钟处理器325还将对流层峰值延迟365确定为估计处理的副产品。对流层峰值延迟和码分级的时钟被发送给相位时钟处理器335。标准时钟处理器320在下文的部分6进行详细描述。相位时钟处理器335可选地连同对流层峰值延迟365和精确轨道350或I⑶轨道 360 一起，使用来自宽巷处理器325的MW偏差345和/或固定宽巷�：�340，以便估计关于每对卫星的窄巷�：群偷ゲ�(singledifference)时钟误差。单差时钟误差和窄巷�：缺蛔楹希曰竦枚杂诿扛鑫佬�(除了参考卫星)的单差相位分级的时钟误差370，其相对于参考卫星而言是单差的。低速码分级的时钟360、对流层峰值延迟365和精确的轨道 350或I⑶轨道360用于估计高速码分级的时钟375。这里，计算工作量与站的数量是线性关系，并且是卫星数量的三次方。快速改变的相位分级的时钟370和码分级的时钟375以例如0. 1秒-0. 2秒的延迟可用。高速相位分级的时钟370和高速码分级的时钟375与MW 偏差340 —起，被发送到调度器355。相位时钟处理器340在下文的部分9进行详细描述。调度器355接收轨道(精确轨道350或I⑶轨道360)、丽偏差340、高速相位分级的时钟370和高速码分级的时钟375。调度器355将这些打包到一起，并将已打包的轨道和时钟以及偏差380转发给消息编码器385，该编码器准备压缩格式的校正消息390，用于向漫游器的传输。到漫游器的传输在卫星链路上占用例如大约10秒-20秒，但也可以使用移动电话或直接互联网连接或其它合适的通信链路来完成。部分3 观测(observation)数据校正器图4是根据本发明一些实施例的数据校正的示意图。可选的观测校正器310对在参考站收集的GNSS信号进行校正，该校正针对由于作用于地球的离心力、回转力和重力而造成的站移位，站的天线相位中心相对于站的天线安装点的位置，站的天线相位中心相对于卫星轨道给出的卫星质量中心的位置，以及取决于站天线和卫星天线的对准方式的这些相位中心的变化。站位移的主要原因是高达500毫米的固体潮(earth tide)、高达100毫米的海洋潮汐负荷，以及高达10毫米的极点潮。所有这些都取决于站位于哪里。可在McCarthy， D. D.，Petit, G.(编)IERS Conventions (2003)，IERS Technical Note No. 32 以及其中引用的参考文献中找到更多描述。由作用于地球的松散体上的天体(主要是月亮)力导致的海洋潮汐还会引起大陆板块被提升和降低。该熟知的效应表现为参考站位置的重复变化。可选地对固体地球潮进行计算以便用于网络处理以及漫游器处理，因为该效应不应当被忽略并且计算工作量较小。第二大效应是大陆的板块变形，这归因于与潮汐一起随时间变化的海洋负荷。用于快速计算站随时间的位移的海洋潮汐负荷参数取决于站的位置。用来得出这些参数的计算工作量相当大。对于给定的位置，可使用任何可通过Onsala空间天文台海洋(Onsala Space Observatory Ocean, http //www. oso. chalmers. se/ ~ loading/, Chalmers Onsala Space0bservatory，2009)提供的在线海洋-潮汐-负荷(ocean-tide-loading)服务获得的熟知模型，对这些参数进行计算。较低准确度的参数，例如来自预先计算的网格的内插， 对于这里讨论的应用而言是足够的。这里提到的最小效应归因于极点潮。这个位移归因于地球的两极运动产生的离心和回转效应引起的板块的升起。地球定向参数用于该计算。这些会被有规律地更新在国际地球自转禾口参考系统月艮务(International Earth Rotation & Reference System Service), http://hpiers. obspm. fr/，巴黎天文台，2009，并且不容易计算。这个较小的效应因而在漫游器处理中可选地被忽略。绝对校正的天线模型用于计算接收机和卫星天线相位中心的偏移和变化。可在 J. Kouba 所著的 “A Guide to Using International GPS Service(IGS) Products，，(Geoodetic Survey Division Natural Resources，力口拿大，2003 年 2 月)找到介绍。通过IGS收集的校正数据可从http://igscb. jpl. nasa. gov/上2009年的antex文件中获得；卫星天线偏移信息例如可在IGS绝对天线文件(IGS absolute antenna file) igs05. atx 中找到。另一个效应是天线相位缠绕(wind-up)。如果接收机天线相对于发送器天线移动，则所记录的数据就显示出相位移动。如果该效应被忽略，则卫星围绕发送轴的360度转体(full turn)将导致在接收机处检测到的载波-相位中一个周期的误差。由于卫星相对于接收机的定向多数时候是已知的，所以该效应可被建模，如mi J. Τ. , Hajj G.A.， Bertiger W. I.以及 LichtenS. Μ·的"Effects of antenna orientation on GPS carrier phase” (MANUSCRIPTA GE0DAETICA，第 18 卷 91-98 页(1993))中所介绍的。站和卫星的相对移动主要原因是绕轨道运行的卫星。如果卫星是遮挡的(这意味着卫星的轨道横跨地球的阴影)，则可能有卫星围绕其发送轴的额外旋转。例如，GPS Block IIA卫星具有正午旋转以及交叉运行(crossing maneuver)的阴影，而GPS Block UR卫星则具有正午旋转和午夜旋转。如果太阳、地球和卫星几乎在同一直线上，则难以计算旋转运行的方向，并且不正确的方向将导致一个周期的载波-相位中的误差。卫星的偏航飞行姿态对相位缠绕和卫星天线校正有影响。可在Kouba，J.的“Asimplified yaw-attitude model for eclipsing GPS satellites，，(GPSS0LUTI0NS (2008))、和 Barlever，Y. E.的"A new model for GPS yaw attitude"(JOURNAL OF GEODESY，第 70 卷 714-723 页(1996))中找到更详细的介绍。在仅使用相位观测结果的情况下，未建模的卫星旋转运行的效应不能与卫星时钟分离开。因此，在相位时钟误差估计中，旋转运行的效应被包含于所估计的卫星时钟误差中。如果漫游器使用那些时钟，则其必须也不对卫星旋转运行进行校正。需要太阳位置以便对卫星体固定的坐标框架进行计算，因为χ轴是通过卫星位置和太阳位置的向量积来定义的。该坐标系统用于计算偏航飞行姿态、卫星的天线校正(偏移和变化，映射到正弦视角中)以及相位缠绕校正。对于固体地球潮来说，也需要月亮的位置。可在例如 Seidelmann, P. K.(编)"Explanatory Supplement to the Astronomical Almanac”(University Science Books, U. S. (1992))中找到如何计算太阳和月亮的位置。还可应用进一步的校正，尽管这些对于市场需要的定位准确度级别而言仅仅是较小的益处。不需要在可选的数据校正器310中考虑作为校正关于相对效应、电离层和对流层延迟的附加效应。相对校正通常被应用于卫星时钟。归因于电离层的主要第一顺序效应通过使用GNSS观测结果的无电离层组合而被消除，并且归因于对流层的效应在一些实施例中被部分地建模并被部分地估计。部分4:形成线性组合4.1基本建模等式对于接收机i和卫星j之间调制类型为m、频带k上的码巧二和载波相位Φ/-观测结果而言，假设以下观测模型将观测结果关联到特定物理量，
权利要求
1.一种处理从多个接收机处对多个历元上的多个卫星的GNSS信号的码观测结果和载波-相位观测结果得到的一组GNSS信号数据的方法，所述GNSS信号具有至少两个载波频率，包括形成在每个历元处针对每个接收机-卫星对的MW(MeIbourne-WiAbena)组合，以获得针对每个历元的MW数据组，以及从针对每个历元的棚数据组估计针对每个卫星的棚偏差和一组WL (宽巷)模糊度， 所述MW偏差可从一个历元到另一个历元地变化，每个WL�：榷杂τ诮邮栈�-卫星链路和卫星-接收机-卫星链路中的一个，其中，针对每个卫星的所述MW偏差被建模成以下之一 (i)单个估计参数和(ii)所估计的偏移加上具有所估计的振幅的谐波变化。
2.根据权利要求1所述的方法，还包括对所述GNSS信号数据施加校正。
3.根据权利要求1所述的方法，还包括在估计针对每个卫星的MW偏差之前，平滑GNSS 信号数据的至少一个线性组合。
4.根据权利要求1所述的方法，还包括施加至少一个MW偏差约束。
5.根据权利要求1所述的方法，还包括施加至少一个整数WL�：仍际�。
6.根据权利要求1所述的方法，还包括在观测图和滤波图之一上使用生成树(ST)以约束所述莫糊度。
7.根据权利要求1所述的方法，还包括在观测图和滤波图之一上使用最小生成树 (MST)以约束所述莫糊度。
8.根据权利要求7所述的方法，其中所述最小生成树基于边权重，每个边权重是从接收机-卫星几何得到的。
9.根据权利要求8所述的方法，其中，所述边权重是关于以下之一来限定的(i)从接收机到卫星的几何距离，(ii)卫星仰角角度，以及(iii)从卫星到接收机到卫星的几何距离，以及(iv)仰角的组合，在所述仰角下，在站处可见单差组合中的两个卫星。
10.根据权利要求7所述的方法，其中，所述最小生成树基于边权重，每个边权重都基于W4莫糊度信息，并且是关于以下之一来限定的(i)WL�：鹊秸牟睿�(ii)m4莫糊度方差，以及(iii)⑴和( )的组合。
11.根据权利要求ι所述的方法，还包括将所述莫糊度中的至少一个固定为整数值。
12.根据权利要求ι所述的方法，还包括确定WL整数�：戎档亩喔龊蜓∽椋纬伤龆喔龊蜓∽榈募尤ㄗ楹希约敖鯳L�：戎械闹辽僖桓龉潭ㄎ∽运黾尤ㄗ楹系闹�。
13.根据权利要求11所述的方法，其中，所述估计的步骤包括引入所固定的莫糊度， 从而估计与所固定的WL�：认嗳莸腗W偏差。
14.根据权利要求13所述的方法，其中，所述估计的步骤包括，将迭代滤波器施加到针对每个历元的所述丽数据，并且其中，引入所固定的莫糊度包括以下之一⑴将所固定的莫糊度作为观测结果置入所述滤波器，(ii)将所固定的莫糊度作为观测结果置入在多个观测更新中的每个更新之后生成的滤波器的副本中，以及(iii)将所述MW组合减去所固定的莫糊度，并且将得到的减少的MW组合置入没有模糊度状态的第二滤波器，以估计针对每个卫星的至少一个MW偏差。
15.根据权利要求1所述的方法，还包括用整数个数的WL周期对至少一个MW偏差进行移位。
16.根据权利要求1所述的方法，还包括用整数个数的WL周期对至少一个MW偏差及其相应的WL模糊度进行移位。
17.根据权利要求1所述的方法，其中，所述导航消息包括轨道信息。
18.一种用于执行根据权利要求1-17之一所述的方法的设备。
19.一种处理从多个接收机处对多个历元上的多个卫星的GNSS信号的码观测结果和载波-相位观测结果得到的一组GNSS信号数据的设备，所述GNSS信号具有至少两个载波频率，包括适于形成在每个历元处针对每个接收机-卫星对的MW(Melbourne-WiAbena)组合，以获得针对每个历元的MW数据组的元件，以及适于从针对每个历元的棚数据组估计针对每个卫星的棚偏差和一组WL (宽巷)�：鹊脑雠锲羁纱右桓隼搅硪桓隼乇浠扛鯳L�：榷杂τ诮邮栈�-卫星链路和卫星-接收机-卫星链路中的一个，其中，针对每个卫星的所述MW偏差被建模成以下之一 (i)单个估计参数和(ii)所估计的偏移加上具有所估计的振幅的谐波变化。
20.根据权利要求19所述的设备，还包括适于对所述GNSS信号数据施加校正的元件。
21.根据权利要求19所述的设备，还包括适于在估计针对每个卫星的MW偏差之前，平滑GNSS信号数据的至少一个线性组合的元件。
22.根据权利要求19所述的设备，还包括适于施加至少一个MW偏差约束的元件。
23.根据权利要求19所述的设备，还包括适于施加至少一个整数莫糊度约束的元件。
24.根据权利要求19所述的设备，还包括适于在观测图和滤波图之一上使用生成树 (ST)以约束所述莫糊度的元件。
25.根据权利要求19所述的设备，还包括在观测图和滤波图之一上使用最小生成树 (MST)以约束所述莫糊度的元件。
26.根据权利要求25所述的设备，其中所述最小生成树基于边权重，每个边权重是从接收机-卫星几何得到的。
27.根据权利要求沈所述的设备，其中，所述边权重是关于以下之一来限定的(i)从接收机到卫星的几何距离，(ii)卫星仰角角度，以及(iii)从卫星到接收机到卫星的几何距离，以及(iv)仰角的组合，在所述仰角下，在站处可见单差组合中的两个卫星。
28.根据权利要求25所述的设备，其中，所述最小生成树基于边权重，每个边权重都基于莫糊度信息，并且是关于以下之一来限定的(i)W4莫糊度到整数的差，(ii)W4莫糊度方差，以及(iii)⑴和( )的组合。
29.根据权利要求19所述的设备，还包括适于将所述莫糊度中的至少一个固定为整数值的元件。
30.根据权利要求19所述的设备，还包括适于确定WL整数�：戎档亩喔龊蜓∽椋纬伤龆喔龊蜓∽榈募尤ㄗ楹希约敖鯳L�：戎械闹辽僖桓龉潭ㄎ∽运黾尤ㄗ楹系闹档脑�。
31.根据权利要求四所述的设备，其中，适于估计的元件还适于引入所固定的莫糊度，从而估计与所固定的WL�：认嗳莸睦銎�。
32.根据权利要求31所述的设备，其中，适于估计的元件还适于将迭代滤波器施加到针对每个历元的所述MW数据，并将所固定的莫糊度作为以下之一来引入(i)所固定的莫糊度作为滤波器中的观测结果，(ii)所固定的莫糊度作为在多个观测更新中的每个更新之后生成的滤波器的副本中的观测结果，以及(iii)所述MW组合被减去所固定的 WL�：炔⒋τ诿挥心：茸刺牡诙瞬ㄆ髦幸员愎兰普攵悦扛鑫佬堑闹辽僖桓鯩W偏差。
33.根据权利要求19所述的设备，还包括适于用整数个数的WL周期对至少一个MW偏差进行移位的元件。
34.根据权利要求19所述的设备，还包括适于用整数个数的WL周期对至少一个MW偏差及其相应的WL�：冉幸莆坏脑�。
35.根据权利要求19所述的设备，其中，所述导航消息包括轨道信息。
36.一种计算机程序，其包括被配置为当在计算机处理单元上执行时实施根据权利要求1所述的方法的指令。
37.一种计算机可读介质，包括根据权利要求36所述的计算机程序。
全文摘要
描述了用于处理从多个接收机处对多个历元上的多个卫星的GNSS信号的码观测结果和载波-相位观测结果得到的一组GNSS信号数据的方法和设备，所述GNSS信号具有至少两个载波频率，包括形成在每个历元处针对每个接收机-卫星对的MW(Melbourne-Wübbena)组合以获得针对每个历元的MW数据组，以及从针对每个历元的MW数据组估计针对每个卫星的MW偏差和一组WL(宽巷)�：龋鯩W偏差可从一个历元到另一个历元地变化，每个WL模糊度对应于接收机-卫星链路和卫星-接收机-卫星链路中的一个，其中，针对每个卫星的所述MW偏差被建模成以下之一(i)单个估计参数和(ii)所估计的偏移加上具有所估计的振幅的谐波变化。
文档编号G01S19/37GK102576081SQ201080041634
公开日2012年7月11日 申请日期2010年9月19日 优先权日2009年9月19日
发明者S·塞格 申请人:天宝导航有限公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
