一種長基線的偽距差分定位方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種偽距差分定位方法，屬于衛(wèi)星導(dǎo)航領(lǐng)域。
【背景技術(shù)】
[0002] 隨著全球衛(wèi)星導(dǎo)航系統(tǒng)（GNSS)的廣泛應(yīng)用，越來越多的電子設(shè)備需要通過衛(wèi)星 來確定其精確位置。GNSS單點定位精度已經(jīng)遠遠不能滿足車輛控制、航空管制、飛機編隊 等方面的需求，急需開發(fā)和研究高精度的定位設(shè)備和定位方法。衛(wèi)星接收機測量信息中包 含著各種誤差：衛(wèi)星時鐘誤差、衛(wèi)星星歷誤差、對流層延時誤差、電離層延時誤差、多路徑誤 差、接收機噪聲等。減少衛(wèi)星接收機測量信息中的測量誤差是提高定位精度的有效方法之 一，差分GNSS是一種應(yīng)用廣泛、能夠有效消除各種測量誤差的方法。差分GNSS的基本工作 原理主要依據(jù)處在同一區(qū)域內(nèi)的不同接收機，它們的測量值中所包含的衛(wèi)星時鐘誤差、衛(wèi) 星星歷誤差、對流層延時誤差、電離層延時誤差這四種誤差近似相等或高度相關(guān)。差分GNSS 定位根據(jù)GNSS基準(zhǔn)站發(fā)送的不同信息可分為：位置差分、偽距差分、載波相位差分三類。偽 距差分是當(dāng)前應(yīng)用最廣泛的一類差分定位模式。傳統(tǒng)的偽距差分GNSS方法假設(shè)用戶端、基 準(zhǔn)站處對同一顆衛(wèi)星的觀測向量是相互平行的，當(dāng)基線小于l〇Km時，這一假設(shè)帶來的模型 誤差很小，但隨著基線長度逐漸增大時模型本身引起的定位誤差也逐漸增大。

【發(fā)明內(nèi)容】

[0003]為了克服現(xiàn)有技術(shù)的不足，本發(fā)明提供一種長基線的偽距差分定位方法，可以滿 足高精度定位的需求。
[0004] 本發(fā)明解決其技術(shù)問題所采用的技術(shù)方案是：由若干顆衛(wèi)星s(1)、用戶端U和基準(zhǔn) 站r組成本發(fā)明的衛(wèi)導(dǎo)偽距差分模型。首先使用基準(zhǔn)站、用戶端與衛(wèi)星S(1)的偽距觀測方程 構(gòu)建相對于同一顆衛(wèi)星的偽距單差方程，其次利用不同衛(wèi)星的偽距單差方程構(gòu)建包含多顆 衛(wèi)星的非線性的偽距雙差方程組，最后采用牛頓迭代最小二乘方法求解出用戶端位置。具 體步驟如下：
[0005] 第一步，假定t時刻某顆衛(wèi)星S(1)在地心直角坐標(biāo)系的坐標(biāo)為（X;，y;，zj， 基準(zhǔn)站r處接收機的地心坐標(biāo)為（Xyyyzj，基準(zhǔn)站r處接收機與衛(wèi)星S(1)之間的偽
,其中，f是基準(zhǔn)站r與衛(wèi)星S(1)的
c代表光速，δ~代表基準(zhǔn) 站r處接收機的時鐘鐘差，St(1)代表衛(wèi)星S(1)的時鐘鐘差，代表衛(wèi)星S(1)信號到基 準(zhǔn)站r的電離層延時，代表衛(wèi)星S(1)信號到基準(zhǔn)站r的對流層延時，代表基準(zhǔn) 站r到衛(wèi)星S(1)偽距測量值的隨機噪聲量；用戶端u的接收機與衛(wèi)星S(1)之間的偽距
，其中，3|3:是用戶端u與衛(wèi)星S(1)的幾何距離，
stu代表用戶端u處接收機的時鐘鐘差，if代表衛(wèi)星S(1)信號到用戶端u的電離層延時，if代表衛(wèi)星S(1)信號到用戶端u的對流層延 時，^代用戶端u到衛(wèi)星S(1)偽距測量值的隨機噪聲量；
[0006] 第二步，計算基準(zhǔn)站、用戶端與衛(wèi)星S(1)的單差偽距，即用戶端相對于衛(wèi)星S(1)的偽
[0007] 同理，可得基準(zhǔn)站、用戶端與衛(wèi)星S(])的偽距單差方程式如下：
[0009] 第三步，將用戶端與基準(zhǔn)站都能收到信號的衛(wèi)星稱為共視衛(wèi)星，在同一時刻共
[0010] 第四步，將用戶端、基準(zhǔn)站的η顆共視衛(wèi)星中仰角最大的衛(wèi)星作為主星S(k)，得到
[0011] 第五步，利用牛頓迭代最小二乘法求解雙差偽距觀測方程組，包括以下步驟：
[0012] a)設(shè)置迭代初始值m= 0,用戶端位置初始值XQ=(Xpypzj;
[0013] b)得到線性化的雙差偽距觀測方程式A·ΔX=B，
r
[0017] c)利用最小二乘法求解線性化的雙差偽距觀測方程式，得到ΔX= (ATA)jATB，其 中Ατ代表A的轉(zhuǎn)置，（ATA) 1代表ATA的逆矩陣；
[0018] d)更新第m+1次的迭代值Xm+1 =Xm+ΔX;
[0019] e)如果第m+1次與第m次迭代值的誤差小于預(yù)先設(shè)定的門限值，則將X"+1作為用 戶端的定位位置結(jié)果，步驟結(jié)束；否則m值增加1，返回步驟b)。
[0020] 本發(fā)明的有益效果是：改進傳統(tǒng)的偽距差分定位模型假設(shè)，即在用戶端、基準(zhǔn)站處 對同一顆衛(wèi)星的觀測向量是相互平行的，有效地解決了傳統(tǒng)平行假設(shè)在基線長度逐漸增大 時模型本身引起的定位誤差也逐漸增大的問題；在本發(fā)明的定位方法下，消除了不同長度 基線下偽距差分定位模型本身的誤差，使其不隨著基線長度的增加而增加，可以維持在很 小的誤差范圍內(nèi)，完全可以滿足高精度定位的需求。
[0021] 此外本發(fā)明也可以應(yīng)用于飛機編隊的相對定位。例如一個飛機編隊有一架長機和 多架架僚機，可以把長機看成用戶端，把每一架僚機看成是移動的基準(zhǔn)站，那么每一架僚機 利用自己的位置、僚機機載接收機的偽距測量值和長機機載接收機的偽距測量值，可以獲 得僚機相對于長機的位置，從而完成飛機編隊的相對定位。
【附圖說明】
[0022] 圖1是本發(fā)明的衛(wèi)導(dǎo)偽距差分模型示意圖；
[0023] 圖2是本發(fā)明方法的流程框架示意圖；
[0024] 圖3是傳統(tǒng)偽距差分模型示意圖；
[0025] 圖4是傳統(tǒng)方法定位誤差示意圖；
[0026] 圖5是本發(fā)明定位誤差示意圖。
【具體實施方式】
[0027] 下面結(jié)合附圖和實施例對本發(fā)明進一步說明，本發(fā)明包括但不僅限于下述實施 例。
[0028] 圖1是本發(fā)明的衛(wèi)導(dǎo)偽距差分模型示意圖，它由若干顆衛(wèi)星S(1)，用戶端u，基準(zhǔn)站 r組成。已知基準(zhǔn)站位置、基準(zhǔn)站與若干顆衛(wèi)星S(1)的偽距、用戶端與若干顆衛(wèi)星S(1)的偽 距，并利用這些位置和偽距測量信息求解用戶端位置。圖2是本發(fā)明方法的流程框架圖。
[0029] 第一步：建立偽距觀測方程式
[0030] 假定t時刻某顆衛(wèi)星S(1)在地心直角坐標(biāo)系的坐標(biāo)為（xpyi，Zl)，基準(zhǔn)站r處接收 機的地心坐標(biāo)為（\，I，zj，基準(zhǔn)站r處接收機與衛(wèi)星S(1)之間的偽距測量值為pf，則基準(zhǔn) 站相對于衛(wèi)星的偽距觀測方程式可表示為：
[0032]其中，是基準(zhǔn)站r與衛(wèi)星S(1)的幾何距離，即
[0034] 方程式⑴中圓括號內(nèi)的上標(biāo)（如i)代表衛(wèi)星編號，而下標(biāo)（如r和u分別代表 基準(zhǔn)站接收機和用戶端接收機），其中代表接收機時鐘鐘差，St(1)代表衛(wèi)星s(1)時鐘 鐘差，if代表衛(wèi)星s(1)信號到基準(zhǔn)站的電離層延時，Tf代表衛(wèi)星s(1)信號到基準(zhǔn)站的對流 層延時^代表基準(zhǔn)站到衛(wèi)星s(1)偽距測量的隨機噪聲量，C代表光速。方程式（1)表示偽 距測量值包含基準(zhǔn)站r與衛(wèi)星S(1)的幾何距離、接收機和衛(wèi)星時鐘鐘差誤差、電離層延時和 對流層延時誤差、偽距測量隨機噪聲量。
[0035] 同理，可得用戶端u的接收機與衛(wèi)星S(1)之間的偽距觀測方程式pf如下：
[0037] 第二步：建立偽距單差觀測方程式。
[0038] 由公式（1