戴邵武,馬長里,代海霞
(1.海軍航空工程學院,山東 煙臺264001;2.北京清華同方工業(yè)有限公司,北京10083)
我國正在運行的北斗一代衛(wèi)星導航系統(tǒng)只是個初級區(qū)域導航定位系統(tǒng)。因此,我國的第二代衛(wèi)星導航系統(tǒng)——北斗二代衛(wèi)星導航系統(tǒng)已經立項并進入組網高峰期,預計在2015年形成覆蓋全球的衛(wèi)星導航定位系統(tǒng)。正在建設的“北斗二代”衛(wèi)星導航系統(tǒng)將由5顆靜止軌道衛(wèi)星(GEO)和30顆非靜止軌道(中軌(MEO)及傾斜地球同步圓軌道(IGSO))衛(wèi)星組成[1]。雖然北斗二代衛(wèi)星定位系統(tǒng)可達到很高的定位精度,但在特殊時刻和特殊地區(qū)由于信號受到遮擋容易出現(xiàn)丟星現(xiàn)象,而慣性導航則不受時間和地點的限制,可自主導航。如何使北斗二代和捷聯(lián)慣導進行組合,使之更好地完成導航定位任務,是當前我國在該領域研究的熱點問題[2]。
“北斗二代”衛(wèi)星導航系統(tǒng)的基本工作原理是[3]:空間段衛(wèi)星接收地面控制系統(tǒng)上行注入的導航電文及參數(shù),并連續(xù)向地面用戶發(fā)播衛(wèi)星導航信號,用戶收到至少四顆衛(wèi)星信號后,進行偽距測量和定位解算,最后得到定位結果。同時,為了保持地面控制系統(tǒng)各站之間的時間同步,以及地面站和衛(wèi)星之間時間同步,通過站間和星地時間比對,完成地面站間和衛(wèi)星與地面站間的時間同步。分布在國土內的監(jiān)控站負責對可視范圍內的衛(wèi)星進行監(jiān)測,采集各類觀測數(shù)據(jù)后將其發(fā)送至主控站,由主控站完成衛(wèi)星軌道精密確定及其它導航參數(shù)的確定、廣域差分信息和完好性信息處理,形成上行注入的導航電文參數(shù)。
對北斗二代/SINS組合導航系統(tǒng)的設計采取偽距、偽距率的組合模式。在組合之前要進行選星,以得到視野內最佳的衛(wèi)星幾何構型。受篇幅所限,不對選星算法做具體研究,可參考文獻[4]。該組合模式中,接收機只提供星歷數(shù)據(jù)和偽距、偽距率即可,省去導航計算處理部分,有精度高、魯棒性好、抗干擾強的特點。組合導航原理圖如圖1所示。
圖1 偽距、偽距率組合導航系統(tǒng)原理圖
捷聯(lián)慣導的誤差方程表示為[5]
系統(tǒng)的狀態(tài)向量定義為
式中:下標E、N分別代表東向和北向;L為載體所處的緯度,λ為載體所處的經度;φE、φN、φU代表東 、北、天方向的姿態(tài)誤差角;εbxεbyεbz代表東 、北、天方向的等效陀螺漂移,▽ax▽ay▽az代表東、北、天方向的等效加速度計漂移。
BeidouII的狀態(tài)向量同樣取為時鐘誤差相應的距離和時鐘頻率誤差相應的距離率,分別設為δ tu2和 δ tru2,其微分方程為
所以BeidouII/SINS組合系統(tǒng)的狀態(tài)向量為
組合導航的狀態(tài)模型為
式中:
1)偽距量測方程
由SINS載體到衛(wèi)星Sj的偽距為ρIj,可表示為[6]
設SINS位置的坐標為[x,y,z]T,將(4)式在[x,y,z]T處展開成泰勒級數(shù),取一次項:
令:
同時接收機測得的偽距為
則偽距量測值可以寫成
取j=1,2,3,4得
又由(x,y,z)與空間直角坐標系的變換關系式:
可得偽距量測方程
式中:
2)偽距率量測方程
對應于捷聯(lián)慣導給出的位置處的偽距率為
接收機測得的偽距率為
則偽距率量測方程如下:
取j=1,2,3,4得
式中:δ﹒x、δ﹒y、δ﹒z為地球直角坐標系中表示的速度誤差,可通過變換坐標系將 δ﹒x、δ﹒y、δ﹒z表示為[7]
則可獲得偽距率量測方程
式中:
式中:
將偽距量測方程和偽距率量測方程合并成組合導航的量測方程為:
整個組合導航系統(tǒng)的仿真總體結構主要包括初始參數(shù)設置、飛行器軌跡仿真、慣導輸出仿真、BeidouII系統(tǒng)的仿真、選星及卡爾曼濾波算法。仿真整體結構圖如圖2所示:
對于運動的載體,由于導航系統(tǒng)是時變的,系統(tǒng)矩陣和量測矩陣的確定都與載體的位置、速度有關,因而在仿真前,要對飛行軌跡進行設計,以提供仿真量測值。本文為了全面研究飛行載體各種機動飛行對仿真性能的影響,設計一個比較完整且機動飛行較多的飛行軌跡(包括平飛、轉彎、爬升和加速等飛行狀態(tài)),仿真時間共1400 s。
圖2 組合導航系統(tǒng)仿真流程圖
平臺初始誤差角取為東向100′,北向100′,天向200′;三維速度、位置誤差分別為1 m/s,50 m,陀螺漂移為0.1(°)/h,等效加速度計零偏為10-4g,陀螺相關時間為0.1(°)/h,陀螺白噪聲漂移為0.1(°)/h。飛機初始位置為北緯32°,東經118°,高度2000 m。初始航向0°,速度350 m/s。偽距測量的誤差涉及到衛(wèi)星時鐘誤差(由衛(wèi)星星歷獲得)、衛(wèi)星星歷誤差、電離層和對流層延遲誤差、多路徑效應誤差和接收機有關的誤差等,上述誤差經過改正及補償后可大大降低偽距測量的誤差。假設偽距測量噪聲是標準差為20 m的白噪聲。對上述飛行軌跡進行仿真分析,仿真后的航跡軌跡、平臺誤差角、速度誤差和位置誤差分別如圖3~6所示:
圖3 航跡軌跡
仿真誤差穩(wěn)態(tài)值如表1所示,結果表明,基于偽距、偽距率的SINS/BeidouII組合導航系統(tǒng)可有效提高導航精度。為我國北斗二代和捷聯(lián)慣導組合導航的研制、開發(fā)提供了思路。
表1 基于偽距、偽距率組合導航系統(tǒng)的穩(wěn)態(tài)誤差值
仿真結果表明,該組合算法可得到很高的導航精度,對環(huán)境的變化也有較強的適應性,有利于工程應用,為北斗二代衛(wèi)星和捷聯(lián)慣導的組合提供了思路。
[1] 石 華.中國北斗二代覆蓋全球[N].環(huán)球時報,2007-03-25(3).
[2] 喬 洋,趙育善,趙金才.慣性導航/雙星組合導航的可行性研究[J].上海航天,2003,(4):34-37
[3] 隨 藝.北斗衛(wèi)星和偽衛(wèi)星組合定位系統(tǒng)研究[D].西北工業(yè)大學,西安:2004
[4] 袁建平,羅建軍等.衛(wèi)星導航原理與應用[M].北京:中國宇航出版社,2004.
[5] 林雪原,劉建業(yè),范勝林.三星定位/SINS偽距組合導航系統(tǒng)的研究[J].中國空間科學技術,2002,(6):6-11.
[6] 林雪原,何 友.雙星定位/捷聯(lián)慣導組合導航技術[J].上海交通大學學報,2005,39(9):1485-1488.
[7] 張友安,林雪原,徐勝紅.綜合導航與制導系統(tǒng)[M].北京:海潮出版社,2005,306-307.