
第31卷第1期 計算機仿真 2014年1月
文章編號:1006—9348(2014)O1—0072-04
衛(wèi)星導(dǎo)航接收機中窄帶干擾抑制算法
藺曉龍,何文濤,徐建華,葉甜春
(中國科學(xué)院微電子研究所,北京100029)
摘要:窄帶干擾能夠嚴(yán)重降低衛(wèi)星導(dǎo)航接收機的性能。為提高衛(wèi)星導(dǎo)航性能,提出一種應(yīng)用于衛(wèi)星導(dǎo)航接收機中的窄帶干
擾抑制算法。首先利用數(shù)字傅里葉變換(DFT)檢測窄帶干擾,傳統(tǒng)DFF實現(xiàn)簡單、便于集成,但計算慢、檢測時間過長。為
解決上述問題,采用最優(yōu)化的多級DFT頻率檢測算法,能顯著減小檢測時間。然后使用級聯(lián)的IIR陷波器抑制窄帶干擾。
最后使用頻率檢測和濾波器更新相結(jié)合的開環(huán)結(jié)構(gòu),用于跟蹤掃頻信號并進行抑制。仿真結(jié)果表明,改進算法能夠有效的
抑制衛(wèi)星導(dǎo)航接收機中的窄帶和掃頻干擾。
關(guān)鍵詞:衛(wèi)星導(dǎo)航接收機;窄帶干擾;數(shù)字傅里葉變換;掃頻干擾
中圖分類號:TN911 文獻標(biāo)識碼:B
Narrow-Band Interference Mitigation Algorithm in GNSS Receiver
LIN Xiao-long,HE Wen—tao,XU Jian-hua,YE Tian-chun
(Institute of Microelectronics of Chinese Academy of Sciences,Beijing 100029,China)
ABSTRACT:Narrow—band interference can severely degrade the performance of the GNSS receivers.To solve this
problem,this article proposed a na ̄owband interference suppression algorithm applied in muhimode GNSS receiver.
Unlike the traditional DFF,which is simple and easy to be integrated,but take too much detection time,an optimized
multi-level DYF ̄equeney detection algorithm which can signiifcantly reduce the detection time was proposed in this
article.After being detected by the multi-level DFT,the narrow-band interference was then mitigated by the cascade
notch filters.Finally,an open loop approach combing the operation of detection and update was introduced to track
and mitigate the sweep interference at the same time.The simulation results show that this algorithm can effectively
suppress narrowband and sweep interferences signals in the GNSS receivers.
KEYWORDS:GNSS receiver;Narrow-band Interference;DFT;Sweep interference
1 引言
GPS、GLONASS、Galileo以及中國的北斗衛(wèi)星導(dǎo)航系統(tǒng)
持續(xù)抑制。
目前有多種窄帶干擾抑制方法,包括:文獻[3][4]使用
F 技術(shù)在頻域?qū)φ瓗Ц蓴_進行檢測和抑制。由于使用
FFT,需要大量緩存數(shù)據(jù),計算量大。文獻[5][6]使用自適
應(yīng)陷波器(ANF)技術(shù),LMS、RLS等多種自適應(yīng)算法被采用,
可使用IIR型或FIR型的濾波器。其中自適應(yīng)算法計算量、
檢測數(shù)目、收斂性需要特別關(guān)注。文獻[7]使用檢測加抑制
的結(jié)構(gòu)對包括窄帶干擾在內(nèi)的多種干擾進行抑制。文獻[8]
使用其它更為復(fù)雜的一些變換與濾波技術(shù),能夠抑制包括多
(BDS)以及其它的增強系統(tǒng)統(tǒng)稱為全球?qū)Ш叫l(wèi)星系統(tǒng)
(GNSS)。地面接收到的GNSS信號強度在一130dBm左右,
遠(yuǎn)在噪聲功率之下 。隨著定位功能越來越多的集成到移
動設(shè)備中,由其它通信設(shè)備泄露到GNSS信號頻段的干擾不
斷增多。這些干擾會嚴(yán)重降低衛(wèi)星導(dǎo)航接收機的性能。其
中對GNSS接收機影響最大,且最常見的干擾類型是窄帶干
擾 。
另外接收機與干擾源之間的相對運動會產(chǎn)生多普勒效
種類型的干擾,考慮了一定的掃頻干擾,但算法及實現(xiàn)復(fù)雜
度都較高,不易于集成。
本文使用檢測加濾波的結(jié)構(gòu)實現(xiàn)衛(wèi)星導(dǎo)航接收機中的
應(yīng),此時即使外界窄帶干擾頻率是固定不變的,GNSS接收機
內(nèi)部,其干擾頻率會隨時間發(fā)生變化,形成掃頻干擾信號。
所以在處理窄帶干擾抑制時需要考慮跟蹤干擾頻率以達(dá)到
基金項目:中國第二代衛(wèi)星導(dǎo)航系統(tǒng)重大專項多模導(dǎo)航型基帶芯片
(2011GFZX03030204)
窄帶干擾抑制。DFT技術(shù)由于其硬件實現(xiàn)簡單,便于集成而
被廣泛采用,但是由于衛(wèi)星接收機中信號檢測帶寬很大,如
2MHz,頻率分辨率又要求很小,如100Hz,傳統(tǒng)DFT計算緩
慢,檢測時間過長而無實際應(yīng)用價值。針對這個問題,本文
提出的多級DFT窄帶干擾頻率檢測算法,使用多等級的檢
收稿日期:2013—04—10
72一
測,從而最小化檢測時間。然后使用級聯(lián)的IIR陷波器陣
列,去除多個輸入通道中的多個窄帶干擾。最后采用DFT檢 波器能去除一個實數(shù)輸入的窄帶干擾,要應(yīng)對多個窄帶干
測加調(diào)整陷波器的開環(huán)方式對掃頻干擾進行抑制。
<b<1表示與IIR陷波器帶寬有關(guān)的參數(shù)。單個二階IIR陷
擾,需要多級級聯(lián)的二階IIR陷波器序列。
2.2 窄帶干擾頻率檢測
2衛(wèi)星導(dǎo)航接收機干擾抑制結(jié)構(gòu)
受多個連續(xù)波干擾的擴頻信號模型:
(£)=5(f)+n(t)+ A cos(2 ̄'fot+ ) (1)
傳統(tǒng)實時數(shù)字傅里葉變換(DFT)表達(dá)式:
L
( )=∑ (n)e (3)
n=1
其中,J表示累加的采樣點數(shù)。DFT結(jié)果的頻率分辨率:
其中 表示采樣頻率, 表示 個采樣點 fo= / =1/ ,
其中:s(t)表示GNSS衛(wèi)星擴頻信號,n(t)表示高斯白噪
聲, 表示連續(xù)波干擾的數(shù)目, , 分別表示不同連續(xù)
波干擾的幅度、頻率、相位。以GPS為例,地表的衛(wèi)星信號功
率在一130dBm左右,帶寬為2MHz。擴頻信號本身具有一定
的抗窄帶干擾能力,但當(dāng)干擾功率超過一110dBm左右時,開
始對接收機產(chǎn)生影響。隨著干擾功率的進一步增加,接收機
將失去鎖定、無法定位。這些窄帶干擾處于有效信號帶寬
內(nèi),不能通過射頻中的帶通濾波器去除,所以需要額外的干
擾抑制處理。
圖1所示為具有窄帶干擾抑制的多模衛(wèi)星導(dǎo)航接收機
結(jié)構(gòu)框圖。射頻對不同GNSS系統(tǒng)的衛(wèi)星信號進行下變頻、
量化等操作,得到數(shù)字中頻信號。干擾檢測模塊首先檢測中
頻信號中的窄帶干擾頻率和功率,然后將這些信息轉(zhuǎn)化為濾
波器參數(shù)置人干擾去除模塊。干擾去除模塊對輸入信號進
行濾波后輸出給接收機基帶進行捕獲、跟蹤及定位解算,最
后輸出定位等信息。
圖1 GNSS接收機中的窄帶干擾抑制
圖1中干擾去除模塊使用級聯(lián)的IIR陷波器陣列;干擾
檢測模塊使用數(shù)字傅里葉變換(DFT)對輸入信號進行頻域
變換,具有實時性且硬件開銷小,適合集成電路實現(xiàn)。
2.1 IIR陷波器陣列
陷波器能夠有效抑制窄帶干擾,去除干擾頻率而保留大
部分衛(wèi)星信號功率。其中無限沖擊響應(yīng)(IIR)類型的陷波器
由于其實現(xiàn)和計算復(fù)雜度低,而被廣泛用于擴頻通信中的窄
帶干擾抑制。對于實數(shù)輸入信號,其傅里葉變換會在正負(fù)頻
率處產(chǎn)生兩個干擾峰,需要二階IIR陷波器,其傳遞函數(shù) :
(2)
其中:gO。=2=fo./f.表示IIR陷波器歸一化中心頻率;0
的累加時間。則以分辨率 檢測帶寬為BW的頻率范圍所
需的檢測時間為:T=BW/A,傳統(tǒng)DFT頻率檢測算法及其檢
測時間與頻率分辨率成平方反比關(guān)系。以GPS衛(wèi)星信號為
例,需要檢測的頻率帶寬為2MHz,要達(dá)到100Hz的頻率分辨
率,則需要200秒的檢測時間,檢測時間太長而無法接受。
多級DFT窄帶干擾頻率檢測算法在保證最終檢測結(jié)果
分辨率不變的情況下,能大幅度減小檢測時間。如圖2所示
其原理為:DFT檢測分為多個等級,每個等級僅有小部分窄
帶干擾所在的頻率分量處會有大的幅度超過門限,從而被記
錄下來由下一級更小的頻率分辨率繼續(xù)檢測,大部分不存在
干擾的頻率范圍被排除掉,不需要下一級的檢測。依次類
推,從而減少總體檢測時間。
圖2 多級DFT窄帶干擾頻率檢測算法
圖2中的多級DFT窄帶干擾頻率檢測算法的級數(shù) 和每
級的頻率分辨率組合[ …. ]選擇不同,最終檢測時
間差別很大。存在最優(yōu)的頻率組合,使檢測時問最小。假設(shè)
使用n級DFT檢測BW的頻率范圍,檢測時間可建模如下:
妻.~ 一
: +j (爭+._…+等) (4)
1 j 2 J n
其中 表示窄帶干擾數(shù)目,求多變量的極值,得方程
鯧:
:
一
二 業(yè)+.『N:0
。一
0T N
-
2 N
fl
—
:—:
一 。 一 0
—
+了
(5)
: + :0
一
。一
,
解遞歸方程組,得到檢測時間最小化條件,即最優(yōu)頻率
組合:
f - n-
{ … (6)
2r,BW,f ̄-1
圖3中,假設(shè)最終分辨率為lOOHz時,傳統(tǒng)DFr需要檢
測時間200s,而使用多級Db3'檢測算法,令等級L=2、3、4
時,分別需要檢測時間為2.052s、0.636s、0.420s,分別是傳統(tǒng)
檢測時間的1.026%、0.318%、0.21%。且隨著檢測等級的
不斷增加,檢測時間不斷減小。
3.2多個窄帶干擾抑制仿真
I
l 一
輸入GPS信號功率一131dBm,中頻為4MHz。同時存在
兩個窄帶干擾信號,第一個干擾頻率為4MHz,功率為一
85dBm,第二個干擾頻率為3MHz,功率為-91 dBm。經(jīng)過2個
IIR陷波器抑制后輸出。分別對輸入輸出信號進行捕獲,結(jié)
l =嘉
由式(6)得到最優(yōu)的頻率組合后,帶入式(4)得到檢測
時間的估計。
2.3掃頻信號跟蹤
窄帶干擾頻率有可能隨時間變化形成掃頻干擾,因此需
要對窄帶干擾的頻率進行跟蹤抑制。本文使用DFT檢測干
擾頻率后更新濾波器參數(shù)的開環(huán)方法對干擾進行跟蹤抑
制。假設(shè)已初始化完成陷波器對一個窄帶干擾進行去除,該
陷波器工作中心頻率為∞。,有效帶寬為BW。對輸入信號在
中心頻率為∞。,帶寬為BW的范圍內(nèi)進行DFT檢測后得到新
的干擾頻率cc, ,則調(diào)整陷波器的參數(shù),使其陷波中心頻率變
為∞.。如此循環(huán),將保持陷波器對該干擾頻率的持續(xù)抑制。
:2 m O
以上的開環(huán)跟蹤算法有一個前提假設(shè)——在DFT檢測
時間內(nèi),干擾頻率是不變的,或者是變化緩慢的。由此可對
該算法所能跟蹤的頻率變化率進行建模,假設(shè)DFT頻率分辨
率為,0,其跟蹤頻率變化的能力:
Js 。印 = T = = -vz (Hz/s)) (7)
3仿真
3.1 最優(yōu)化的多級DFT窄帶干擾檢測速度仿真
仍以GPS信號為例,帶寬BW=2MHz,假設(shè)存在10個窄
帶干擾,不同等級和最終分辨率的多級DFT頻率檢測時間如
圖3。
圖3 多級DFT頻率檢測算法檢測時間
.---——
74..-.——
果如圖4、5、6所示。
圖4是未進行窄帶干擾抑制,直接對輸入信號進行捕獲
的結(jié)果,可以看到頻率和碼相位平面上沒有明顯相關(guān)峰值,
無法捕獲信號。
頻率域0 0 碼相位域
圖4輸入信號捕獲結(jié)果(2個窄帶干擾)
圖5是陷波器1對4MHz附近的窄帶干擾抑制后,對其
輸出信號進行捕獲。由于GPS信號處于4MHz,該頻率處的
窄帶干擾對信號影響最大,經(jīng)過陷波器1的抑制后,GPS信
號質(zhì)量得到提高,捕獲結(jié)果中已經(jīng)可以看到相關(guān)峰,但是峰
值與噪底區(qū)分不大——信噪比低,原因是還有一個窄帶干擾
沒有抑制。此時的信噪比相比無干擾情況時損失了4.5dB。
圖6是陷波器2對3MHz附近的窄帶干擾抑制后,對其
輸出信號進行捕獲。此時兩個窄帶干擾都已經(jīng)被抑制了,所
以GPS信號的信噪比進一步提升。最后的信噪比相比無干
擾情況時只損失了0.5dB。
3.3跟蹤掃頻信號仿真
由式(7)可知,檢測加調(diào)整濾波器參數(shù)的開環(huán)跟蹤算法
的跟蹤能力僅與DFT檢測頻率分辨率的平方成正比關(guān)系。
信號初始頻率4kHz,DFT頻率分辨率10Hz情況下,跟蹤能
力仿真如圖7所示。
由式(7)可知,在頻率分辨率為10Hz時,其跟蹤能力極
限為50Hz/s,圖4的仿真結(jié)果表明對于40Hz/s及以下的頻
率斜升干擾信號,跟蹤良好;當(dāng)頻率變化率達(dá)到50Hz/s時,
已無法正確跟蹤,與式(7)理論結(jié)果符合。
l2000
tO000
8000
倒6000
驟
4oo0
2O00
O
30
頻率域0 0 碼相位域
圖5陷波器1輸出信號捕獲結(jié)果
頻率域0 0 碼相位域
圖6陷波器2輸出信號捕獲結(jié)果
不問斜率的掃描信號蝶蹤站果
o 2 4 6 8 1o
時間/s
圖7跟蹤頻率隨時間變化信號
4
結(jié)束語
本文針對衛(wèi)星導(dǎo)航接收機中窄帶干擾抑制問題,使用頻
域檢測和時域濾波組合的結(jié)構(gòu)。提出了最優(yōu)化的多級DFT
O 0枷 ㈣ O 0 湖 鯽 O 0 枷 猢 OO , *聰
窄帶干擾頻率檢測算法,該算法利用傳統(tǒng)DFT易于實現(xiàn)的優(yōu)
點,同時克服了傳統(tǒng)方法檢測時間長的缺點,仿真結(jié)果表明
該算法的檢測時間是傳統(tǒng)DFT檢測時間的1%甚至更少。
對多個-85dBm左右的窄帶干擾的抑制情況進行了仿真,結(jié)
果表明濾波后信噪比損失小于ldB。利用DFT檢測干擾頻
率并調(diào)整陷波器系數(shù)的開環(huán)方式對輸入的掃頻信號進行跟
蹤,該結(jié)構(gòu)跟蹤性能僅與DFT檢測頻率分辨率有關(guān),通過調(diào)
整DFT檢測頻率分辨率可適應(yīng)不同應(yīng)用場景。
參考文獻:
[1]C J Hegarty.GNSS signals—An overview[C].Frequency Control
Symposium(FCS).2012:1-7.
[2]A Rusu-Casandra,I Marghescu,E S Lohan.Impact of nariDw.
band interference on unambiguousacquisition approaches in Galileo
[C].Localization and GNSS(ICL—GNSS),2011:127—132.
[3]Liu Yong-Jun,Hu Han—Ying.The Research on GPS Frequency
Domain Anti—Jamming Algorithms[C].Wireless Communications,
Networking andMobile Computing.2009:1-3.
[4]Zhang Jie,E Lohan.Effect and mitigation of na ̄owband interfer-
enee on Galileo E1 signal acquisition and tracking accuracy[C].
Localization and GNSS(ICL—GNSS),201 1:36—41.
[5] Chien Ying-Ren,et a1.A Novel Continuous Wave Interference
Detectable Adaptive Notch Filter for GPS Receivers[C].Global
Telecommunications Conference.2010:1—6.
[6]M Nouri,S A Aghdam,V T Vakili.An optimal method for nar.
rowband interferencemitigation in the GPS[C].Modeling,Simula—
tion and Applied Optimization(ICMSAO),201 1 4th International
Conference,2011:1—4.
[7]C Chang hung-Liang.Novel multiplexing technique in anti-jam—
ming GNSS receiver[C].American Control Conference(ACC),
2011:3453—3458.
[8] Yang Xiao—bo.An Improved GPS Receiver Anti—jammer Algo—
rithm Based on Space-time Adaptive Processing[C].Computer
Distirbuted Control and Intelligent EnvironmentalMonitoring(CD—
CIEM).2012:106—109.
[作者簡介]
藺曉龍(1985一),男(漢族),河南平頂山人,博士研
究生,研究方向為衛(wèi)星導(dǎo)航接收機信號處理算法;
何文濤(1982一),男(漢族),湖北黃岡人,博士,助
理研究員,研究方向為衛(wèi)星導(dǎo)航接收機信號處理算
法及超大規(guī)模集成電路設(shè)計;
徐建華(1963一),男(漢族),浙江金華人,研究員,博士生導(dǎo)師,研究
方向為GPS接收機、FPGA平臺開發(fā)、芯片系統(tǒng)設(shè)計等方面;
葉甜春(1965一),男(漢族),四川綿陽人,研究員,博士生導(dǎo)師,研究
方向為集成電路設(shè)計與應(yīng)用、納米加工制造技術(shù)、新型器件等方面。
一
75一

本文發(fā)布于:2023-11-03 01:50:59,感謝您對本站的認(rèn)可!
本文鏈接:http://www.newhan.cn/zhishi/a/169894745927207.html
版權(quán)聲明:本站內(nèi)容均來自互聯(lián)網(wǎng),僅供演示用,請勿用于商業(yè)和其他非法用途。如果侵犯了您的權(quán)益請與我們聯(lián)系,我們將在24小時內(nèi)刪除。
本文word下載地址:衛(wèi)星導(dǎo)航接收機中窄帶干擾抑制算法.doc
本文 PDF 下載地址:衛(wèi)星導(dǎo)航接收機中窄帶干擾抑制算法.pdf
| 留言與評論(共有 0 條評論) |