李文剛, 王屹偉

(西安電子科技大學通信工程學院, 陜西 西安 710071)

0 引 言

高動態(tài)條件下,全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GNSS)的信號將會帶有較大的多頻勒頻率偏移,如果載體的運動具有加速度,其多普勒頻率偏移也將會處于變化之中,這對于GNSS信號的處理十分不利[1-2]。在GNSS信號被捕獲之后,需要對GNSS信號進行導(dǎo)航數(shù)據(jù)的比特同步,才能夠?qū)NSS信號進行跟蹤和解調(diào)等后續(xù)處理[3]。以北斗B1頻點為例,其擴頻碼的周期為1 ms,中低軌道的北斗衛(wèi)星數(shù)據(jù)傳送速率為50 bit/s,因而每一個導(dǎo)航電文持續(xù)的時間為20 ms,如果捕獲過程中相干累加時間為1 ms,那么存在20個可能的導(dǎo)航數(shù)據(jù)的比特跳變位置[4]?，F(xiàn)代GNSS接收機比特同步過程如圖1所示,將一個數(shù)據(jù)比特周期即20 ms的擴頻碼分為20個1 ms,即20個候選比特跳變位置,比較候選位置是否發(fā)生比特跳變,如果在該候選位置處發(fā)生比特跳變,則直方圖計數(shù)加1。經(jīng)過若干個數(shù)據(jù)比特周期的檢測,如果直方圖的峰值超過門限值,則判斷該候選位置為比特跳變位置。目前的研究成果主要是基于這種方法提高其在低信噪比條件下的性能[5-7],而在高動態(tài)條件下,GNSS信號的頻率會產(chǎn)生實時的變化,其變化的規(guī)律和載體運動的狀態(tài)有關(guān),依照傳統(tǒng)比特同步方法和相關(guān)的改進方法,GNSS接收機無法進行穩(wěn)定而可靠地進行比特同步。

圖1 傳統(tǒng)比特同步方法Fig.1 Conventional bit synchronization method

新一代的GNSS民用信號普遍添加了正交或同向的導(dǎo)頻信號,并將導(dǎo)頻信號和數(shù)據(jù)信號進行了分離[8]:全球定位系統(tǒng)(global positioning system, GPS)的L1C、L2C、L5頻點,伽利略導(dǎo)航系統(tǒng)(Galileo satellite navigation system)的E1、E5頻點均采用這導(dǎo)頻信號與數(shù)據(jù)信號分離并同步播發(fā)的形式。我國北斗導(dǎo)航系統(tǒng)的B1C、B2頻點也將采用這種結(jié)構(gòu)[4,6]。導(dǎo)頻信號和數(shù)據(jù)信號的同時出現(xiàn)改變了GNSS接收機的結(jié)構(gòu),同時也使得GNSS信號處理的方式更加多樣化:導(dǎo)頻信號和數(shù)據(jù)信號可以以多種方式實現(xiàn)聯(lián)合捕獲[9],對導(dǎo)頻信號和數(shù)據(jù)信號可以進行聯(lián)合跟蹤[10]。由于導(dǎo)頻信號沒有調(diào)制導(dǎo)航數(shù)據(jù),因而不存在比特跳變,因而可以進行任意時間長度的相干累加,GNSS接收機對信號的檢測方式更加靈活,對高動態(tài)和弱信號的檢測更加靈敏[11]。在研究新一代GNSS信號時發(fā)現(xiàn):導(dǎo)頻信號不存在比特跳變的特性不僅僅能夠延長對GNSS信號的相干累加時間,提高GNSS接收機的捕獲靈敏度,這一特性同時可以用于GNSS數(shù)據(jù)信號的比特同步,尤其是在高動態(tài)條件下這一特性顯得尤為重要。參考文獻[12]提出并總結(jié)了“差分相干累加”的比特同步方法,然而該方法在高動態(tài)條件下不能夠克服多普勒速率的影響。結(jié)合新一代GNSS信號的結(jié)構(gòu),提出一種導(dǎo)頻輔助下的比特同步方法:將當前時刻的數(shù)據(jù)信號相干累加結(jié)果和導(dǎo)頻信號相干累加的結(jié)果進行共軛相乘以復(fù)現(xiàn)當前的導(dǎo)航數(shù)據(jù)比特,并對結(jié)果進行相干累加和非相干累加,對可能的數(shù)據(jù)比特跳變位置處的非相干累加結(jié)果進行比較,判斷比特跳變的實際位置。

1 導(dǎo)頻輔助下的比特同步方法

新一代GNSS信號一般按照數(shù)據(jù)信號和導(dǎo)頻信號功率1∶1或1∶3正交或同相調(diào)制而成,射頻信號經(jīng)過下變頻成為中頻信號后,對中頻信號經(jīng)過捕獲后進行比特同步和跟蹤等后續(xù)處理。

1.1 高動態(tài)條件下的中頻信號捕獲模型

以新體制北斗B2頻點為例分析,B2頻點的參數(shù)如表1所示,整個北斗B2頻點屬于兩路交替二進制偏移載波(alternate binary offset carrier, AltBOC)[13],B2頻點可以分為B2a和B2b兩個信號,每一個信號可以看作正交相移鍵控(quadrature phase shift keying,QPSK)調(diào)制[14]。B2a以功率比為1∶1在QPSK兩個正交的相位上調(diào)制數(shù)據(jù)信號和導(dǎo)頻信號,其擴頻碼長度為10 230 chip、速率為10.23 Mcps、周期為1 ms。

表1 北斗B2頻點調(diào)制方式

當接收機對接收到的B2a信號下變頻,得到含有數(shù)據(jù)分量和導(dǎo)頻分量的中頻信號為

rIF(n)=AndncD(n-τ)cos(2π(fIF+fd)n+φn)+

AncP(n-τ)sin(2π(fIF+fd)n+φn)+ηIF(n)

(1)

式中,An為當前時刻信號幅度;dn=±1為當前時刻數(shù)據(jù)比特符號;τ為信號時延;cD(n)為數(shù)據(jù)通道擴頻碼和副載波乘積;fIF為中頻載波頻率;fd為中頻載波的多普勒頻偏;φn為載波相位;cP(n-τ)為導(dǎo)頻通道擴頻碼和副載波乘積;ηIF(n)是均值為0的高斯白噪聲。對接收信號的數(shù)據(jù)信號進行捕獲,設(shè)相干累加時間為Tc=1 ms,忽略自相關(guān)函數(shù)和頻率偏移對相關(guān)累加值幅度的影響,捕獲后相關(guān)器輸出的1 ms相干累加值為

(2)

(3)

式中,f0為初始多普勒頻偏,由接收機載體相對于衛(wèi)星的速度引起;α為多普勒頻偏速率,由接收機載體相對于衛(wèi)星的加速度引起,忽略自相關(guān)函數(shù)和頻率偏移對相關(guān)累加值幅度的影響,假設(shè)An=A恒定,R(Δτ)≈1,將式(3)代入式(2),得到

(4)

同理,在相同條件下,對接收信號的導(dǎo)頻信號進行相同處理,相關(guān)器輸出的1 ms相干累加值為

(5)

1.2 導(dǎo)頻輔助下的比特同步方法

由式(4)可以看出,由載體加速度多普勒頻偏速率α會導(dǎo)致多普勒頻偏時刻變化,另外接收機載體運動狀態(tài)往往具有不確定性,尤其是載體具有較大加速度情況下,此時α較大,比特同步的性能會受到嚴重影響。不能僅考慮消除初始速度引起的多普勒初始頻偏f0對比特同步造成的影響[12]。注意到導(dǎo)頻信號沒有調(diào)制導(dǎo)航電文,提出利用導(dǎo)頻信號1 ms相干累加值與數(shù)據(jù)信號1 ms的相干累加值進行實時差分的方法,即

(6)

共軛相乘的結(jié)果中第一項為信號項,其余3項為噪聲項?？梢姷谝豁椥盘栱椡耆皇茴l率偏移影響,說明共軛相乘的結(jié)果完全消除了相干累加值中頻率偏移的影響,反映了當前時刻數(shù)據(jù)比特符號dn的值。設(shè)m為n之后20個可能的比特跳變位置中的某一個候選位置,分別以每一個候選位置為起點,按照式(6)對每一時刻相關(guān)器輸出的1 ms輸出進行處理,并將結(jié)果進行20 ms的相干累加得

(7)

當m為比特跳變位置,dn+m+l(l=0,1,2,…,19)必然同號,但是相隔20 ms的數(shù)據(jù)比特不一定變號,即m+20不一定是下一個比特跳變位置,因而有必要根據(jù)實際情況對N個比特(N×20 ms)相干累加的結(jié)果進行非相干累加得

(8)

按照式(8)經(jīng)過N個比特的非相干累加,求出Rn+m,m=0,1,2,…,19中的最大值,即可估算出數(shù)據(jù)比特跳變的位置

(9)

在實際中由于噪聲功率的不確定性,不能判斷經(jīng)過N個比特的非相干累加能夠準確判斷比特跳變的位置,可以定義檢測量ρ,并根據(jù)ρ是否超過某一判決門限ρthres用來判斷該次非相干累加后是否進行判決。定義檢測量ρ為

(10)

即根據(jù)式(9)檢測的出的最大相干累加結(jié)果與相鄰1個候選位置處的相干累加結(jié)果進行比較,若在某次非相干累加結(jié)束后,檢驗出ρ≥ρthres,則對數(shù)據(jù)比特跳變位置進行判決。所提方法的示意圖如圖2所示。

圖2 所提方法示意圖Fig.2 Schematic diagram of the proposed method

1.3 檢測量概率分布分析

假設(shè)m為比特跳變位置,m+20亦為比特跳變時,按照式(7)進行20 ms的相干累加后即對n之后20個候選位置進行判決(即N=1的情況),注意到式(6)中前兩個噪聲項的二階矩為

(11)

(12)

(13)

式中,K0(·)為第二類零階修正的Bessel級數(shù)。在一定的載噪比條件下,信號的幅度和噪聲功率的關(guān)系為

(14)

由于高斯分布的和仍然是高斯分布,高斯分布項與常量項之和包絡(luò)服從萊斯分布,因而式(7)的包絡(luò)近似為萊斯分布,其概率密度函數(shù)為

(15)

(16)

由于m為比特跳變的位置,m+20亦為比特跳變時,因而dn+m+l(l=0,1,2,…,19)同號,此時

(17)

對于m+m′(m′=1,2,3,…,19)這些不是比特跳變的位置的候選位置,由于m+20位置發(fā)生比特跳變,dn+m+m′+l(l=0,1,2,…,19)不同號,此時

(18)

根據(jù)概率密度函數(shù)(17),對其進行積分,可以求出概率分布函數(shù)

(19)

其中

(20)

稱之為一階Marcum-Q分布函數(shù)[17]。對于任何候選位置,其檢測概率均服從萊斯分布,與候選位置距離實際比特跳變位置的距離md有關(guān)。

1.4 正確檢測概率與錯誤檢測概率分析

為分析所提方法的檢測概率,在相隔20 ms的兩個數(shù)據(jù)比特發(fā)生跳變,且對每一個比特跳變位置進行20 ms時長的相干累加的前提下,對當前檢測的候選位置m進行如下對立假設(shè):

(1) 假設(shè)H0:m是實際比特跳變的位置,m+20為下一個比特跳變的位置;

(2) 假設(shè)H1:m不是實際比特跳變的位置,實際比特跳變的位置為m+md,m+md+20為下一個比特跳變的位置。

在H0條件下判斷m為實際比特跳變位置即為正確檢測的概率,由于rn+m(m=0,1,2,…,19)之間相互獨立,因而計算正確檢測概率可以得到:

Pd=P{rn+m>rn+m+m′,m=1,2,…,19|H0}=

(21)

在H0條件下沒有判斷m為實際比特跳變位置為即漏檢概率:

Pm=1-Pd

(22)

在H1條件下判斷m為實際比特跳變位置即為虛警概率,虛警概率可以計算得到:

Pfa=P{rn+m>rn+m+m′,m=1,2,…,19|H1,md}=

(23)

對于虛警概率Pfa,其大小的不僅取決于載噪比的大小,還與該候選位置距離實際比特跳變位置的距離md有關(guān),md越小則虛警概率越大。

2 性能測試與分析

2.1 工作特征曲線分析

在一定條件下以虛警概率為自變量、正確檢測概率為因變量所繪制的曲線稱為接收機的工作特征曲線,工作特征曲線能夠用來評價接收機在一定條件下的工作性能。令接收機處于相對靜止的情況下,載噪比C/N0分別等于35 dBHz、40 dBHz和45 dBHz,對N=1個即20 ms的數(shù)據(jù)進行觀測,進行5 000次實驗并統(tǒng)計在虛警概率Pfa一定的情況下的平均正確檢測概率Pd。虛警概率與候選位置到實際比特跳變位置的距離md有關(guān),折中考慮,仿真中取md=11。

如圖3是本文所提方法的工作特征曲線。由于本文第1.4節(jié)對于所提方法性能的理論分析是近似的,因而當載噪比為35 dBHz時,理論性能與實際仿真得出的性能具有差異。當載噪比上升至40 dBHz、45 dBHz時,理論性能與實際性能的差異已經(jīng)可以忽略不計,說明而理論性能與實際性能的差異隨著載噪比的上升而減小,實驗因而證明了本文對接收機性能的理論分析的合理性。

圖3 工作特征曲線圖Fig.3 Receiver operating characteristic curve

2.2 高動態(tài)條件下的比特同步性能分析

為了驗證所提方法在高動態(tài)條件下的比特同步性能,對相同條件下對傳統(tǒng)比特同步方法、參考文獻[12]中的方法、本文所提方法進行仿真,并對數(shù)據(jù)進行100 ms的觀測后進行判決。仿真考慮的3種動態(tài)情形如表2所示。

表2 仿真參數(shù)設(shè)置

對于情形1,當載噪比為35 dBHz時,其檢測示意圖如圖4。傳統(tǒng)直方圖法的檢測示意圖如圖4所示。由于噪聲和初始多普勒頻偏的干擾,傳統(tǒng)方法經(jīng)過N=5即100 ms的檢測與統(tǒng)計,無法根據(jù)檢測的結(jié)果從候選比特跳變位置中判斷出比特跳變的實際位置,而所提方法如圖5所示。經(jīng)過5個導(dǎo)航數(shù)據(jù)比特的檢測,其檢測峰值處于第11個候選位置。其余候選位置隨著遠離第11個候選位置其檢測結(jié)果也逐漸遞減,其檢測結(jié)果具有一定的規(guī)律性。

圖4 傳統(tǒng)方法檢測結(jié)果Fig.4 Detection result of traditional method

圖5 所提方法檢測結(jié)果Fig.5 Detection result of proposed method

由圖6～圖8的仿真結(jié)果可以看出,在高動態(tài)環(huán)境下,初始多普勒頻和多普勒頻偏速率均會嚴重影響傳統(tǒng)比特同步的性能,在3種情景下,傳統(tǒng)方法均處于失效狀態(tài),其比特同步的正確檢測概率也不隨信噪比明顯變化,且低于0.2,此時傳統(tǒng)檢測方法基本失效。由仿真的結(jié)果同時可以看出,參考文獻[12]中的方法(對照方法)具有在高動態(tài)條件下比特同步的能力,當多普勒頻偏速率增大時,其性能會受到一定的影響。所提的導(dǎo)頻輔助下的比特同步方法不僅能夠在3種動態(tài)情形下正常工作,且在3種情景下變化趨勢和性能也基本相同:當載噪比增大,正確檢測概率趨近于100%,載噪比由48 dBHz降低到24 dBHz的條件下,所提算法的正確檢測概率大約下降0.40～0.43。說明載體的速度和加速度引起的初始多普勒頻和多普勒頻偏速率的改變對所提方法的影響較低,影響所提方法的主要因素為載噪比的大小,所提方法在高動態(tài)條件具有較為可靠的比特同步性能。

圖6 情景1方法性能Fig.6 Method performance in scene 1

圖7 情景2方法性能Fig.7 Method performance in scene 2

圖8 情景3方法性能Fig.8 Method performance in scene 3

2.3 檢測時間分析

文中第2.1節(jié)和第2.2節(jié)對信號進行100 ms觀測后進行判決。實際情況中不知道載噪比的大小,需要按照式(10)進行若干數(shù)據(jù)比特的非相干累加后,按照判決門限進行判決。改變判決門限ρthres的大小,同時改變載噪比,統(tǒng)計此時所提算法需要的判決時間,結(jié)果如圖9所示:判決時間與載噪比和所設(shè)判決門限有關(guān),當載噪比較低時而判決門限設(shè)置較高時,檢測時間會驟增到800 ms以上。實際中應(yīng)根據(jù)載噪比的大小合理設(shè)置判決門限,以免造成漏檢或虛警。

圖9 檢測時間仿真結(jié)果Fig.9 Detection time simulation result

3 結(jié) 論

提出一種導(dǎo)頻輔助的比特同步方法。由理論推導(dǎo)和仿真實驗可以證明,該方法能夠非常有效的降低高動態(tài)條件下的初始多普勒頻偏和多普勒頻偏速率對接收機造成的影響,其在高動態(tài)條件下的比特同步的性能是傳統(tǒng)比特同步方法無法達到的。相比于傳統(tǒng)方法,其運算復(fù)雜度會有一定程度的增加,但該方法具有工程利用的價值,且易于實現(xiàn)。該方法可以推廣到數(shù)據(jù)信號和導(dǎo)頻信號功率按照1∶1正交調(diào)制的GNSS頻點(北斗B1C頻點),也可以應(yīng)用到數(shù)據(jù)信號和導(dǎo)頻信號同相調(diào)制的GNSS頻點。

[1] 趙琳,高帥和,丁繼成.基于FFT的高動態(tài)GPS信號捕獲方法優(yōu)化[J].系統(tǒng)工程與電子技術(shù),2011,33(1):151-156.

ZHAO L, GAO X H, DING J C.Improved acquisition for high dynamic GPS signals based on FFT[J]. Systems Engineering and Electronics, 2011, 33(1):151-156.

[2] 馬永奎, 張一, 張中兆,等. 改進的高動態(tài)高靈敏度GPS信號捕獲算法[J]. 系統(tǒng)工程與電子技術(shù), 2009, 31(2):265-269.

MA Y K, ZHANG Y, ZHANG Z Z, et al. Modified method of high dynamic and high sensitivity GPS signal acquisition[J]. Systems Engineering and Electronics, 2009, 31(2):265-269.

[3] MIAO Z Y, Lü Y L, XU D J, et al. Analysis of a variational Bayesian adaptive cubature Kalman filter tracking loop for high dynamic conditions[J]. GPS Solutions, 2016, 21(1):1-12.

[4] BETZ J W. Engineering satellite-based navigation and timing: Global navigation satellite systems, signals, and receivers[M].New Jersey: Wiley, 2016:157-405.

[5] LEI D, LU W, ZHANG Y, et al. Balance differential coherent bit synchronization algorithm for GNSS receiver[J]. IEICE Trans.on Communications, 2015, E98.B(6):1133-1140.

[6] SOKOLOVIC V S, MARKOVIC G B, DIKIC G D. Improvement of signal synchronization in GPS software receiver[J]. Military Technical Courier, 2015, 63(4):11-24.

[7] QI J, SONG P. Research on the bit synchronization algorithm and multiplexing technology in GNSS receiver[J].China Communications(English version), 2014,11(s2):30-36.

[8] QUAN Y, LAU L, ROBERTS G W, et al. Measurement signal quality assessment on all available and new signals of multi-GNSS (GPS, GLONASS, Galileo, BDS, and QZSS) with real data[J]. Journal of Navigation, 2015, 69(2):313-334.

[9] 田豐,唐小妹,歐鋼.數(shù)據(jù)和導(dǎo)頻通道的功率加權(quán)聯(lián)合捕獲算法[J].國防科技大學學報,2016,38(2):112-116.

TIAN F,TANG X M,OU G. Power-weighted combining acquisition algorithm with data and pilot signals[J].Journal of National University of Defense Technology,2016,38(2):112-116.

[10] WARD P W, FUCHSER T D. Stability criteria for GNSS receiver tracking loops[J]. Navigation, 2015, 61(4):293-309.

[11] ELMOWAFY A. Pilot evaluation of integrating GLONASS, Galileo and BeiDou with GPS in Araim[J]. Artificial Satellites, 2016, 51(1):31-44.

[12] 李新山.基于差分相干的微弱GPS信號捕獲與比特同步研[D]. 成都:電子科技大學, 2014.

LI X S. Research on weak GPS signal acqucistion and bit synchronization based on differential coherent[D]. Chengdu: University of Electronic Science and Technology of China,2014.

[13] HUANG X M, ZHU X W, TANG X M, et al. Asymmetric AltBOC modulation and its generalised form for BeiDou B2 signals[J]. Electronics Letters, 2016, 52(2):146-148.

[14] ZHU Y, CUI X, LU M. Dual binary phase-shift keying tracking method for Galileo E5 AltBOC(15,10) signal and its thermal noise performance[J]. IET Radar, Sonar & Navigation, 2015, 9(6):669-680.

[15] YANG C, NGUYEN T, BLASCH E, et al. Post-correlation semi-coherent integration for high-dynamic and weak GPS signal acquisition[C]∥Proc.of the IEEE/ION position, location and navigation symposium, 2008:1341-1349.

[16] PAGOT J B, THEVENON P, JULIEN O. Signal quality monitoring for new GNSS signals[C]∥Proc.of the 29th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS+ 2016), 2016:1750-1763.

[17] SIMON M K. Probability distributions involving Gaussian random variables[M]. US: Springer, 2006:1-16.