用以估計mw偏差的gnss信號處理的制作方法

專利名稱：用以估計mw偏差的gnss信號處理的制作方法
技術領域：
本發明涉及GNSS信號處理，并且具體涉及包含精確衛星數據的GNSS信號處理。

發明內容
描述了用于處理從多個接收機處對多個歷元上的多個衛星的GNSS信號的碼觀測結果和載波-相位觀測結果得到的一組GNSS信號數據的方法和設備，所述GNSS 信號具有至少兩個載波頻率，包括形成在每個歷元處針對每個接收機-衛星對的麗(Melbourne-WUbbena)組合，以獲得針對每個歷元的麗數據組，以及從針對每個歷元的MW數據組估計針對每個衛星的MW偏差和一組WL (寬巷)模糊度，所述MW偏差可從一個歷元到另一個歷元地變化，每個WL模糊度對應于接收機-衛星鏈路和衛星-接收機-衛星鏈路中的一個，其中，針對每個衛星的所述MW偏差被建模成以下之一 (i)單個估計參數和 ( )所估計的偏移加上具有所估計的振幅的諧波變化。在一些實施例中，根據1的方法，校正被施加到所述GNSS信號數據。在一些實施例中，在估計針對每個衛星的MW偏差之前，GNSS信號數據的至少一個線性組合被平滑。在一些實施例中，施加至少一個MW偏差約束。在一些實施例中，施加至少一個整數WL模糊度約束。在一些實施例中，在觀測圖和濾波圖之一上使用生成樹(ST)以約束所述莫糊度。 在一些實施例中，在觀測圖和濾波圖之一上使用最小生成樹(MST)以約束所述莫糊度。在一些實施例中，所述最小生成樹基于邊權重，每個邊權重是從接收機-衛星幾何得到的。在一些實施例中，所述邊權重是關于以下之一來限定的(i)從接收機到衛星的幾何距離，(ii)衛星仰角角度，以及(iii)從-衛星到接收機到衛星的幾何距離，以及(iv) 仰角的組合，在所述仰角下，在站處可見單差組合中的兩個衛星。在一些實施例中，所述最小生成樹基于邊權重，每個邊權重都基于WL模糊度信息，并且是關于以下之一來限定的 (i)m4莫糊度到整數的差，(ii)m4莫糊度方差，以及(iii)⑴和(ii)的組合。在一些實施例中，所述莫糊度中的至少一個被固定為整數值。在一些實施例中，確定WL整數模糊度值的多個候選組，形成所述多個候選組的加權組合，以及將所述WL 模糊度中的至少一個固定為取自所述加權組合的值。在一些實施例中，所述估計步驟包括， 引入所固定的WL模糊度，從而估計與所固定的WL模糊度相容的棚偏差。在一些實施例中，所述估計步驟包括，將迭代濾波器施加到針對每個歷元的所述MW數據，并且其中，引入所固定的莫糊度包括以下之一 (i)將所固定的莫糊度作為觀測結果置入所述濾波器，(ii)將所固定的莫糊度作為觀測結果置入在多個觀測更新中的每個更新之后生成的濾波器副本中，以及(iii)將所述MW組合減去所固定的莫糊度，并且將得到的減少的 MW組合置入沒有模糊度狀態的第二濾波器，以估計針對每個衛星的至少一個MW偏差。在一些實施例中，用整數個數的WL周期對至少一個MW偏差進行移位。在一些實施例中，至少一個MW偏差及其相應的WL模糊度被移位整數個數的WL周期。在一些實施例中，所述導航消息包括軌道信息。一些實施例提供了用于執行一個或多個上述方法的設備。一些實施例提供了包括被配置為當在計算機處理單元上執行時實施一個或多個上述方法的指令的計算機程序。一些實施例提供了包括此類計算機程序的有形計算機可讀介質。


根據下面參考附圖所描述的實施例，本發明的所述和其他方面和特征將變得更容易理解，其中圖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漫游器處理相比，根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的性能；圖48是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖49是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖50是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖51是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖52是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖53是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖M是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖55是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖56是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖57是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖58是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖59是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖60是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖；圖61是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖； 禾口圖62是根據本發明一些實施例的使用模糊度固定的GNSS漫游器處理的示意圖。
具體實施例方式部分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整數模糊度值的多個候選組，形成所述多個候選組的加權組合，以及將所述WL模糊度中的至少一個固定為取自所述加權組合的值。
13.根據權利要求11所述的方法，其中，所述估計的步驟包括引入所固定的莫糊度， 從而估計與所固定的WL模糊度相容的MW偏差。
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 (寬巷)模糊度的元件，所述棚偏差可從一個歷元到另一個歷元地變化，每個WL模糊度對應于接收機-衛星鏈路和衛星-接收機-衛星鏈路中的一個，其中，針對每個衛星的所述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整數模糊度值的多個候選組，形成所述多個候選組的加權組合，以及將所述WL模糊度中的至少一個固定為取自所述加權組合的值的元件。
31.根據權利要求四所述的設備，其中，適于估計的元件還適于引入所固定的莫糊度，從而估計與所固定的WL模糊度相容的麗偏差。
32.根據權利要求31所述的設備，其中，適于估計的元件還適于將迭代濾波器施加到針對每個歷元的所述MW數據，并將所固定的莫糊度作為以下之一來引入(i)所固定的莫糊度作為濾波器中的觀測結果，(ii)所固定的莫糊度作為在多個觀測更新中的每個更新之后生成的濾波器的副本中的觀測結果，以及(iii)所述MW組合被減去所固定的 WL模糊度并處于沒有模糊度狀態的第二濾波器中以便估計針對每個衛星的至少一個MW偏差。
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(寬巷)模糊度，所述MW偏差可從一個歷元到另一個歷元地變化，每個WL模糊度對應于接收機-衛星鏈路和衛星-接收機-衛星鏈路中的一個，其中，針對每個衛星的所述MW偏差被建模成以下之一(i)單個估計參數和(ii)所估計的偏移加上具有所估計的振幅的諧波變化。
文檔編號G01S19/37GK102576081SQ201080041634
公開日2012年7月11日 申請日期2010年9月19日 優先權日2009年9月19日
發明者S·塞格 申請人:天寶導航有限公司