徐 浩,曾芳玲

(國防科技大學(xué)電子對抗學(xué)院,合肥 230037)

0 引言

全球定位系統(tǒng)簡稱GPS,提供的定位服務(wù)分為兩種,分別是由C/A碼調(diào)制的信號提供的標(biāo)準(zhǔn)定位服務(wù)(SPS),和由P碼調(diào)制的信號提供的高精度定位服務(wù)(PPS)[1]。C/A碼周期短,捕獲簡單,而P碼的周期長,捕獲困難。

傳統(tǒng)的P碼捕獲方法,是通過C/A碼間接捕獲P碼的方法[2]。但C/A碼易受外來干擾,一旦被干擾便無法再進(jìn)行P碼的捕獲,因此研究P碼的直接捕獲技術(shù)在GPS應(yīng)用特別是軍事應(yīng)用中是十分迫切的。目前,P碼的直接捕獲方法研究主要集中在頻域。主要包括基于FFT的循環(huán)相關(guān)捕獲法、XFAST(擴(kuò)展復(fù)制重疊)、均值法等方法[3]。

文中以上述均值捕獲方法和XFAST方法為基礎(chǔ),提出一種時(shí)頻融合的捕獲算法,其基本思想是對本地序列進(jìn)行擴(kuò)展復(fù)制重疊處理,在不影響相關(guān)度的基礎(chǔ)上擴(kuò)大搜索范圍;并對處理后的本地序列和接收序列進(jìn)行均值處理,以提高捕獲速度,同時(shí)在時(shí)域上對載波多普勒頻移進(jìn)行并行搜索。另外對這種方法進(jìn)行了進(jìn)一步的優(yōu)化,并通過仿真分析證實(shí)了優(yōu)化后的方法在捕獲速度上有了進(jìn)一步的提高。

1 P碼頻域捕獲算法原理

P碼的碼周期很長,相關(guān)捕獲的時(shí)間比較長,難以滿足實(shí)時(shí)性的要求[4]。所以P碼的快速直接捕獲方法主要集中在頻域,主要是采用頻域相乘代替時(shí)域的相關(guān)運(yùn)算,主要的原理是是基于信號的互相關(guān)、線性卷積以及循環(huán)卷積之間的關(guān)系[5]。對于兩個(gè)信號x(n)和h(n),它們之間的信號相關(guān)和信號FFT之間存在如下的關(guān)系:

?h(-k)=

IFFT{FFT[x(k)·FFT[h(-k)]}=

IFFT{FFT[x(k)·FFT*[h(k)]}

根據(jù)上述關(guān)系,現(xiàn)只需要進(jìn)行2次FFT和1次IFFT運(yùn)算,就能得到與信號相關(guān)同樣的結(jié)果,明顯的減少了運(yùn)算量和運(yùn)算時(shí)間。雖然采用FFT捕獲P碼相對于頻域相關(guān)捕獲來說減少了運(yùn)算量,提高了捕獲的速度,但是P碼的長周期使得直接采用FFT捕獲用時(shí)依然較長,因此在FFT的基礎(chǔ)上進(jìn)行一些處理,來進(jìn)一步提高捕獲速度。

1.1 擴(kuò)展復(fù)制重疊捕獲方法

擴(kuò)展復(fù)制重疊算法以FFT變換為基礎(chǔ)[6],其主要思想是將本地偽碼分成M段,每一段的長度包含L個(gè)碼元,然后將這M段子碼段對應(yīng)位置進(jìn)行算術(shù)相加,得到1個(gè)新的重疊的碼段,這個(gè)碼段的長度為L,但包含了原來長為M×L的碼段的信息。圖1以M=4,L=4對主要的原理進(jìn)行了解釋。

圖1 擴(kuò)展復(fù)制原理

根據(jù)上述原理,將本地P碼分成M段,每段對應(yīng)位置相加得到新的重疊碼段,然后進(jìn)行循環(huán)相關(guān),找到這個(gè)新的重疊碼和接收碼的最大相關(guān)結(jié)果G1。若相關(guān)結(jié)果大于門限值,此時(shí)便得到了接收碼在本地重疊偽碼中相對于段首的相對位置m。再將接收信號和每段子碼段在相對m的位置分別進(jìn)行頻域相關(guān)運(yùn)算,得到m個(gè)相關(guān)結(jié)果,最大相關(guān)峰值G2出現(xiàn)碼段記為d,那么偽碼偏移的位置為(d-1)L+m。

擴(kuò)展復(fù)制重疊減少了FFT的點(diǎn)數(shù),減少了運(yùn)算量,相對于直接FFT方法,捕獲時(shí)間提高了很多。但是由于子碼段之間的疊加,增加了背景噪聲,而且可能產(chǎn)生峰值的多值性,因此,重疊子段數(shù)不能無限制的增加。表1給出了該算法的仿真測試結(jié)果,其中接收信號的參數(shù)為信噪比-12 dB,中頻頻率為20.46 MHz,l(l=2L)為每次FFT的點(diǎn)數(shù),M為擴(kuò)展段數(shù);這里G1、G2分別表示上述過程的兩個(gè)相關(guān)比例峰值;T1、T2分別表示1次搜索到和無峰值條件下搜索1 000 000碼片所用時(shí)間。

表1 XFAST捕獲算法處理結(jié)果

經(jīng)統(tǒng)計(jì)分析表明,M值越大,峰值不存在的情況下搜索時(shí)間相對越短,但峰值也越模糊,碼片誤差大,捕獲越不準(zhǔn)確,因此XFAST在提高捕獲速度上也是有限的。

1.2 均值捕獲方法

均值捕獲算法以循環(huán)相關(guān)理論為基礎(chǔ),并且通過分段補(bǔ)零的方法很好的實(shí)現(xiàn)了頻域的碼多普勒頻移的補(bǔ)償[7]。其原理如圖2所示。

圖2 均值算法原理

根據(jù)上述均值算法原理,取L長的接收碼序列,進(jìn)行補(bǔ)L個(gè)零處理,得到長度為2L的新序列,然后對該序列每i個(gè)碼元取均值后得到長為2k的序列G(2L=i×2k);取長為2L的本地樣本,每i個(gè)碼元取均值得到長為2k的序列F(2L=i×2k);對兩個(gè)新序列G和F按頻域捕獲原理,求得2k個(gè)相關(guān)結(jié)果,取最大的相關(guān)峰值G1進(jìn)行門限判斷,若超過門限則得到最大相關(guān)峰值的相對位置m。取出相應(yīng)位置為m處的相應(yīng)碼長的本地偽碼和相應(yīng)位置相應(yīng)長度的接收碼,按頻域捕獲原理求相關(guān)值,并得到最大峰值G2的位置m1,那么得到最終的偽碼偏移為(m-1)i+m1。

均值算法減少了FFT運(yùn)算的點(diǎn)數(shù),提高了處理的速度,但是平均的碼段越長,它的碼相位偏移越大,峰值捕獲的性能越差,因此平均的碼段長度不能過長,這就在時(shí)間域上限制了檢索的范圍。表2給出了相同仿真參數(shù)下的分析結(jié)果,其中i為均值的點(diǎn)數(shù),G1、G2、T1、T2表示的意義與表1相同。

表2 均值捕獲算法處理結(jié)果

經(jīng)統(tǒng)計(jì)分析i值越大,單次運(yùn)行時(shí)間長,峰值不存在的情況下搜索時(shí)間相對較短,但峰值也越模糊,碼片誤差越大,捕獲越不準(zhǔn)確。因此均值點(diǎn)數(shù)受到限制,捕獲速度提高的水平也受到限制。

2 時(shí)頻融合優(yōu)化算法

2.1 時(shí)頻融合算法

XFAST捕獲算法可以擴(kuò)大搜索的范圍,理論上重疊的段數(shù)越大,捕獲速度越快,但段數(shù)過多會(huì)影響捕獲成功率[8]。均值捕獲方法可以提高相關(guān)處理速度,理論上均值的點(diǎn)數(shù)越多捕獲速度越快,但點(diǎn)數(shù)過多會(huì)使峰值模糊,捕獲成功率下降。為此,文中提出一種時(shí)頻融合捕獲法,結(jié)合以上兩種方法的優(yōu)勢,以達(dá)到更優(yōu)更快的捕獲效果。其基本思想是對本地序列進(jìn)行擴(kuò)展復(fù)制重疊處理,其目的是在不影響相關(guān)度的基礎(chǔ)上擴(kuò)大搜索范圍,然后對處理之后的本地信號以及接收信號進(jìn)行均值處理,來提高捕獲速度,同時(shí)在時(shí)域上進(jìn)行載波的多普勒頻移的并行搜索。其捕獲主要過程如圖3所示。

圖3 時(shí)頻融合捕獲流程

該算法的原理可行性分析如下:

對于一個(gè)長為L的接收P碼序列可表示為A={a1,a2,…,aL},對該序列補(bǔ)L個(gè)零,并且每i個(gè)碼元取均值,得到序列:

相應(yīng)的取長為2L×M的本地偽碼樣本記為B={b1,b2,…,b2L×M},將其分成M段,每段長為2L,然后將這M段對應(yīng)位置相加后,每i個(gè)碼元取均值,可得到如下序列:

假定偽碼的偏移量為0,根據(jù)P碼良好的自相關(guān)特性,只有碼元對齊的情況下兩個(gè)碼元{a1b1,a2b2,…,aLbL}的相關(guān)值不為0,其余相關(guān)值均為無窮小量。那么將序列A′和B′進(jìn)行相關(guān)運(yùn)算,也只有偏移量為0的情況下,存在{a1b1,a2b2,…,aLbL}項(xiàng),才存在最大峰值,這和原序列的運(yùn)算是一致的,是尋找碼片偏移的原理。然后根據(jù)展開后每段本地樣本和接收信號的相關(guān)值,找到具體的碼元頻偏段數(shù),從而最終確定碼片偏移位置。

具體的處理過程為:

1)將接收信號中頻處理后,以適當(dāng)?shù)乃俾?、載波頻率進(jìn)行A/D采樣和緩存,生成一個(gè)長為L的樣本:G′={g1,g2,…,gL}(L=i×k)。

2)取L長的接收信號,并對其進(jìn)行補(bǔ)L個(gè)零操作,這樣得到一個(gè)長為2L的樣本,將其記為:G″={g1,g2,…,g2L}(gL+1,gL+2,…,g2L=0),對該樣本每i個(gè)碼元求均值,得到新的2k長的樣本記為:

G=gq[p]{gq=(g(p-1)i+1+g(p-1)i+2+…+gpi)÷i,
p=1,…,2k}，2L=i×2k

3)取一定偏移量的長為2L×M的本地?cái)U(kuò)頻偽碼樣本,將其分成M段,并且將這M段碼元對應(yīng)位置相加,這樣得到一個(gè)長為2L的樣本。然后在對這個(gè)樣本進(jìn)行類似第二步中的均值運(yùn)算處理,得到長為2k的樣本。

4)對G進(jìn)行FFT運(yùn)算,并取共軛;對F進(jìn)行FFT運(yùn)算;將兩者處理后的結(jié)果相乘后進(jìn)行FFT反變換,對結(jié)果取絕對值,即得到相關(guān)結(jié)果。對相關(guān)結(jié)果進(jìn)行采樣并對2k長度的各峰值進(jìn)行比較,若超過門限則繼續(xù)執(zhí)行,否則,返回重新產(chǎn)生后的本地偽碼序列。

5)去除模糊度,主要包括去除重疊模糊度和均值模糊度。其方法是將接收序列和本地序列進(jìn)行均值展開處理,然后進(jìn)行相關(guān)運(yùn)算,找到最大峰值的位置,從而去除均值模糊度;在去除均值模糊度的基礎(chǔ)上,將本地序列進(jìn)行重疊段數(shù)展開,將每一段和接收序列進(jìn)行相關(guān)運(yùn)算,找到最大峰值的段數(shù),從而去除重疊模糊度。具體的過程同第一節(jié)中相同。

2.2 時(shí)頻融合捕獲算法的進(jìn)一步優(yōu)化

時(shí)頻融合捕獲算法結(jié)合均值法和擴(kuò)展復(fù)制重疊兩種方法的優(yōu)點(diǎn),進(jìn)一步提高了捕獲速度和捕獲精度。但在對多普勒頻移的搜索過程中,采用在不同的本地載波偏移下分別進(jìn)行捕獲的驗(yàn)證方法,獲得最小相關(guān)峰值衰減時(shí)的頻移,完成多普勒補(bǔ)償。但是這種方法增加了捕獲的搜索量,尤其在高動(dòng)態(tài)的情況下使捕獲的速度變得很慢,因此需要做進(jìn)一步的優(yōu)化。

圖4 優(yōu)化算法的捕獲過程

根據(jù)時(shí)域復(fù)指數(shù)相乘等于頻域偏移的原理,通過移動(dòng)頻域數(shù)據(jù)序列對多普勒的搜索進(jìn)行補(bǔ)償。該方法不需要在各個(gè)本地載波頻率偏移下分別進(jìn)行FFT變換,簡化了計(jì)算過程,節(jié)省了頻域變換的計(jì)算量,提高了捕獲速度,具體的捕獲過程如圖4所示。

其主要的處理方法是將2.1節(jié)中的接收序列G進(jìn)行FFT變換,然后根據(jù)多普勒搜索的范圍,按照頻率搜索精度,將接收信號的FFT變換序列進(jìn)行N次移位,就可以產(chǎn)生N個(gè)不同的多普勒補(bǔ)償后接收信號的FFT變換序列。將補(bǔ)償后的序列取共軛后同本地處理后的序列進(jìn)行相乘,并對乘積進(jìn)行IFFT變換,最后進(jìn)行相關(guān)峰值的比較判定。

3 P碼直捕算法的仿真分析與比較

利用Matlab平臺(tái),對優(yōu)化后P碼的直接捕獲算法的捕獲結(jié)果進(jìn)行仿真,并對不同捕獲方法多次捕獲過程進(jìn)行統(tǒng)計(jì)分析。仿真參數(shù)設(shè)定為信噪比-12 dB,中頻頻率20.46 MHz,截短碼長為2 048。其優(yōu)化算法的捕獲結(jié)果如圖5所示。

圖5 捕獲的相關(guān)峰值和歸一化峰值

表3、表4、表5是各算法捕獲結(jié)果分析。

表3 時(shí)頻融合算法捕獲結(jié)果

表4 優(yōu)化算法捕獲結(jié)果

表5 不同捕獲方法在同等參數(shù)下的性能比較

根據(jù)捕獲性能的統(tǒng)計(jì)分析,時(shí)頻融合捕獲方法雖然單次運(yùn)行時(shí)間T1相對擴(kuò)展復(fù)制重疊法和均值法較長,但是這種相差很小,只是幾秒的差別。然而從總體的捕獲時(shí)間T2來看,時(shí)頻融合算法較擴(kuò)展復(fù)制重疊法捕獲時(shí)間縮短了一個(gè)數(shù)量級,較均值法也縮短了50%左右,捕獲速度有了很大的提高。從多普勒誤差來看,視頻融合法較擴(kuò)展復(fù)制重疊法誤差要小一個(gè)數(shù)量級,精度有了很大提高。并且通過對優(yōu)化算法和時(shí)頻融合算法的對比,可以發(fā)現(xiàn)優(yōu)化算法的單次運(yùn)行時(shí)間T1和總體捕獲時(shí)間T2較時(shí)頻融合算法有了更進(jìn)一步的提高,取得了優(yōu)化的效果,達(dá)到了快捕的目的。

4 結(jié)束語

文中主要對P碼的直接捕獲算法進(jìn)行研究和分析,提出了一種結(jié)合均值法和擴(kuò)展復(fù)制重疊法兩種方法的時(shí)頻融合的捕獲算法,對其進(jìn)行了Matlab模擬仿真試驗(yàn)。并且對算法進(jìn)行了進(jìn)一步的優(yōu)化,通過模擬實(shí)驗(yàn)結(jié)果的統(tǒng)計(jì)分析,證實(shí)了優(yōu)化后的算法在捕獲速度和精度上都有了很大的提高,滿足了現(xiàn)在P碼直接捕獲的快速要求。