[發明專利]一種基于三頻模糊度解算的整周模糊度計算方法

在審

• 申請號: 202110408310.7
• 申請日: 2021-04-16
• 公開（公告）號: CN112987059A
• 公開（公告）日: 2021-06-18
• 發明（設計）人: 梁曉東;匡宇龍;周延棟;雷孟飛;熊用;姜水軍;曾創
• 申請（專利權）人: 湖南聯智科技股份有限公司
• 主分類號: G01S19/44
• 分類號: G01S19/44;G06F17/11
• 代理公司: 長沙七源專利代理事務所(普通合伙) 43214
• 代理人: 張勇;周曉艷

1.一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，包括以下步驟：

步驟A：通過對原始觀測量的多頻測量值進行線性組合，得到超寬巷組合觀測量、寬巷組合觀測量和窄巷組合觀測量；

步驟B：通過三頻模糊度解算算法對各組合觀測量進行逐級模糊度解算，得到超寬巷模糊度整數值、寬巷模糊度整數值和窄巷組合觀測量偽距的平滑值；

步驟C：將步驟B中得到的窄巷組合觀測量偽距的平滑值結合LAMBDA算法進行搜索，確定窄巷模糊度整數值；

步驟D：結合超寬巷模糊度整數值、寬巷模糊度整數值和通過LAMBDA算法確定的窄巷模糊度整數值，得到原始整周模糊度。

2.根據權利要求1所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟A中，每個歷元的原始觀測量包括雙差偽距觀測量和以米為單位的雙差載波相位觀測量，每個歷元中第i個頻點對應的雙差偽距觀測量和以米為單位的雙差載波相位觀測量分別通過表達式1）和表達式2）計算：

1）；

其中，代表衛星與接收機之間的雙差空間矢量，代表第i個頻點的雙差電離層誤差項，代表雙差對流層誤差項，代表第i個頻點的偽距觀測量的隨機觀測誤差；

2）；

其中，為第i個頻點對應的波長，為第i個頻點的原始整周模糊度，代表第i個頻點的載波相位觀測量的隨機觀測誤差。

3.根據權利要求2所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟A中，超寬巷組合觀測量、寬巷組合觀測量和窄巷組合觀測量均包括以周期為單位的雙差載波相位觀測量和組合波長，和分別通過表達式3）和表達式4）計算：

3）；

4）；

其中，、、分別代表三個頻點的線性組合系數，、、分別代表三個頻點對應的波長，g代表設備群時延，為電離層延遲放大系數，為組合觀測量的整周模糊度，代表組合觀測量的隨機觀測誤差；下標Z取EWL、WL和NL三者中的一個，EWL、WL和NL分別代表與超寬巷組合觀測量、寬巷組合觀測量和窄巷組合觀測量相關的變量。

4.根據權利要求3所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟A中，通過表達式5）計算電離層延遲放大系數：

5）。

5.根據權利要求4所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟B中，進行逐級模糊度解算時先將超寬巷組合觀測量的模糊度固定后得到超寬巷組合觀測量偽距的平滑值，再將超寬巷組合觀測量偽距的平滑值用于寬巷組合觀測量的模糊度固定；將寬巷組合觀測量的模糊度固定后得到寬巷組合觀測量偽距的平滑值，最后將寬巷組合觀測量偽距的平滑值用于求解窄巷組合觀測量偽距的平滑值。

6.根據權利要求5所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟B中通過表達式7）計算模糊度整數值：

7）；

其中，，代表的隨機觀測誤差。

7.根據權利要求6所述的一種基于三頻模糊度解算的整周模糊度計算方法，其特征在于，所述步驟B中通過表達式8）計算超寬巷組合觀測量、寬巷組合觀測量和窄巷組合觀測量偽距的平滑值：

8）。