引言

    隨著科技的發(fā)展，因為應用的廣泛GPS是被我們所熟知的一項技術(shù)。GPS就是通過接受衛(wèi)星信號，進行定位或者導航的終端。而接受信號就必須用到天線。 GPS衛(wèi)星信號分為L1和L2，頻率分別為1575.42MHZ和1228MHZ，其中L1為開放的民用信號，信號為圓形極化。信號強度為166-DBM左右，屬于比較弱的信號。 這些特點決定了要為GPS信號的接受準備專門的天線。

　　GPS衛(wèi)星定位解算是根據(jù)偽距、偽距增量等測量值，計算接收機的位置P、速度V和時間T等信息的過程。目前GPS實時定位解算中為常用的兩種方法為迭代二乘算法（ILS）和擴展卡爾曼濾波（EKF）。為了準確計算接收機的三維位置以及時間未知數(shù)的值，解算過程需要至少4顆衛(wèi)星的測量信息。但是當GPS信號出現(xiàn)遮擋時，接收機只能接收到3顆衛(wèi)星的測量信息，解算方程就不夠4個，ILS不再適用。引入EKF，利用隨著時間推移的多組數(shù)據(jù)進行實時定位解算，但定位也很難滿足用戶的需求。

　　為了解決上述問題，本文提出了一種改進的EKF算法。利用在垂直地面方向上的位置變化緩慢這一運動特性，建立了改進EKF算法的系統(tǒng)模型，并通過理論分析得到了濾波器參數(shù)，利用真實的衛(wèi)星數(shù)據(jù)進行驗證。需要特別指出的是，由于本文提出的改進EKF算法利用的是在垂直地面方向上的位置變化緩慢的特征，故而該算法的適用場合為車載等地面用戶的定位解算，不適合在垂直地面方向上高速運動的情形。

　　1 定位解算的系統(tǒng)模型

　　衛(wèi)星定位解算的系統(tǒng)模型包括狀態(tài)模型和觀測模型兩部分。令向量yt、xt分別表示系統(tǒng)模型的測量值和系統(tǒng)狀態(tài)參量：

　　1.1 測量模型

　　系統(tǒng)的測量模型描述了系統(tǒng)測量值與系統(tǒng)狀態(tài)參量之間的關(guān)系。偽距與系統(tǒng)狀態(tài)參量的關(guān)系可表示為：

　　1.2狀態(tài)模型

　　系統(tǒng)的狀態(tài)模型描述了系統(tǒng)狀態(tài)參量的時間更新過程。更新過程的表達式為：

　　式（9）中的T為采樣時間間隔。

　　式（8）中的wt表示系統(tǒng)狀態(tài)轉(zhuǎn)移的噪聲模型：

　　2 改進的EKF

　　根據(jù)第1節(jié)中建立的系統(tǒng)模型，系統(tǒng)測量值yt為2參數(shù)，系統(tǒng)狀態(tài)參量xt為8參數(shù)，為了得到準確解，故而需要至少4顆有效衛(wèi)星的測量值。為了能夠滿足用戶需求的定位，本文提出利用一般車輛行人等用戶在垂直地面方向上的位置變化緩慢這一運動特性，增加如下方程：

　　式（17）中，R1為從用戶在ECEF坐標系下的狀態(tài)參量xt到用戶在NED坐標系下的狀態(tài)參量mt的轉(zhuǎn)換矩陣，利用坐標系旋轉(zhuǎn)的原理可以求得。A2為用戶在NED坐標系下的狀態(tài)參量mt的一步轉(zhuǎn)移矩陣。它們的表達式如下：

　　綜上所述，得到改進的EKF算法的計算過程如下：

　　根據(jù)理論分析可以得出，改進的EKF相對于普通的EKF來說，增加了一項先驗信息，故而使得只有3顆有效衛(wèi)星時的定位解算有了4個解算方程，因此可以相對準確地解算出4個未知數(shù)的值；而在有效衛(wèi)星數(shù)不少于4顆時，增加的方程（15）仍然滿足，故而也不會影響定位。

　　3 算法仿真和分析

　　在本文的算法驗證實驗中，使用了Leika GPS1200測量型接收機在車載動態(tài)下采集的GPS信號的偽距和偽距增量的測量值，其中偽距測量值經(jīng)過偽距增量測量值的平滑。算法的驗證是在PC機上使用Matlab完成的。

　　在算法驗證中，設置采樣時間間隔為T=1 s。接收機的運動狀態(tài)為從靜止狀態(tài)變化到在地面上高速運動又逐步靜止。測試數(shù)據(jù)長度為2 500 s。參考軌跡是接收機使用動態(tài)實時差分（RTK）方法測定的，方差至毫米級。分別在有效衛(wèi)星不少于4顆和只有3顆的情況下，比較了幾種定位解算算法的定位結(jié)果。

　　3.1 有效衛(wèi)星數(shù)不少于4顆時的定位結(jié)果

　　以有效衛(wèi)星數(shù)為4顆為例，分別使用ILS、EKF和改進的EKF算法進行定位，得到以ECEF坐標系下XYZ三維坐標表示的定位結(jié)果，如圖1所示。

　　由圖 1可知，接收機在前400 s時基本保持靜止，然后開始運動至1 600 s，逐步靜止至測試結(jié)束。三種解算方法都可以較好地進行定位解算。為了更好地分析定位結(jié)果，畫出圖 1中的XYZ三維坐標上的定位誤差的均方根，如圖 2所示。

　　為了更好地比較ILS、普通EKF與改進EKF算法，計算圖 2中三種算法在XYZ三個坐標方向上的定位結(jié)果的均方根誤差的平均值，如表1所示。

　　從圖2和表1可以更清晰地看出，在有效衛(wèi)星不少于4顆時，三種算法在以ECEF坐標系表示的XYZ方向上的定位誤差的均方根均不超過20 m，平均定位誤差不超過10 m，定位解算都很好。其中改進的EKF與普通EKF定位相當，而ILS的定位均方根誤差抖動很小。這是因為，每當有新的測量數(shù)據(jù)時，ILS算法都會經(jīng)過多次迭代計算直至結(jié)果收斂，而EKF與改進的EKF算法均只是利用新息進行計算，計算量要遠遠小于ILS，故而定位結(jié)果會稍有抖動。

　　3.2 有效衛(wèi)星數(shù)只有3顆時的定位結(jié)果

　　在有效衛(wèi)星數(shù)只有3顆時，ILS無法進行定位。分別使用普通EKF和改進的EKF算法進行定位，得到在ECEF坐標系下的XYZ三維坐標上的定位結(jié)果，如圖 3所示。

　　為了更好地分析定位結(jié)果，畫出圖 3中的XYZ三維坐標上的定位誤差的均方根，如圖 4所示。

　　為了更好地比較普通EKF算法與改進EKF算法，計算圖 3中兩種算法XYZ三維坐標上的定位結(jié)果的均方根誤差的平均值，如表 2所示。

　　由圖 3、圖 4和表 2可知，在只有3顆有效衛(wèi)星的情況下，普通EKF的定位結(jié)果會出現(xiàn)明顯的偏移，而本文提出的改進的EKF算法的定位結(jié)果與參考軌跡擬合得很好，各個方向上的平均定位誤差均不超過10 m，定位很好。這說明本文提出的改進EKF算法在只有3顆有效衛(wèi)星時，極大地提高了EKF的定位。與理論結(jié)果吻合。

　　3.3 有效衛(wèi)星為3顆時改進EKF與有效衛(wèi)星為4顆時普通EKF 的定位結(jié)果

　　為了進一步分析改進EKF算法的性能，將改進EKF算法在只有3顆有效衛(wèi)星時的定位誤差與普通EKF算法在有4顆有效衛(wèi)星時的定位誤差進行比較，如圖 5所示。

　　計算圖 5中的兩種情況下XYZ三維坐標上的平均誤差的均方根，如表3所示。

　　由圖 5和表 3可知，在只有3顆有效衛(wèi)星時，本文提出的改進EKF算法所增加的先驗信息與增加一顆衛(wèi)星的測量信息的效果基本一致。故而說明在出現(xiàn)GPS信號阻塞的情況下，只有3顆有效衛(wèi)星時，此改進EKF算法的定位解算非常好。

　　綜上所述，利用實測的衛(wèi)星數(shù)據(jù)進行驗證，結(jié)果表明：在3顆有效衛(wèi)星的情況下，在靜止狀態(tài)和車載狀態(tài)下，改進的EKF算法的定位均明顯優(yōu)于普通的EKF，且與4顆有效衛(wèi)星時EKF的定位相當；在不少于4顆有效衛(wèi)星的情況下，在靜止狀態(tài)和車載狀態(tài)下，改進EKF算法的定位與普通EKF相當。由此說明，本文提出增加的先驗信息在已有足夠測量信息時，并不會影響定位，而只有在3顆有效衛(wèi)星時，測量信息對于定位的改進與一顆有效衛(wèi)星的測量信息效果基本一致，從而大大提高了EKF算法的魯棒性。

　　總結(jié)

　　本文介紹當GPS信號發(fā)生阻塞時，可見衛(wèi)星數(shù)會不足4顆，這時迭代二乘算法（ILS）與普通擴展卡爾曼濾波（EKF）都不再適用。針對這一問題提出了一種改進的EKF算法。該算法利用在垂直地面方向上的位置變化緩慢這一運動特性，建立了改進EKF算法的系統(tǒng)模型。通過理論分析得到了濾波器參數(shù)，利用真實的GPS衛(wèi)星數(shù)據(jù)進行驗證。實驗結(jié)果表明，在可見衛(wèi)星數(shù)不少于4顆時，此改進的EKF算法定位與普通的EKF算法基本相同；在GPS信號阻塞只有3顆可見衛(wèi)星時，此改進的EKF算法的定位明顯優(yōu)于普通的EKF算法。本文提出的改進EKF算法適合車載等地面運用，不適合在垂直地面方向上高速運動的情形。