基于GNSS的通信星座系統(tǒng)的時間同步研究

黃 琳 馬文杰 李國軍 霍德聰 蒙 薇 楊 芳

(航天東方紅衛(wèi)星有限公司,北京 100094)

1 引言

2 基于GNSS的通信衛(wèi)星星座系統(tǒng)時統(tǒng)技術(shù)方案

全網(wǎng)時間同步技術(shù)是衛(wèi)星星座系統(tǒng)的核心技術(shù)之一。就同步內(nèi)容而言,對不同時鐘的秒脈沖相位差進(jìn)行調(diào)整是一個物理過程,稱為物理同步;若不進(jìn)行相位差調(diào)整,而只是廣播相位差以使有關(guān)部分進(jìn)行同步則是一個數(shù)學(xué)過程,稱為數(shù)學(xué)同步。物理同步和數(shù)學(xué)同步構(gòu)成了時間同步的兩大基本內(nèi)容[1]。與導(dǎo)航星座系統(tǒng)重點關(guān)注星間數(shù)學(xué)同步[2-3]的情況不同,本文更關(guān)注通信衛(wèi)星星座系統(tǒng)的空間段與地面段的時間同步問題,即,各星上時鐘與地面系統(tǒng)時鐘的秒脈沖相位同步。

當(dāng)前,國外著名的通信衛(wèi)星星座系統(tǒng)的星地時間同步處理,基本上都依賴全球布設(shè)的地面站系統(tǒng)來完成。比如,具有星間鏈路的銥星(Iridium)系統(tǒng)有11個全球分布的地面站;不具有星間鏈路的全球星(Globalstar)系統(tǒng)需要300 多個地面站。然而,由于我國地面站的分布區(qū)域有限,地面站只能在衛(wèi)星過境的短時期內(nèi)觀測到星地時差,為了實現(xiàn)高精度同步,必然對硬件性能(比如,時鐘)和星地時差預(yù)報模型的估計精度提出過高要求,從而導(dǎo)致星座系統(tǒng)總成本和工程實現(xiàn)的困難程度大大上升。顯然,這種以地面站介入控制為主的時統(tǒng)方法不太適合我國的星座系統(tǒng)的星地時間同步處理。

為了提高星座系統(tǒng)時統(tǒng)操作的可靠性,且考慮到當(dāng)前絕大多數(shù)衛(wèi)星都裝備某種全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)(比如GPS 或GLONASS)接收機(jī),因此,基于GNSS系統(tǒng)實現(xiàn)星座系統(tǒng)的時間同步是一個可能選擇。比如,G PS系統(tǒng)作為一個高精度、高穩(wěn)定度時間基準(zhǔn),已成功用于相距數(shù)千千米之遙的地面站間時鐘同步[4]。本文將引入星上時間管理單元自主完成星地時統(tǒng)處理的思想,提出一個基于GNSS系統(tǒng)實現(xiàn)近地全球通信衛(wèi)星星座系統(tǒng)全網(wǎng)星地同步的技術(shù)方案,重點設(shè)計一個半實物仿真方案,并利用數(shù)學(xué)、半物理仿真對該技術(shù)方案的正確性與可行性進(jìn)行了驗證。

2.1 方案概述與時統(tǒng)關(guān)鍵技術(shù)

時統(tǒng)方案主要內(nèi)容包括:地面站和各星座星上的GNSS 接收機(jī)每隔一定周期測得地面系統(tǒng)時鐘、各星上時鐘相對于GNSS系統(tǒng)時鐘的相位差。而地面系統(tǒng)時鐘相對GNSS系統(tǒng)時鐘的時差觀測量經(jīng)過星地、星間鏈路,被發(fā)送給任一星上時間管理單元,并與各星座星上測得的時差量一起進(jìn)行折算,得到各星時鐘相對地面系統(tǒng)時鐘的時差觀測量。各星上時間管理單元中的卡爾曼濾波器(KF),周期性估計出星地時差預(yù)報模型中的幾個變量(時差量、頻差和頻漂率)、并周期性利用新估計結(jié)果替換掉時差預(yù)報模型的原有變量估計。同時,時間管理單元基于星上時差預(yù)報模型進(jìn)行時差廣播,并進(jìn)行自主調(diào)相時機(jī)的判斷,若調(diào)相時機(jī)到來,則自動發(fā)出指令驅(qū)動秒脈沖發(fā)生器進(jìn)行相應(yīng)相位差的調(diào)整,從而將各星上時鐘相對地面系統(tǒng)時鐘的秒脈沖相位差控制在同步指標(biāo)要求的范圍內(nèi)。可見,該星座時統(tǒng)技術(shù)方案具有很大程度上的自主性。

星上時間管理單元自主時統(tǒng)處理的關(guān)鍵,在于時間管理單元具備準(zhǔn)確跟蹤真實的星地時差Δαk變化趨勢的能力,星上采用如下二次模型來跟蹤真實星地時差的時變趨勢[5]。

其中,Δα0表示起始預(yù)報時刻t0的星地時差量;α1和α2分別表示星上時鐘相對地面系統(tǒng)時鐘的頻差和頻漂率估計值。然而,依據(jù)上述模型對星地時差進(jìn)行跟蹤,其跟蹤誤差將隨預(yù)報時間間隔tk-t0的增加而快速增加(tk為k個預(yù)報時刻)。為了保證長期準(zhǔn)確跟蹤時差,需要利用不時獲得的星地時差觀測量,通過在線估計器,對模型變量進(jìn)行在線估計,并周期性地替換原有的模型變量估計(需要改變起始預(yù)報時刻t0)。

2.2 星地時差觀測

地面站和星座各星上的GNSS 接收機(jī)(其內(nèi)部時鐘被屏蔽,實際采用地面站或星上時鐘信號)各自測得地面系統(tǒng)時鐘、各星上時鐘相對于GNSS系統(tǒng)時鐘的相位差。其原理為GNSS 定位原理[5]

在得到Δαjsys觀測量之后,記地面系統(tǒng)時鐘相對GNSS系統(tǒng)時鐘的時差為Δα0sys,星座星的時鐘相對GNSS系統(tǒng)時鐘的時差為Δαjsys(j =1,2,…,n),則每顆星座星時鐘相對地面系統(tǒng)時鐘的時差觀測量為

2.3 時差預(yù)報模型變量估計

時差模型變量估計是一個狀態(tài)估計問題,其狀態(tài)量為

離散的狀態(tài)方程和觀測方程為

其中

Δτ表示離散周期;wk表示時鐘的隨機(jī)過程噪聲;vk表示隨機(jī)時差觀測噪聲;xk表示k時刻的狀態(tài)量;zk表示k時刻的觀測量。

上述狀態(tài)估計問題通常可用文獻(xiàn)[6]給出的KF 來解決。

2.4 自主調(diào)相處理

當(dāng)星地秒脈沖相位差滿足如下條件時

星上時間管理單元判斷調(diào)相時機(jī)已到。其中,Δαlim和Δαgate分別稱為調(diào)相指標(biāo)和調(diào)相閾值(其值大小反映調(diào)相的敏感度)。由于當(dāng)前真實時差量未知,所以需要利用時差預(yù)報量(或估計量)。

采用名義調(diào)相殘差最小策略,需要調(diào)整的相位差為

其中,f0為時鐘基頻;round(·)表示按照四舍五入原則取整操作;相應(yīng)的名義調(diào)相殘差為

3 半實物仿真方案設(shè)計

3.1 任務(wù)場景設(shè)想

一個近地通信衛(wèi)星星座,是典型的Walker(24/3/1)星座,衛(wèi)星軌道高度1 000km,并考慮一個地面站,以10°最低觀測仰角計算,該系統(tǒng)可以保障任何時候都至少有1顆衛(wèi)星在地面站的服務(wù)范圍內(nèi),同軌面內(nèi)的相鄰星間時刻存在穩(wěn)定持續(xù)的鏈路,而任意兩個異軌道面之間也有適宜的星間鏈路,因此,該星座系統(tǒng)可以保證星座任意衛(wèi)星與地面站之間處于時時刻刻的連通狀態(tài)。為了實現(xiàn)星座系統(tǒng)任務(wù),要求星座各衛(wèi)星與地面站的秒脈沖相位差保持在一定范圍內(nèi)。假設(shè)兩者同步精度要求為100ns。

3.2 半實物仿真方案設(shè)想

考慮到該時統(tǒng)方案下星間鏈路不參與星間時差觀測,而只是將地面GPS時差觀測數(shù)據(jù)傳輸?shù)浇匦亲餍l(wèi)星,可以說,各近地衛(wèi)星的星地時統(tǒng)是相對獨立的,因此,對一顆近地衛(wèi)星的星地時統(tǒng)進(jìn)行考察,可以在很大程度上反映出整個近地星座系統(tǒng)的星地時統(tǒng)能力。考慮到所需硬件成本和數(shù)量的要求很高,本文僅提出一個半實物仿真原理方案,如圖1所示。

上述仿真平臺架構(gòu)方案共包括三部分內(nèi)容。

1)星上(單星)相關(guān)時統(tǒng)部分:1臺GPS 接收機(jī)和1臺時間管理單元(外接1臺時鐘)。為了彌補(bǔ)現(xiàn)有時間管理單元數(shù)據(jù)處理能力的不足,利用1臺工控機(jī)完成時間管理單元中的時差數(shù)據(jù)在線處理工作。

2)地上(單地面站)相關(guān)時統(tǒng)部分:1臺GPS 接收機(jī)和1臺時間管理單元(外接1臺時鐘)。同樣,利用1臺工控機(jī)完成時間管理單元中的時差數(shù)據(jù)在線處理工作。

圖1 時間同步半實物仿真原理圖Fig.1 Plan of GPS-based time synchronization hardware-in-the-loop simulation

3)輔助部分:1臺GPS系統(tǒng)仿真器和1 套仿真平臺供電設(shè)備以及1臺真實時差觀測設(shè)備。此外,地面工控機(jī)還承擔(dān)仿真管控的作用,也是輔助仿真的設(shè)備之一。

考慮到現(xiàn)有條件,星上和地上所用的GPS 接收機(jī)、時間管理單元和工控機(jī)都采用了同類型產(chǎn)品,不同之處在于星上時間管理單元采用外接銣鐘,而地上則外接了更高性能的銫鐘。有關(guān)設(shè)備的配置情況見表1。

根據(jù)原理圖1,將表1所示設(shè)備連接起來的半實物仿真平臺如圖2所示。

表1 仿真設(shè)備Table1 Chosen equipments for the hardwarein-the-loop simulation

圖2 時間同步半實物仿真設(shè)備連接圖Fig.2 Hardware-in-the-loop simulation of GPS-based time synchronization

4 半物理仿真結(jié)果與分析

4.1 仿真設(shè)置

半物理仿真按照圖3所示進(jìn)行參數(shù)設(shè)置。其中,星地秒脈沖同步(調(diào)相)指標(biāo)和調(diào)相閾值分別為

周期設(shè)置:系統(tǒng)仿真周期、時差觀測周期、時差預(yù)報模型變量替換周期和時差預(yù)報周期皆設(shè)置為4s。

圖3 仿真參數(shù)設(shè)置Fig.3 Parameters set ting of simulation

至于數(shù)學(xué)仿真,還需要仿真出地面站時鐘相對GPS系統(tǒng)時鐘的時差觀測量,以及星上時鐘相對GPS系統(tǒng)時鐘的時差觀測量。由于真實參數(shù)和模型未知,根據(jù)星地時鐘產(chǎn)品的名義性能參數(shù)和GPS系統(tǒng)時鐘的大致性能了解,并依據(jù)式(1)所示的簡單模型仿真出這些時差觀測量信息。隨機(jī)觀測噪聲選擇白噪聲,均方差假設(shè)為50ns。

4.2 仿真結(jié)果與分析

數(shù)學(xué)仿真結(jié)果如圖4～圖6所示。其中,圖4(上)給出了基于GPS系統(tǒng)的星地時差觀測量,受到如圖4(下)所示的100ns 量級的隨機(jī)觀測噪聲的干擾影響。圖5(上)給出了在線時差濾波器的星地秒脈沖相位偏差估計的誤差,結(jié)果顯然是一致性收斂的,圖5(下)給出了星上時差預(yù)報模型的誤差。不難看出,這兩圖趨勢幾乎一致,原因在于,時差預(yù)報模型變量替換周期等于觀測更新周期,即,每個更新時刻的濾波估計值都用于替換預(yù)報模型的上一個參數(shù)估計。圖6給出了不采用(上)/采用(下)自主調(diào)相措施下的星地秒脈沖相位差。不難看出,若無自主調(diào)相措施,星地秒脈沖不可能控制在指標(biāo)要求的范圍內(nèi)。

半實物仿真給出的星地秒脈沖相位差估計和觀測量結(jié)果如圖7所示,可以看出:

1)同數(shù)學(xué)仿真一樣,星地秒脈沖相位差的絕對值一直被控制在100ns 以內(nèi),這說明星上自主時統(tǒng)處理成功實現(xiàn)了對調(diào)相時機(jī)的判斷、調(diào)相量的計算以及調(diào)相指令的發(fā)送,而時間管理單元的秒脈沖發(fā)生器也成功地進(jìn)行了相應(yīng)的調(diào)整。

2)星地秒脈沖相位差的觀測量隨機(jī)分布在線性增長的真實相位差附近有限區(qū)域內(nèi),分布區(qū)域大小顯示隨機(jī)觀測誤差高于100ns(尚不包括無數(shù)據(jù)輸出的異常情況)。不難理解,單純的星地時差觀測無法準(zhǔn)確跟蹤星地秒脈沖相位差,也無法據(jù)此達(dá)到相位同步指標(biāo)要求。上述結(jié)果反過來說明,濾波器很好地估計出真實相位差及其變化趨勢,故而時間管理單元才能準(zhǔn)確判斷出自主調(diào)相時機(jī),計算調(diào)相量并發(fā)出調(diào)相指令。

3)星上自主時統(tǒng)處理流程成功實現(xiàn)了單星的時統(tǒng)目標(biāo)。由于這種時統(tǒng)方案下各星同步相對獨立,只要星間鏈路能夠發(fā)揮數(shù)據(jù)傳輸?shù)墓δ?本文所提的時統(tǒng)方法可望實現(xiàn)整個星座系統(tǒng)的星地時統(tǒng)目 標(biāo)。

圖4 星地時差觀測量(上)與隨機(jī)觀測誤差(下)Fig.4 Observations(upper)and random errors(lower)

圖5 時差估計(上)與時差預(yù)報誤差(下)Fig.5 Estimation errors(upper)and prediction errors(lower)

需要指出的是:

1)半實物仿真中采用的G PS 仿真器的時鐘性能相當(dāng)于高穩(wěn)晶振產(chǎn)品,遠(yuǎn)低于地面銫鐘和星上銣鐘的性能,這導(dǎo)致星地相位差的隨機(jī)觀測誤差大大增大(將近200ns),而實際的GPS導(dǎo)航系統(tǒng)時鐘相對仿真中所用的銫鐘和銣鐘來說更準(zhǔn)確、更穩(wěn)定,故而相位差的觀測精度將有所提高。換句話說,半實物仿真的硬件配置條件相比真實情況來說更惡劣,但是,這個自主時統(tǒng)方案仍然達(dá)到了星地同步指標(biāo)要求。

2)半物理仿真結(jié)果與數(shù)學(xué)仿真結(jié)果存在一些微小差別,主因在于數(shù)學(xué)仿真對星/地時鐘、GPS系統(tǒng)時間等參數(shù)的設(shè)置與半實物仿真中的情況有差別,由于真實參數(shù)/模型不可精確得知,數(shù)學(xué)仿真采用名義上的參數(shù)和簡單模型必然存在較大差別。但這兩類仿真結(jié)果都很好地說明這套時統(tǒng)方案的可行性與有效性。

圖6 未校準(zhǔn)的星地秒脈沖偏差(上)與時統(tǒng)后的偏差(下)Fig.6 Phase errors of PPS before(upper)/after(low er)calibration

圖7 半實物仿真結(jié)果采樣Fig.7 Phase errors of one satellite PPS relative to the ground PPS

5 結(jié)論

本文為近地衛(wèi)星星座系統(tǒng)提出一個利用全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)實現(xiàn)全網(wǎng)星地秒脈沖相位同步的時統(tǒng)方法,并提出一個半實物仿真方案對此加以驗證。數(shù)學(xué)仿真與半實物仿真的結(jié)果都表明,近地衛(wèi)星星座系統(tǒng)基于GPS 類GNSS系統(tǒng)可以實現(xiàn)高精度的星地時間同步,證實了星上時間管理單元自主時統(tǒng)處理流程的正確性。

隨著我國自主研發(fā)的GNSS系統(tǒng)的建立,本文所提出的星座系統(tǒng)時統(tǒng)方法可以避免系統(tǒng)安全性方面的隱患,也有望大大降低此類星座系統(tǒng)的維護(hù)成本。

References)

[1]郭海榮.導(dǎo)航衛(wèi)星原子鐘時頻特性分析理論與方法研究[D].中國人民解放軍信息工程大學(xué),2006

[2]帥平.導(dǎo)航星座的自主導(dǎo)航技術(shù)——衛(wèi)星自主時間同步[J].飛行器測控學(xué)報,2004:23(4)11-15

[3]惠衛(wèi)華,盧曉春,華宇,等.衛(wèi)星星座時間同步方法的研究和仿真[C]//西安:2005年全國時間頻率學(xué)術(shù)交流會議,2005:47-51

[4]聶桂根.高精度GPS 測時與時間傳遞的誤差分析及應(yīng)用研究[D].武漢大學(xué),2002

[5]童寶潤.時間統(tǒng)一系統(tǒng)[M].北京:國防工業(yè)出版社,2003

[6]Stein S R,Evans J.The application of Kalman filters and ARIM Amodels to the study of time prediction errors of clocks for use in the defense communication system[C]//44thAnnual Symposium on Frequency Control.1990:630-635