
第31卷第1期 計(jì)算機(jī)仿真 2014年1月
文章編號(hào):1006—9348(2014)O1—0072-04
衛(wèi)星導(dǎo)航接收機(jī)中窄帶干擾抑制算法
藺曉龍,何文濤,徐建華,葉甜春
(中國(guó)科學(xué)院微電子研究所,北京100029)
摘要:窄帶干擾能夠嚴(yán)重降低衛(wèi)星導(dǎo)航接收機(jī)的性能。為提高衛(wèi)星導(dǎo)航性能,提出一種應(yīng)用于衛(wèi)星導(dǎo)航接收機(jī)中的窄帶干
擾抑制算法。首先利用數(shù)字傅里葉變換(DFT)檢測(cè)窄帶干擾,傳統(tǒng)DFF實(shí)現(xiàn)簡(jiǎn)單、便于集成,但計(jì)算慢、檢測(cè)時(shí)間過(guò)長(zhǎng)。為
解決上述問(wèn)題,采用最優(yōu)化的多級(jí)DFT頻率檢測(cè)算法,能顯著減小檢測(cè)時(shí)間。然后使用級(jí)聯(lián)的IIR陷波器抑制窄帶干擾。
最后使用頻率檢測(cè)和濾波器更新相結(jié)合的開(kāi)環(huán)結(jié)構(gòu),用于跟蹤掃頻信號(hào)并進(jìn)行抑制。仿真結(jié)果表明,改進(jìn)算法能夠有效的
抑制衛(wèi)星導(dǎo)航接收機(jī)中的窄帶和掃頻干擾。
關(guān)鍵詞:衛(wèi)星導(dǎo)航接收機(jī);窄帶干擾;數(shù)字傅里葉變換;掃頻干擾
中圖分類號(hào):TN911 文獻(xiàn)標(biāo)識(shí)碼: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以及中國(guó)的北斗衛(wèi)星導(dǎo)航系統(tǒng)
持續(xù)抑制。
目前有多種窄帶干擾抑制方法,包括:文獻(xiàn)[3][4]使用
F 技術(shù)在頻域?qū)φ瓗Ц蓴_進(jìn)行檢測(cè)和抑制。由于使用
FFT,需要大量緩存數(shù)據(jù),計(jì)算量大。文獻(xiàn)[5][6]使用自適
應(yīng)陷波器(ANF)技術(shù),LMS、RLS等多種自適應(yīng)算法被采用,
可使用IIR型或FIR型的濾波器。其中自適應(yīng)算法計(jì)算量、
檢測(cè)數(shù)目、收斂性需要特別關(guān)注。文獻(xiàn)[7]使用檢測(cè)加抑制
的結(jié)構(gòu)對(duì)包括窄帶干擾在內(nèi)的多種干擾進(jìn)行抑制。文獻(xiàn)[8]
使用其它更為復(fù)雜的一些變換與濾波技術(shù),能夠抑制包括多
(BDS)以及其它的增強(qiáng)系統(tǒng)統(tǒng)稱為全球?qū)Ш叫l(wèi)星系統(tǒng)
(GNSS)。地面接收到的GNSS信號(hào)強(qiáng)度在一130dBm左右,
遠(yuǎn)在噪聲功率之下 。隨著定位功能越來(lái)越多的集成到移
動(dòng)設(shè)備中,由其它通信設(shè)備泄露到GNSS信號(hào)頻段的干擾不
斷增多。這些干擾會(huì)嚴(yán)重降低衛(wèi)星導(dǎo)航接收機(jī)的性能。其
中對(duì)GNSS接收機(jī)影響最大,且最常見(jiàn)的干擾類型是窄帶干
擾 。
另外接收機(jī)與干擾源之間的相對(duì)運(yùn)動(dòng)會(huì)產(chǎn)生多普勒效
種類型的干擾,考慮了一定的掃頻干擾,但算法及實(shí)現(xiàn)復(fù)雜
度都較高,不易于集成。
本文使用檢測(cè)加濾波的結(jié)構(gòu)實(shí)現(xiàn)衛(wèi)星導(dǎo)航接收機(jī)中的
應(yīng),此時(shí)即使外界窄帶干擾頻率是固定不變的,GNSS接收機(jī)
內(nèi)部,其干擾頻率會(huì)隨時(shí)間發(fā)生變化,形成掃頻干擾信號(hào)。
所以在處理窄帶干擾抑制時(shí)需要考慮跟蹤干擾頻率以達(dá)到
基金項(xiàng)目:中國(guó)第二代衛(wèi)星導(dǎo)航系統(tǒng)重大專項(xiàng)多模導(dǎo)航型基帶芯片
(2011GFZX03030204)
窄帶干擾抑制。DFT技術(shù)由于其硬件實(shí)現(xiàn)簡(jiǎn)單,便于集成而
被廣泛采用,但是由于衛(wèi)星接收機(jī)中信號(hào)檢測(cè)帶寬很大,如
2MHz,頻率分辨率又要求很小,如100Hz,傳統(tǒng)DFT計(jì)算緩
慢,檢測(cè)時(shí)間過(guò)長(zhǎng)而無(wú)實(shí)際應(yīng)用價(jià)值。針對(duì)這個(gè)問(wèn)題,本文
提出的多級(jí)DFT窄帶干擾頻率檢測(cè)算法,使用多等級(jí)的檢
收稿日期:2013—04—10
72一
測(cè),從而最小化檢測(cè)時(shí)間。然后使用級(jí)聯(lián)的IIR陷波器陣
列,去除多個(gè)輸入通道中的多個(gè)窄帶干擾。最后采用DFT檢 波器能去除一個(gè)實(shí)數(shù)輸入的窄帶干擾,要應(yīng)對(duì)多個(gè)窄帶干
測(cè)加調(diào)整陷波器的開(kāi)環(huán)方式對(duì)掃頻干擾進(jìn)行抑制。
<b<1表示與IIR陷波器帶寬有關(guān)的參數(shù)。單個(gè)二階IIR陷
擾,需要多級(jí)級(jí)聯(lián)的二階IIR陷波器序列。
2.2 窄帶干擾頻率檢測(cè)
2衛(wèi)星導(dǎo)航接收機(jī)干擾抑制結(jié)構(gòu)
受多個(gè)連續(xù)波干擾的擴(kuò)頻信號(hào)模型:
(£)=5(f)+n(t)+ A cos(2 ̄'fot+ ) (1)
傳統(tǒng)實(shí)時(shí)數(shù)字傅里葉變換(DFT)表達(dá)式:
L
( )=∑ (n)e (3)
n=1
其中,J表示累加的采樣點(diǎn)數(shù)。DFT結(jié)果的頻率分辨率:
其中 表示采樣頻率, 表示 個(gè)采樣點(diǎn) fo= / =1/ ,
其中:s(t)表示GNSS衛(wèi)星擴(kuò)頻信號(hào),n(t)表示高斯白噪
聲, 表示連續(xù)波干擾的數(shù)目, , 分別表示不同連續(xù)
波干擾的幅度、頻率、相位。以GPS為例,地表的衛(wèi)星信號(hào)功
率在一130dBm左右,帶寬為2MHz。擴(kuò)頻信號(hào)本身具有一定
的抗窄帶干擾能力,但當(dāng)干擾功率超過(guò)一110dBm左右時(shí),開(kāi)
始對(duì)接收機(jī)產(chǎn)生影響。隨著干擾功率的進(jìn)一步增加,接收機(jī)
將失去鎖定、無(wú)法定位。這些窄帶干擾處于有效信號(hào)帶寬
內(nèi),不能通過(guò)射頻中的帶通濾波器去除,所以需要額外的干
擾抑制處理。
圖1所示為具有窄帶干擾抑制的多模衛(wèi)星導(dǎo)航接收機(jī)
結(jié)構(gòu)框圖。射頻對(duì)不同GNSS系統(tǒng)的衛(wèi)星信號(hào)進(jìn)行下變頻、
量化等操作,得到數(shù)字中頻信號(hào)。干擾檢測(cè)模塊首先檢測(cè)中
頻信號(hào)中的窄帶干擾頻率和功率,然后將這些信息轉(zhuǎn)化為濾
波器參數(shù)置人干擾去除模塊。干擾去除模塊對(duì)輸入信號(hào)進(jìn)
行濾波后輸出給接收機(jī)基帶進(jìn)行捕獲、跟蹤及定位解算,最
后輸出定位等信息。
圖1 GNSS接收機(jī)中的窄帶干擾抑制
圖1中干擾去除模塊使用級(jí)聯(lián)的IIR陷波器陣列;干擾
檢測(cè)模塊使用數(shù)字傅里葉變換(DFT)對(duì)輸入信號(hào)進(jìn)行頻域
變換,具有實(shí)時(shí)性且硬件開(kāi)銷(xiāo)小,適合集成電路實(shí)現(xiàn)。
2.1 IIR陷波器陣列
陷波器能夠有效抑制窄帶干擾,去除干擾頻率而保留大
部分衛(wèi)星信號(hào)功率。其中無(wú)限沖擊響應(yīng)(IIR)類型的陷波器
由于其實(shí)現(xiàn)和計(jì)算復(fù)雜度低,而被廣泛用于擴(kuò)頻通信中的窄
帶干擾抑制。對(duì)于實(shí)數(shù)輸入信號(hào),其傅里葉變換會(huì)在正負(fù)頻
率處產(chǎn)生兩個(gè)干擾峰,需要二階IIR陷波器,其傳遞函數(shù) :
(2)
其中:gO。=2=fo./f.表示IIR陷波器歸一化中心頻率;0
的累加時(shí)間。則以分辨率 檢測(cè)帶寬為BW的頻率范圍所
需的檢測(cè)時(shí)間為:T=BW/A,傳統(tǒng)DFT頻率檢測(cè)算法及其檢
測(cè)時(shí)間與頻率分辨率成平方反比關(guān)系。以GPS衛(wèi)星信號(hào)為
例,需要檢測(cè)的頻率帶寬為2MHz,要達(dá)到100Hz的頻率分辨
率,則需要200秒的檢測(cè)時(shí)間,檢測(cè)時(shí)間太長(zhǎng)而無(wú)法接受。
多級(jí)DFT窄帶干擾頻率檢測(cè)算法在保證最終檢測(cè)結(jié)果
分辨率不變的情況下,能大幅度減小檢測(cè)時(shí)間。如圖2所示
其原理為:DFT檢測(cè)分為多個(gè)等級(jí),每個(gè)等級(jí)僅有小部分窄
帶干擾所在的頻率分量處會(huì)有大的幅度超過(guò)門(mén)限,從而被記
錄下來(lái)由下一級(jí)更小的頻率分辨率繼續(xù)檢測(cè),大部分不存在
干擾的頻率范圍被排除掉,不需要下一級(jí)的檢測(cè)。依次類
推,從而減少總體檢測(cè)時(shí)間。
圖2 多級(jí)DFT窄帶干擾頻率檢測(cè)算法
圖2中的多級(jí)DFT窄帶干擾頻率檢測(cè)算法的級(jí)數(shù) 和每
級(jí)的頻率分辨率組合[ …. ]選擇不同,最終檢測(cè)時(shí)
間差別很大。存在最優(yōu)的頻率組合,使檢測(cè)時(shí)問(wèn)最小。假設(shè)
使用n級(jí)DFT檢測(cè)BW的頻率范圍,檢測(cè)時(shí)間可建模如下:
妻.~ 一
: +j (爭(zhēng)+._…+等) (4)
1 j 2 J n
其中 表示窄帶干擾數(shù)目,求多變量的極值,得方程
鯧:
:
一
二 業(yè)+.『N:0
。一
0T N
-
2 N
fl
—
:—:
一 。 一 0
—
+了
(5)
: + :0
一
。一
,
解遞歸方程組,得到檢測(cè)時(shí)間最小化條件,即最優(yōu)頻率
組合:
f - n-
{ … (6)
2r,BW,f ̄-1
圖3中,假設(shè)最終分辨率為lOOHz時(shí),傳統(tǒng)DFr需要檢
測(cè)時(shí)間200s,而使用多級(jí)Db3'檢測(cè)算法,令等級(jí)L=2、3、4
時(shí),分別需要檢測(cè)時(shí)間為2.052s、0.636s、0.420s,分別是傳統(tǒng)
檢測(cè)時(shí)間的1.026%、0.318%、0.21%。且隨著檢測(cè)等級(jí)的
不斷增加,檢測(cè)時(shí)間不斷減小。
3.2多個(gè)窄帶干擾抑制仿真
I
l 一
輸入GPS信號(hào)功率一131dBm,中頻為4MHz。同時(shí)存在
兩個(gè)窄帶干擾信號(hào),第一個(gè)干擾頻率為4MHz,功率為一
85dBm,第二個(gè)干擾頻率為3MHz,功率為-91 dBm。經(jīng)過(guò)2個(gè)
IIR陷波器抑制后輸出。分別對(duì)輸入輸出信號(hào)進(jìn)行捕獲,結(jié)
l =嘉
由式(6)得到最優(yōu)的頻率組合后,帶入式(4)得到檢測(cè)
時(shí)間的估計(jì)。
2.3掃頻信號(hào)跟蹤
窄帶干擾頻率有可能隨時(shí)間變化形成掃頻干擾,因此需
要對(duì)窄帶干擾的頻率進(jìn)行跟蹤抑制。本文使用DFT檢測(cè)干
擾頻率后更新濾波器參數(shù)的開(kāi)環(huán)方法對(duì)干擾進(jìn)行跟蹤抑
制。假設(shè)已初始化完成陷波器對(duì)一個(gè)窄帶干擾進(jìn)行去除,該
陷波器工作中心頻率為∞。,有效帶寬為BW。對(duì)輸入信號(hào)在
中心頻率為∞。,帶寬為BW的范圍內(nèi)進(jìn)行DFT檢測(cè)后得到新
的干擾頻率cc, ,則調(diào)整陷波器的參數(shù),使其陷波中心頻率變
為∞.。如此循環(huán),將保持陷波器對(duì)該干擾頻率的持續(xù)抑制。
:2 m O
以上的開(kāi)環(huán)跟蹤算法有一個(gè)前提假設(shè)——在DFT檢測(cè)
時(shí)間內(nèi),干擾頻率是不變的,或者是變化緩慢的。由此可對(duì)
該算法所能跟蹤的頻率變化率進(jìn)行建模,假設(shè)DFT頻率分辨
率為,0,其跟蹤頻率變化的能力:
Js 。印 = T = = -vz (Hz/s)) (7)
3仿真
3.1 最優(yōu)化的多級(jí)DFT窄帶干擾檢測(cè)速度仿真
仍以GPS信號(hào)為例,帶寬BW=2MHz,假設(shè)存在10個(gè)窄
帶干擾,不同等級(jí)和最終分辨率的多級(jí)DFT頻率檢測(cè)時(shí)間如
圖3。
圖3 多級(jí)DFT頻率檢測(cè)算法檢測(cè)時(shí)間
.---——
74..-.——
果如圖4、5、6所示。
圖4是未進(jìn)行窄帶干擾抑制,直接對(duì)輸入信號(hào)進(jìn)行捕獲
的結(jié)果,可以看到頻率和碼相位平面上沒(méi)有明顯相關(guān)峰值,
無(wú)法捕獲信號(hào)。
頻率域0 0 碼相位域
圖4輸入信號(hào)捕獲結(jié)果(2個(gè)窄帶干擾)
圖5是陷波器1對(duì)4MHz附近的窄帶干擾抑制后,對(duì)其
輸出信號(hào)進(jìn)行捕獲。由于GPS信號(hào)處于4MHz,該頻率處的
窄帶干擾對(duì)信號(hào)影響最大,經(jīng)過(guò)陷波器1的抑制后,GPS信
號(hào)質(zhì)量得到提高,捕獲結(jié)果中已經(jīng)可以看到相關(guān)峰,但是峰
值與噪底區(qū)分不大——信噪比低,原因是還有一個(gè)窄帶干擾
沒(méi)有抑制。此時(shí)的信噪比相比無(wú)干擾情況時(shí)損失了4.5dB。
圖6是陷波器2對(duì)3MHz附近的窄帶干擾抑制后,對(duì)其
輸出信號(hào)進(jìn)行捕獲。此時(shí)兩個(gè)窄帶干擾都已經(jīng)被抑制了,所
以GPS信號(hào)的信噪比進(jìn)一步提升。最后的信噪比相比無(wú)干
擾情況時(shí)只損失了0.5dB。
3.3跟蹤掃頻信號(hào)仿真
由式(7)可知,檢測(cè)加調(diào)整濾波器參數(shù)的開(kāi)環(huán)跟蹤算法
的跟蹤能力僅與DFT檢測(cè)頻率分辨率的平方成正比關(guān)系。
信號(hào)初始頻率4kHz,DFT頻率分辨率10Hz情況下,跟蹤能
力仿真如圖7所示。
由式(7)可知,在頻率分辨率為10Hz時(shí),其跟蹤能力極
限為50Hz/s,圖4的仿真結(jié)果表明對(duì)于40Hz/s及以下的頻
率斜升干擾信號(hào),跟蹤良好;當(dāng)頻率變化率達(dá)到50Hz/s時(shí),
已無(wú)法正確跟蹤,與式(7)理論結(jié)果符合。
l2000
tO000
8000
倒6000
驟
4oo0
2O00
O
30
頻率域0 0 碼相位域
圖5陷波器1輸出信號(hào)捕獲結(jié)果
頻率域0 0 碼相位域
圖6陷波器2輸出信號(hào)捕獲結(jié)果
不問(wèn)斜率的掃描信號(hào)蝶蹤站果
o 2 4 6 8 1o
時(shí)間/s
圖7跟蹤頻率隨時(shí)間變化信號(hào)
4
結(jié)束語(yǔ)
本文針對(duì)衛(wèi)星導(dǎo)航接收機(jī)中窄帶干擾抑制問(wèn)題,使用頻
域檢測(cè)和時(shí)域?yàn)V波組合的結(jié)構(gòu)。提出了最優(yōu)化的多級(jí)DFT
O 0枷 ㈣ O 0 湖 鯽 O 0 枷 猢 OO , *聰
窄帶干擾頻率檢測(cè)算法,該算法利用傳統(tǒng)DFT易于實(shí)現(xiàn)的優(yōu)
點(diǎn),同時(shí)克服了傳統(tǒng)方法檢測(cè)時(shí)間長(zhǎng)的缺點(diǎn),仿真結(jié)果表明
該算法的檢測(cè)時(shí)間是傳統(tǒng)DFT檢測(cè)時(shí)間的1%甚至更少。
對(duì)多個(gè)-85dBm左右的窄帶干擾的抑制情況進(jìn)行了仿真,結(jié)
果表明濾波后信噪比損失小于ldB。利用DFT檢測(cè)干擾頻
率并調(diào)整陷波器系數(shù)的開(kāi)環(huán)方式對(duì)輸入的掃頻信號(hào)進(jìn)行跟
蹤,該結(jié)構(gòu)跟蹤性能僅與DFT檢測(cè)頻率分辨率有關(guān),通過(guò)調(diào)
整DFT檢測(cè)頻率分辨率可適應(yīng)不同應(yīng)用場(chǎng)景。
參考文獻(xiàn):
[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.
[作者簡(jiǎn)介]
藺曉龍(1985一),男(漢族),河南平頂山人,博士研
究生,研究方向?yàn)樾l(wèi)星導(dǎo)航接收機(jī)信號(hào)處理算法;
何文濤(1982一),男(漢族),湖北黃岡人,博士,助
理研究員,研究方向?yàn)樾l(wèi)星導(dǎo)航接收機(jī)信號(hào)處理算
法及超大規(guī)模集成電路設(shè)計(jì);
徐建華(1963一),男(漢族),浙江金華人,研究員,博士生導(dǎo)師,研究
方向?yàn)椋牵校咏邮諜C(jī)、FPGA平臺(tái)開(kāi)發(fā)、芯片系統(tǒng)設(shè)計(jì)等方面;
葉甜春(1965一),男(漢族),四川綿陽(yáng)人,研究員,博士生導(dǎo)師,研究
方向?yàn)榧呻娐吩O(shè)計(jì)與應(yīng)用、納米加工制造技術(shù)、新型器件等方面。
一
75一

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