王天雄,陳仕進(jìn),張玉鵬,閆勝虎
(中國電子科技集團(tuán)公司第五十四研究所,河北 石家莊 050081)
時間同步是時分多址(TDMA)網(wǎng)絡(luò)正常運行的基礎(chǔ),網(wǎng)內(nèi)各個節(jié)點的時間同步精度達(dá)到微秒量級,即可保證節(jié)點間的無碰撞同步信息傳輸,常用的時間同步算法均能滿足需求,例如:參考廣播方法[1]、RTT時間同步算法[2-3]和PTP時間同步算法[4]。文獻(xiàn)[5-6]針對自組織網(wǎng)絡(luò)研究了分布式同步算法,雖然時間同步精度不高,但是無需中心節(jié)點,網(wǎng)絡(luò)抗毀性能高。若要利用TDMA數(shù)據(jù)鏈路進(jìn)行節(jié)點間的測距,則需要提高網(wǎng)內(nèi)時間同步精度,例如:30 m的測距誤差就要求時間同步精度達(dá)到100 ns以內(nèi)。為了提高時間同步精度,文獻(xiàn)[7]設(shè)計了基于卡爾曼濾波的高精度RTT時間同步算法,將時間同步誤差控制在50 ns以下。文獻(xiàn)[8]將卡爾曼濾波和RTT-PTP算法相結(jié)合,在Link22數(shù)據(jù)鏈結(jié)構(gòu)的基礎(chǔ)上降低了時間偏差率誤差。
隨著多戰(zhàn)斗單元協(xié)同作戰(zhàn)理論和體系的發(fā)展,數(shù)據(jù)鏈路的時間同步面臨著更嚴(yán)峻的挑戰(zhàn)。例如:戰(zhàn)斗單元執(zhí)行協(xié)同探測、協(xié)同偵查等任務(wù)時需要各個戰(zhàn)斗單元間納秒級的時間同步[9-11],而衛(wèi)星導(dǎo)航信號一旦被干擾,外部時鐘源無法進(jìn)行授時,這就需要靠數(shù)據(jù)鏈路完成納秒級的精密時間同步,為其他作戰(zhàn)載荷提供時統(tǒng),支撐協(xié)同作戰(zhàn)任務(wù)。
本文通過研究TDMA網(wǎng)絡(luò)的RTT同步技術(shù),分析同步誤差產(chǎn)生的原因。設(shè)計了基于線性調(diào)頻(LFM)脈沖的時間同步算法,利用LFM脈沖壓縮的抗多普勒優(yōu)勢[12-13],通過采用大時間帶寬積的LFM脈沖提高了時隙到達(dá)時間的估計精度。仿真和實際平臺驗證表明該算法行之有效。
RTT時間同步方法是通過測量網(wǎng)絡(luò)內(nèi)節(jié)點之間的信號往返時間來估計時間同步誤差的。在估計時間誤差的同時,還能夠估計出節(jié)點間的信號傳輸時延,因此,常用于TDMA網(wǎng)絡(luò)的節(jié)點間測距。
RTT時間同步原理如圖1所示。其中,tm為主節(jié)點(時鐘參考節(jié)點)時間,tc為從節(jié)點時間,Δt為主從節(jié)點的時鐘誤差。通過RTT算法可以將時間tc與時間tm同步。
圖1 RTT時間同步原理
從節(jié)點在其本地時間tc1時發(fā)送RTT詢問時隙,主節(jié)點在其本地時間tm1時收到來自從節(jié)點的RTT詢問時隙,時隙傳輸時延為tp1。主節(jié)點在其本地時間tm2發(fā)送RTT應(yīng)答時隙,從節(jié)點在其本地時間tc2時收到來自主節(jié)點的RTT應(yīng)答時隙,時隙傳輸時延為tp2。
假設(shè)主從節(jié)點交互時隙的時間非常短,主從節(jié)點間的距離可以近似認(rèn)為不變,即可以取tp1=tp2=tp,此時,以上各參數(shù)滿足以下方程:
(1)
從而有:
Δt=(tm1+tm2-tc1-tc2)/2,
(2)
(3)
分析上述RTT時間誤差估計流程可以看出,在不考慮硬件一致性等原因引起的系統(tǒng)誤差的前提下,僅從算法本身考慮,主要存在以下估計誤差:在估計tm1,tc2的時候存在估計誤差εm,εs;由主從節(jié)點的時鐘相對漂移引入的隨機誤差εc;主從節(jié)點交互脈沖時由于主從節(jié)點的相對位置變化會引起tp1≠tp2,因此會造成產(chǎn)生估計誤差εtp。
目前減小上述誤差的途徑主要有3條:在硬件上提高用作系統(tǒng)時鐘的晶振穩(wěn)定度,減小εc的方差;在軟件上使用卡爾曼濾波等估計算法,減小由εc,εm,εs引起的估計誤差方差;在幀結(jié)構(gòu)設(shè)計上,盡量減小主從節(jié)點的時隙交互時間,減小由εtp引起的估計誤差。
然而,由于εm,εs的方差較大,即使采用了相應(yīng)的估計算法,時間同步誤差也只能控制在幾十納秒量級。要想進(jìn)一步提高時間同步精度,需要從根本上減小εm,εs的方差。當(dāng)信號帶寬受限時,增加信號持續(xù)時間、提高信號的時間帶寬積是提高信號到達(dá)時間估計精度的有效手段。對于常規(guī)的偽隨機序列調(diào)制信號,提高信號持續(xù)時間就意味著增加序列碼長。但在多普勒頻偏條件下,對偽隨機序列進(jìn)行長碼匹配濾波,其累積能量損失會隨著其碼長的增加而增大,累積能量信噪比的降低直接影響信號的捕獲概率和到達(dá)時間的估計準(zhǔn)確度[14-16]。若采用差分能量累積算法,則會額外損失信噪比,要想達(dá)到與相干能量累積相同的效果,需要進(jìn)一步加長偽隨機序列,增加額外時間開銷。若采用分段匹配加FFT的能量累積算法,額外信噪比損失較小,但會大大增加系統(tǒng)的復(fù)雜度,且原有的時域一維峰值搜索將變?yōu)闀r頻二維峰值搜索,增加了算法的復(fù)雜度。因此,采用偽隨機長碼序列調(diào)制信號進(jìn)行信號到達(dá)時間估計存在很大的局限性。
線性調(diào)頻信號(LFM)對多普勒頻偏有天然的抗性,文獻(xiàn)[17-18]指出,通過脈沖壓縮方法對LFM信號進(jìn)行能量累積時,由多普勒頻偏引起的能量損失只與頻偏帶寬比有關(guān),與信號持續(xù)時間無關(guān)。因此,可以對大時間帶寬積的LFM信號直接進(jìn)行脈沖壓縮,并利用能量峰值位置準(zhǔn)確估計脈沖到達(dá)時間,算法復(fù)雜度和實現(xiàn)復(fù)雜度均較低。
線性調(diào)頻信號可以表示為:
(4)
式中,A為幅度,f0為初始頻率,k為調(diào)頻率。設(shè)線性調(diào)頻信號的持續(xù)時間為τ,其理想的時頻分布如圖2所示。
圖2 線性調(diào)頻信號的理想時頻分布
從圖2中能夠看出,線性調(diào)頻信號的時寬為T=τ,其頻率隨時間線性變化,因此是非平穩(wěn)信號,其瞬時頻率為:
f=f0+kt。
(5)
脈沖壓縮是雷達(dá)術(shù)語,在通信領(lǐng)域,即為匹配濾波。因此,LFM脈沖的捕獲和到達(dá)時間估計,是通過搜索信號匹配濾波響應(yīng)峰值來實現(xiàn)的[12]。
線性調(diào)頻信號對多普勒頻移不敏感,在較大多普勒頻移的情況下,仍能通過匹配濾波得到很好的能量累積,但匹配濾波器響應(yīng)會產(chǎn)生附加時延,延時的大小正比于多普勒頻移,延時會使脈沖到達(dá)時間的估計產(chǎn)生誤差。
LFM多普勒頻移影響示意圖如圖3所示。LFM脈沖附加了多普勒頻移之后,相對于原有頻率附加了頻率fd。實線部分為原始頻率變化范圍,虛線部分為頻移后的頻率變化范圍,因為只有頻率在f1~f2部分的波形可以通過脈沖壓縮收集能量,所以引入多普勒頻移后,將會產(chǎn)生以下不利影響:
① 可壓縮段變小了,將會損失一定的累積能量,損失部分可以表示為Δtd/T;
② 匹配濾波器輸出響應(yīng)產(chǎn)生附加時延Δtd=Tfd/(f1-f2)。
圖3 LFM多普勒頻移影響示意圖
通過上述分析可知,當(dāng)實際的多普勒頻移與LFM信號帶寬比較小時,產(chǎn)生的能量損失可以忽略不計。但產(chǎn)生的附加時延必需予以補償。正負(fù)斜率脈沖多普勒補償算法可以有效地修正附加時延。
正負(fù)斜率脈沖補償算法要求發(fā)送端連續(xù)發(fā)送正斜率脈沖和負(fù)斜率脈沖2個LFM脈沖,即2個脈沖信號的調(diào)頻斜率幅值相等,方向相反。在接收端需要2個匹配濾波器分別匹配正斜率脈沖和負(fù)斜率脈沖。當(dāng)發(fā)射正斜率脈沖時,多普勒頻移使脈壓后的波形包絡(luò)在時間軸上發(fā)生移動,偏移量為Δtd1=-2πfd/k。同理,發(fā)射負(fù)斜率脈沖時,多普勒頻移使脈壓后的波形包絡(luò)在時間軸上發(fā)生移動,偏移量為Δtd2=2πfd/k,其原理如圖4所示。
圖4 多普勒對LFM正負(fù)斜率脈沖的影響
為了利用LFM脈沖進(jìn)行信號到達(dá)時間估計的優(yōu)勢。TDMA網(wǎng)絡(luò)時幀結(jié)構(gòu)如圖5所示。一個時幀被分為2部分,前半部分為LFM脈沖部分,用于節(jié)點間進(jìn)行信號往返時間估計,每個節(jié)點發(fā)送一個LFM脈沖串#n′,每個脈沖串包含若干個正負(fù)斜率脈沖對;后半部分為通信波形部分,用于節(jié)點間進(jìn)行正常通信和傳輸LFM脈沖串到達(dá)時間的估計結(jié)果。
圖5 系統(tǒng)時幀結(jié)構(gòu)圖
將LFM脈沖集中設(shè)計在時幀的前區(qū),是為了盡可能地減小節(jié)點間信號往返交互的時間,從而減小誤差εtp的影響。
基于上述時幀結(jié)構(gòu),主從節(jié)點工作流程如圖6所示。
圖6 主從節(jié)點工作流程
主節(jié)點流程描述為:
① 系統(tǒng)初始化后,在tm=0時產(chǎn)生幀起始脈沖,為通信時隙和LFM脈沖時隙提供時間基準(zhǔn);
③ 在tm為tm2時發(fā)送RTT應(yīng)答脈沖串;
由于LFM脈沖串本身沒有調(diào)制信息,因此,從LFM信號本身無法區(qū)分脈沖的信源,所以,在主從節(jié)點時間不同步時,從節(jié)點無法辨認(rèn)主節(jié)點發(fā)送的LFM脈沖串,也就無法進(jìn)行時間的初始化校準(zhǔn)。為了克服上述問題,從節(jié)點借助通信時隙來完成時間的初始化校準(zhǔn),因此,從節(jié)點的時間同步分為時間粗同步和時間細(xì)同步2個步驟。
在時間粗同步階段,從節(jié)點工作流程描述為:
① 系統(tǒng)初始化后從節(jié)點保持靜默,接收主節(jié)點通信時隙,從中恢復(fù)時間的粗定時;
② 在粗同步的時間下,根據(jù)脈沖串到達(dá)時間區(qū)分脈沖串的信源,并進(jìn)行本節(jié)點各LFM脈沖時隙和通信時隙的發(fā)送。
在時間細(xì)同步階段,從節(jié)點工作流程描述為:
① 在tc=0時產(chǎn)生幀起始脈沖,為通信時隙和LFM脈沖調(diào)制提供時間基準(zhǔn);
② 在tc為tc1時發(fā)送RTT詢問脈沖串;
⑤ 按式(3)計算Δt,并修正本地時間;
本文針對上述TDMA網(wǎng)絡(luò)時間同步方法,以4節(jié)點網(wǎng)絡(luò)進(jìn)行了仿真和實現(xiàn)驗證。
仿真條件1:校時幀周期為40 ms,其中LFM脈沖部分占3 ms,通信波形部分占37 ms,每個LFM脈沖串包含8個正負(fù)斜率脈沖對,每個脈沖脈寬40 μs;信號帶寬20 MHz;多普勒頻偏20 kHz;系統(tǒng)時鐘晶振穩(wěn)定度10-8;節(jié)點相對運動速度不大于30 m/s。主從節(jié)點的時間同步誤差如圖7所示。
圖7 條件1下主從節(jié)點的同步誤差圖
從圖7可以看出,在接收信噪比達(dá)到-14 dB時,時間同步誤差開始收斂,當(dāng)信噪比達(dá)到-8 dB時,同步誤差標(biāo)準(zhǔn)差不大于5.5 ns。
仿真條件2:在條件1的基礎(chǔ)上分別將LFM脈沖帶寬和脈沖寬度降低到10 MHz,20 μs,仿真結(jié)果如圖8所示。從圖8可以看出,由于降低LFM的帶寬或脈寬減小了脈沖的時間帶寬積,致使時間同步靈敏度下降,-8 dB時,時間同步誤差標(biāo)準(zhǔn)差增大到10 ns左右。
仿真條件3:在仿真條件1的基礎(chǔ)上將多普勒頻偏增加到200 kHz,仿真結(jié)果如圖9所示。
圖8 條件2下主從節(jié)點時間同步誤差圖
圖9 條件3下主從節(jié)點時間同步誤差圖
比較圖7和圖9可知,LFM信號具有良好的抗多普勒性能,在200 kHz多普勒頻偏下,時間同步靈敏度和時間同步精度均未有明顯的損失。
仿真條件4:在條件1的基礎(chǔ)上分別將時鐘穩(wěn)定度調(diào)制為10-7,節(jié)點間相對運動速度調(diào)整為300 m/s,校時幀周期調(diào)整為1 s,仿真結(jié)果如圖10所示。
圖10 條件4下主從節(jié)點時間同步誤差圖
從圖10可以看出,時鐘穩(wěn)定度、相對運動速度和校時幀周期直接影響時間同步精度,為了減小時間同步精度損失,應(yīng)選用高穩(wěn)定度的時鐘和減小校時周期。
根據(jù)上述方法設(shè)計FPGA程序?qū)崿F(xiàn),并加載到現(xiàn)有的通信終端設(shè)備上,實驗條件是仿真條件1的基礎(chǔ)上去掉節(jié)點間的相對運動速度,利用示波器來觀察實際的主從節(jié)點幀起始脈沖的同步情況。
實驗結(jié)果顯示,從節(jié)點的幀起始脈沖在主節(jié)點幀起始脈沖附近波動,誤差介于±4 ns左右之間,因為沒有模擬節(jié)點間的相對運動,所以實測結(jié)果要好于仿真結(jié)果。
本文利用LFM信號抗多普勒性能強的特點,基于LFM信號設(shè)計了TDMA網(wǎng)絡(luò)的RTT時間同步算法。LFM信號的高時間帶寬積提高了接收端脈沖壓縮的能量累積量,減小了脈沖到達(dá)時間的估計誤差;節(jié)點間集中交互LFM脈沖的幀結(jié)構(gòu),減小了脈沖交互時間,降低了節(jié)點間相對運動引入的時間同步誤差。仿真和實際平臺驗證結(jié)果表明,該方法在大多普勒頻移環(huán)境下仍能夠大幅提高TDMA網(wǎng)絡(luò)內(nèi)節(jié)點的時間同步精度,可以為其他任務(wù)載荷提供高精度授時,特別適用于小規(guī)模協(xié)同作戰(zhàn)。