李利軍,黃際彥,劉征

(1.中國鐵建電氣化局集團(tuán)第二工程有限公司,山西 太原 030023;2.電子科技大學(xué),四川 成都 610054)

基于循環(huán)相關(guān)函數(shù)的北斗衛(wèi)星導(dǎo)航抗干擾算法

李利軍1,黃際彥2,劉征2

(1.中國鐵建電氣化局集團(tuán)第二工程有限公司,山西 太原 030023;2.電子科技大學(xué),四川 成都 610054)

當(dāng)前國內(nèi)外衛(wèi)星導(dǎo)航信號(hào)空域抗干擾技術(shù),或者需要信號(hào)到達(dá)角(DOA)信息,或者需要參考信號(hào)先驗(yàn)信息,或者會(huì)對(duì)導(dǎo)航信號(hào)和干擾同時(shí)抑制,無法形成最優(yōu)的輸出信干噪比(SINR),或者僅對(duì)白噪聲干擾有效,無法抑制人造干擾。針對(duì)上述問題,本文提出了基于循環(huán)平穩(wěn)特性的北斗抗干擾盲波束形成新算法。本文首先分析了人造干擾和北斗信號(hào)的循環(huán)平穩(wěn)特性,基于此分析構(gòu)造了新的盲波束形成算法的代價(jià)函數(shù),并采用了拉格朗日方法推導(dǎo)了算法的解。該方法無需DOA和參考信號(hào)先驗(yàn)信息,可以形成陣列增益,得到最優(yōu)輸出SINR,可以同時(shí)抑制寬帶白噪聲干擾和人造調(diào)制干擾。理論分析和仿真結(jié)果證明了該方法的有效性。

北斗衛(wèi)星導(dǎo)航系統(tǒng);人造干擾;循環(huán)平穩(wěn)特性;輸出信干噪比

0 引 言

北斗衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)與美國的GPS、俄羅斯的GLONASS、歐盟的Galileo系統(tǒng)組成了全球四大衛(wèi)星導(dǎo)航系統(tǒng),是中國獨(dú)立研發(fā)、自主運(yùn)行的衛(wèi)星導(dǎo)航系統(tǒng)[1-2]。北斗衛(wèi)星導(dǎo)航系統(tǒng)作為一種無線電導(dǎo)航定位系統(tǒng),在設(shè)計(jì)之初并未重點(diǎn)考慮抗干擾性能,不可避免的具有易受電磁干擾的特性。導(dǎo)航衛(wèi)星信號(hào)的發(fā)射功率和接收信噪比低,接收機(jī)跟蹤環(huán)路易失鎖等特點(diǎn)使得干擾一方很容易對(duì)衛(wèi)星信號(hào)進(jìn)行有意的針對(duì)性的欺騙和干擾。研究表明,正常接收條件下衛(wèi)星接收機(jī)天線的信噪比低于-20 dB.因此,使用一個(gè)很小功率的干擾源即可大范圍內(nèi)的使得衛(wèi)星接收機(jī)不能正常工作。由此可見,北斗衛(wèi)星導(dǎo)航系統(tǒng)的脆弱性會(huì)對(duì)經(jīng)濟(jì)與軍事利益產(chǎn)生巨大的損害,發(fā)展抗干擾技術(shù)的意義重大[2-3]。

就衛(wèi)星接收機(jī)抗干擾技術(shù)而言,主要分為時(shí)-頻域?yàn)V波技術(shù)[4-5]和空域?yàn)V波技術(shù)[6-8]。時(shí)-頻域?yàn)V波技術(shù)只能處理窄帶干擾,對(duì)寬帶干擾、欺騙式干擾、多徑和偽衛(wèi)星的遠(yuǎn)近效應(yīng)效果不佳,因此發(fā)展了空時(shí)域聯(lián)合濾波技術(shù)?；诓ㄊ纬傻目沼?yàn)V波技術(shù)是目前導(dǎo)航信號(hào)抗干擾技術(shù)的主要研究方向。傳統(tǒng)的導(dǎo)航信號(hào)空域抗干擾技術(shù)主要分為三類:一是MVDR類算法[6]。該類算法需要知道衛(wèi)星信號(hào)的DOA信息,該信息在實(shí)際環(huán)境中較難獲得,對(duì)陣列誤差敏感;二是MMSE、LMS和RLS等算法[7]。該類算法依靠鎖相環(huán)捕獲的導(dǎo)航信號(hào)作為參考信號(hào),然而在干擾環(huán)境中,鎖相環(huán)無法捕獲到參考信號(hào);三是功率倒置算法[8-9],該算法無需DOA和參考信號(hào)信息,但只能抑制干擾,無法形成信號(hào)陣列增益。

為了在無先驗(yàn)信息情況下,抑制干擾,同時(shí)得到陣列增益,研究人員發(fā)展了基于衛(wèi)星信號(hào)周期平穩(wěn)特性的OCAB盲波束形成算法[10]。盡管該算法能夠有效抑制寬帶白噪聲干擾,然而其未考慮人造調(diào)制干擾與衛(wèi)星信號(hào)具有相同的周期平穩(wěn)特性,在人造干擾環(huán)境中,性能不佳。

針對(duì)上述問題,本文提出了基于循環(huán)平穩(wěn)特性的北斗抗干擾盲波束形成算法。該方法無需DOA和參考信號(hào)先驗(yàn)信息,可以形成陣列增益,形成最優(yōu)輸出SINR,可以同時(shí)抑制寬帶白噪聲干擾和人造干擾。本文首先分析了人造干擾和北斗衛(wèi)星信號(hào)的循環(huán)平穩(wěn)特性,指出OCAB算法中的周期平穩(wěn)特性是循環(huán)平穩(wěn)特性的一個(gè)特例,并證明了人造干擾和衛(wèi)星信號(hào)具備相同的周期平穩(wěn)特性和不同的循環(huán)頻率。這也是OCAB算法在人造干擾條件下失效的主要原因?；谘h(huán)平穩(wěn)特性的分析,本文構(gòu)造了新的盲波束形成算法的代價(jià)函數(shù),并采用了拉格朗日方法推導(dǎo)了算法的解。仿真結(jié)果驗(yàn)證了本文算法的有效性。

1 基于多天線陣列的北斗信號(hào)模型

北斗二代衛(wèi)星導(dǎo)航系統(tǒng)采用碼分多址(CDMA)的多載波衛(wèi)星導(dǎo)航系統(tǒng),提供了多種的導(dǎo)航信號(hào)以滿足不同用戶的需求。北斗二代信號(hào)包括B1、B2和B3三個(gè)頻點(diǎn)的信號(hào),每個(gè)頻點(diǎn)信號(hào)都是由I、Q兩個(gè)支路的“測距碼+導(dǎo)航電文+載波”構(gòu)成。測距碼分I支路的普通測距碼(C碼)和Q支路精密測距碼(P碼)。本文以B3頻點(diǎn)的I支路為例說明,本文算法可直接擴(kuò)展到其他頻點(diǎn)的I、Q支路。B3頻點(diǎn)I支路輸出信號(hào)為:

S(t)=AICIDI(t)cos(2πf0t+φ) ,

(1)

式中:下角標(biāo)I表示I支路;A為信號(hào)振幅;C為測距碼;D為測距碼上調(diào)制的數(shù)據(jù)碼;f0為載波頻率;Φ為載波初相。空域?yàn)V波技術(shù)采用多天線陣列進(jìn)行衛(wèi)星導(dǎo)航抗干擾處理。假設(shè)共有M個(gè)陣元,第一個(gè)陣元位置為(0,0,0),第i個(gè)陣元位置為(xi,yi,zi),衛(wèi)星信號(hào)S(n)的方位角和俯仰角為(θd,φd),共有N個(gè)干擾,第j個(gè)干擾Ij(n)的來波方向?yàn)?θj,φj),則天線陣列接收信號(hào)可建模為

(2)

式中:n(n)為測量噪聲;ad和aj分別為衛(wèi)星信號(hào)和干擾的陣列流型。

(3)

式中,c為電波傳播速度。空域波束形成的目的就是找到每個(gè)陣元的權(quán)值,在陣列輸出端進(jìn)行加權(quán)疊加,從而在期望信號(hào)方向形成波束,在干擾處形成零陷,其陣列輸出信號(hào)為

y(n)=wHx(n),

(4)

2 基于周期平穩(wěn)的OCAB盲波束形成算法

為了在無先驗(yàn)信息情況下,抑制干擾,同時(shí)得到陣列增益,研究人員發(fā)展了基于衛(wèi)星信號(hào)周期平穩(wěn)特性的OCAB盲波束形成算法,并運(yùn)用在GPS衛(wèi)星系統(tǒng)抗干擾中。GPS和北斗信號(hào)結(jié)構(gòu)形式類似,由式(1)可知,衛(wèi)星信號(hào)由數(shù)據(jù)碼、擴(kuò)頻碼和載波組成,圖1示出了C碼的結(jié)構(gòu)圖。

圖1 衛(wèi)星信號(hào)結(jié)構(gòu)

由圖1可知,一個(gè)衛(wèi)星信號(hào)的數(shù)據(jù)碼包含多個(gè)周期的C碼(測距碼)。設(shè)當(dāng)前處理碼片為S(n),參考碼片為S(n-L),其中L為C碼周期的整數(shù)倍。如果當(dāng)前處理碼片與參考碼片同屬于一個(gè)數(shù)據(jù)碼,S(n)與S(n-L)將完全是相同的信號(hào),這就是C碼的周期平穩(wěn)特性。OCAB算法首先構(gòu)造它們的互協(xié)方差矩陣[10]:

Rxxl=E{x(n)x(n-L)H}

(5)

其中,Rs=E[S(n)S(n-L)H];Rj=E[Ij(n)Ij(n-L)H];Rn=E[n(n)n(n-L)H]=0.

明顯的,由于L為C碼周期的整數(shù)倍,Rs是相干積累可以得到最大互相關(guān)值,而此時(shí)Rn=0.OCAB算法利用上述周期平穩(wěn)特性構(gòu)造代價(jià)函數(shù)[10]:

(6)

式中:w為波束形成要求解的陣列權(quán)向量;Rxx=E{x(n)x(n)H}。該代價(jià)函數(shù)的目的是利用GPS信號(hào)的周期平穩(wěn)特性得到衛(wèi)星信號(hào)的最大增益,同時(shí)抑制干擾Ij(n)和噪聲n(n)。

3 基于循環(huán)平穩(wěn)特性的北斗抗干擾盲波束形成新算法

3.1 北斗信號(hào)循環(huán)平穩(wěn)特性分析

盡管OCAB算法利用了衛(wèi)星信號(hào)的周期平穩(wěn)特性對(duì)衛(wèi)星信號(hào)進(jìn)行了相干積累,得到了最大的陣列增益,同時(shí)又最小化了噪聲,然而其未考慮大多數(shù)人造干擾的周期平穩(wěn)特性,在存在人造干擾情況下,性能惡化。若干擾是寬帶白噪聲干擾,此時(shí)干擾Ij(n)與Ij(n-L)相互獨(dú)立,且均值為0,則

Rj=E[Ij(n)Ij(n-L)H]

=E[Ij(n)]E[Ij(n-L)H]=0.

(7)

此時(shí),OCAB算法性能最優(yōu)。若存在人造干擾時(shí),Rj≠0,即意味著人造干擾與衛(wèi)星信號(hào)具備相同的周期平穩(wěn)特性,采用基于公式(6)的OCAB算法將會(huì)同時(shí)將人造干擾最大化,從而導(dǎo)致算法失效。下面以調(diào)幅干擾(AM)的循環(huán)平穩(wěn)特性為例,說明上述結(jié)論。考慮實(shí)隨機(jī)AM干擾:

x(t) =a(t)cos(2πf0t+θ)

(8)

信號(hào)的循環(huán)平穩(wěn)特性可由循環(huán)平穩(wěn)相關(guān)函數(shù)描述:

(9)

式中:τ為循環(huán)周期;β為循環(huán)頻率。對(duì)比式(5)和式(9)可知OCAB算法的周期平穩(wěn)特性是循環(huán)平穩(wěn)特性當(dāng)循環(huán)頻率β為零時(shí)的特例。

由于式(8)中a(t)為實(shí)平穩(wěn)信號(hào),則有[11]:

β≠0

(10)

將式(8)和式(10)代入式(9),可得實(shí)隨機(jī)AM干擾的循環(huán)相關(guān)函數(shù):

(11)

從上式中可知,對(duì)于AM干擾來說,當(dāng)循環(huán)頻率β=0時(shí),式(5)中的Rj≠0.對(duì)于其他形式的調(diào)制干擾來說,采用相同的分析方法可以得到相似的結(jié)論。這也說明了人造干擾與衛(wèi)星信號(hào)具備相同的周期平穩(wěn)特性,采用基于公式(5)和(6)的OCAB算法將會(huì)同時(shí)將人造干擾最大化,從而導(dǎo)致算法失效。

采用相同的方法也可以證明衛(wèi)星信號(hào)也具備循環(huán)平穩(wěn)特性,其循環(huán)相關(guān)函數(shù)為

(12)

式中:ε為C碼的碼率,對(duì)于北斗二代B3頻點(diǎn)來說ε=10.23 MHz.圖2為北斗二代衛(wèi)星基帶信號(hào)與其自身的循環(huán)相關(guān)函數(shù)圖和單頻干擾與其自身的循環(huán)相關(guān)函數(shù)圖。圖中單頻干擾定義為cos(2πf0t),f0=3.5 MHz.

圖2 北斗信號(hào)與單頻干擾環(huán)相關(guān)函數(shù)

從圖中可知,衛(wèi)星信號(hào)與單頻干擾的周期平穩(wěn)特性(β=0)一致,基于周期平穩(wěn)特性的OCAB對(duì)人造干擾無效。但是衛(wèi)星信號(hào)與單頻干擾循環(huán)頻率不一樣,因此采用循環(huán)相關(guān)函數(shù)可以抑制大部分人造干擾。

3.2 基于循環(huán)平穩(wěn)特性的北斗抗干擾盲波束形成算法

基于上述分析,本文提出了基于循環(huán)平穩(wěn)特性的北斗抗干擾盲波束形成新算法。首先重構(gòu)代價(jià)函數(shù)

(13)

采用拉格朗日方法,上述代價(jià)函數(shù)變?yōu)?/p>

J(w,c)=wHRxxw+cHRuuc-

(14)

其中,γ為拉格朗日參數(shù)。令wJ(w,c)=0和cJ(w,c)=0,最小化J(w,c),得:

(15)

由上式可解得:

(16)

將式(16)代入式(15),可得:

(17)

(18)

由上式可知,加權(quán)矩陣w的解為以下矩陣最大奇異值所對(duì)應(yīng)的奇異矢量:

(19)

將加權(quán)矩陣w代入y(k)=wHx(k)中,即可將衛(wèi)星信號(hào)S(k)提取出來。由上述代價(jià)函數(shù)和求解過程中可知,盲波束形成技術(shù)實(shí)際上就是利用了衛(wèi)星信號(hào)在L和β處的循環(huán)對(duì)稱相關(guān)矩陣Rss(β,L)不為0,而干擾信號(hào)和噪聲在τ和α處的循環(huán)對(duì)稱相關(guān)矩陣為0的特性，如圖2所示,最大化衛(wèi)星信號(hào),而最小化干擾信號(hào)和噪聲。

在求解濾波器權(quán)值w時(shí),需要求解協(xié)方差陣Rxx、Rxu和Ruu,由于實(shí)際應(yīng)用和仿真中,它們可以采用多個(gè)采樣點(diǎn)信號(hào)的估計(jì)矩陣,對(duì)于具有N個(gè)采樣點(diǎn)信號(hào)組成的數(shù)據(jù)塊xN和參考數(shù)據(jù)塊uN有:

(20)

由于數(shù)據(jù)碼的存在,數(shù)據(jù)塊可能受數(shù)據(jù)碼的跳變影響。當(dāng)輸入信號(hào)數(shù)據(jù)塊或者參考數(shù)據(jù)塊被數(shù)據(jù)碼跳變點(diǎn)隔開,則在協(xié)方差陣估計(jì)過程中可能會(huì)使部分?jǐn)?shù)據(jù)被抵消掉,引起估計(jì)值的不準(zhǔn)確,對(duì)于該問題,可采用文獻(xiàn)[10]中的滑動(dòng)平均方法解決。

4 仿真結(jié)果

參考北斗二代B3頻點(diǎn)參數(shù)對(duì)仿真環(huán)境進(jìn)行設(shè)置,參數(shù)如表1所示。

表1 北斗二代B3頻點(diǎn)仿真參數(shù)

仿真中,MVDR算法[6]、OCAB算法[10]及功率倒置算法[8]都將與本文算法進(jìn)行比較。OCAB算法、功率倒置算法和本文算法都屬于盲波束形成算法,無需知道衛(wèi)星信號(hào)到達(dá)角的先驗(yàn)信息。而MVDR算法需要知道衛(wèi)星信號(hào)到達(dá)角的先驗(yàn)信息,當(dāng)該先驗(yàn)信息無誤差時(shí)性能最優(yōu),而存在誤差時(shí),性能惡化。多天線陣列采用9天線圓陣。衛(wèi)星到達(dá)方位角和俯仰角為(240°,50°),干擾方向?yàn)?40°,75°)。干擾類型分別建模為寬帶白噪聲干擾、寬帶掃頻干擾和單頻干擾。圖3示出了不同寬帶白噪聲干擾下的輸出信干噪比(SINR),從圖中可以看出,OCAB、本文算法和MVDR算法的SINR高于功率倒置算法。這是因?yàn)楣β实怪盟惴ㄖ荒芤种茝?qiáng)干擾,無法形成陣列增益。另外OCAB、本文算法和MVDR算法的性能一致,說明了OCAB和本文算法在寬帶白噪聲情況下,性能已經(jīng)達(dá)到最優(yōu)。這與式(7)的分析結(jié)果一致。圖4和圖5示出了在人造干擾(單頻干擾和掃頻干擾)下的四種算法的SINR對(duì)比。從圖中可看出本文算法和MVDR算法SINR一致,優(yōu)于其他兩種算法,即本文算法性能已達(dá)最優(yōu)。功率倒置算法能夠有效抑制干擾,但是無法形成陣列增益,而OCAB算法在人造干擾情況下,性能嚴(yán)重惡化,這是由于人造干擾與衛(wèi)星信號(hào)具備完全一致的周期平穩(wěn)特性,基于式(5)和式(6)代價(jià)函數(shù)的OCAB算法會(huì)同時(shí)放大衛(wèi)星信號(hào)及人造干擾。而本文算法考慮了衛(wèi)星信號(hào)與干擾信號(hào)具備不同的循環(huán)頻率這一特性,在人造干擾情況下,性能也達(dá)到了最優(yōu)。

圖3 寬帶白噪聲干擾下的輸出信干噪比

圖4 單頻干擾下的輸出信干噪比

圖5 掃頻干擾下的輸出信干噪比

5 結(jié)束語

基于北斗衛(wèi)星信號(hào)的循環(huán)平穩(wěn)特性,本文提出了一種新的北斗衛(wèi)星抗干擾盲波束形成算法。該算法利用衛(wèi)星信號(hào)與人造干擾之間不同的循環(huán)頻率抑制干擾,無需DOA和參考信號(hào)先驗(yàn)信息,可以形成陣列增益,得到最優(yōu)輸出SINR,能夠同時(shí)抑制寬帶白噪聲干擾和人造干擾。仿真結(jié)果驗(yàn)證了本文算法的有效性。本文中雖然使用的是北斗二代B3頻點(diǎn)仿真環(huán)境,但可以推廣到其他頻點(diǎn)和衛(wèi)星導(dǎo)航系統(tǒng)的抗干擾應(yīng)用中,如GPS、GLONASS和Galileo.

[1] 張乃通,張中兆,李英濤. 衛(wèi)星移動(dòng)通信[M]. 北京:電子工業(yè)出版社,2000.

[2] MISRA P, BURKE B, PRATT M. GPS performance in navigation[J]. Proceedings of the IEEE. 1999, 87(1): 65-85.

[3] 王曉君,張偉,杜萌萌. 北斗衛(wèi)星導(dǎo)航系統(tǒng)的應(yīng)用及其抗干擾技術(shù)[J].河北工業(yè)科技, 2014,31(3):234-238.

[4] 徐海生.基于信道編碼的衛(wèi)星導(dǎo)航信號(hào)抗干擾技術(shù)研究[D].西安:西安電子科技大學(xué),2009.

[5] 韓明. 基于混合擴(kuò)頻的衛(wèi)星導(dǎo)航系統(tǒng)抗干擾技術(shù)研究[D]. 西安:西安電子科技大學(xué),2009.

[6] CAPON J. High resolution frequency-wave number spectrum analysis[J]. Proceedings of the IEEE, 1969, 57(8): 1408-1418.

[7] KAY S M. Fundamentals of statistical signal processing: Estimation theory[M]. Englewood Cliffs, NJ: Prentice-Hall, 1993.

[8] COMPTON R T. The power-inversion adaptive array: concept and performance[J]. IEEE Transactions on Aerospace and Electronic Systems, 1979, 15(6):803-814.

[9] PARKINSON B W, SPILKER J J. Global positioning system: Theory and applications[M]. Boston: Artech House, 1996.

[10]SUN W, AMIN M G. A self-coherence anti-jamming GPS receiver[J]. IEEE Transaction on Signal Processing, 2005, 53(10): 3910-3915.

[11]GARDNER W A. Exploitation of spectral redundancy in cyclostationary signals[J].IEEE Signal Processing Magazine, 1991, 8(2): 14-36.

Beamforming Technique in Anti-Jamming Compass Receiver by Exploitation of Cyclostationarity

LI Lijun1,HUANG Jiyan2,LIU Zheng2

(1.ChinaRailwayConstructionEldctrificationBureauGroupCo.Ltd,Taiyuan030023,China; 2.UniversityofElectronicScienceandTechnologyofChina,Chengdu610054,China)

Although many anti-jamming methods in space domain have been proposed for satellite systems in the literature, they require the prior knowledge of DOA and reference signals, or cannot output optimal SINR since those algorithms suppressed both the desired signal and interference, or cannot mitigate the man-made interference. To counter the problems above, a novel blind beamforming method based on cyclostationarity of satellite signal is proposed in this paper for anti-jamming Compass Receiver. First, the cyclostationarities of compass signal and man-made interference are analyzed in the paper. An objective function based on the analysis results is constructed for the proposed method. Finally, the Lagrange method is used to obtain the solution. Compared with the existing methods, the proposed algorithm need not the prior knowledge of DOA and reference signals, can obtain the optimal array gains, and can suppress both the white noise interference and man-made intergerence. Theoretical analysis and simulation results are provided to show the effectiveness of the proposed method.

Compass Navigation Satellite System; man-made interference; cyclostationarity; signal to interference and noise ratio

10.13442/j.gnss.1008-9268.2017.03.009

2017-01-10

國家自然科學(xué)基金項(xiàng)目(批準(zhǔn)號(hào):61201275)； 中國鐵建電氣化局集團(tuán)開放基金項(xiàng)目資助

P228.4

A

1008-9268(2017)03-0042-06

李利軍 (1973-),男,高級(jí)工程師，現(xiàn)從事鐵路四電工程、城市軌道交通機(jī)電設(shè)備安裝工程的施工技術(shù)研究及施工管理。

黃際彥 (1981-),男,博士,副教授,主要研究方向?yàn)闊o線定位、陣列信號(hào)處理、通信中的信號(hào)處理技術(shù)等。

聯(lián)系人: 黃際彥 E-mail: huangjiyan@uestc.edu.cn