基于FPGA的窄帶干擾抑制算法在GPS中的應(yīng)用

GPS傳輸過程中信號電平非常低，地面用戶接收到的GPS信號功率約為-160 dBW，很容易受到各種有意或無意信號的干擾[1]。按干擾的頻譜寬度相對于擴(kuò)頻信號帶寬的大小，可分為寬帶干擾和窄帶干擾。窄帶干擾產(chǎn)生容易，功率譜密度高，幾個(gè)窄帶干擾就可覆蓋一定的寬度，所以相對于寬帶壓制性干擾而言，窄帶干擾更有利于干擾GPS通信。

頻域抑制的方法能夠很好地抑制窄帶干擾。由于GPS信號的擴(kuò)頻特性，信號功率要遠(yuǎn)遠(yuǎn)低于噪聲功率，期望信號的頻譜在頻域上為平坦的白噪聲特性。當(dāng)存在窄帶干擾信號的時(shí)候，在干擾的頻域就會出現(xiàn)比較高的頻域特性，這樣就可以在頻域采用相應(yīng)的方法進(jìn)行干擾抑制。文中提出了一種基于FPGA的重疊FFT頻域窄帶干擾抑制算法的實(shí)現(xiàn)方案。測試結(jié)果表明，該方案對于窄帶干擾具有較好的抑制作用。

1 重疊FFT頻域抑制窄帶干擾的原理

重疊FFT頻域干擾抑制算法的原理是窄帶干擾相對于擴(kuò)頻信號能量集中在很窄的頻帶內(nèi)，在頻域上表現(xiàn)為很窄的尖峰，因此可先將混合信號變換到頻域，檢測出干擾的頻譜位置，將這些譜線置零或進(jìn)行衰減，最后反變換到時(shí)域進(jìn)行解擴(kuò)[2]。重疊FFT頻域抑制窄帶干擾的原理如圖1所示。

圖1 頻域抑制窄帶干擾流程圖Fig.1 Flow chart of NBI suppression in the frequency domain

1.1 加窗處理

實(shí)際信號是無限長的，因此數(shù)據(jù)截?cái)嗖豢杀苊?。由于?shù)據(jù)截?cái)嗟挠绊懀肍FT方法估計(jì)的信號存在頻譜泄露。在不加窗（相當(dāng)于加矩形窗）的情況下，離散傅里葉變換的第一旁瓣只有-13 dB，存在嚴(yán)重的頻譜泄露。通常采用對分段數(shù)據(jù)加窗的方法減小頻譜泄露。從時(shí)域來看，加窗本質(zhì)上是對輸入數(shù)據(jù)進(jìn)行加權(quán)處理，窗函數(shù)從中心向兩端逐步衰減，保證了數(shù)據(jù)段兩端的平滑，從而減小頻譜泄露。然而加窗卻會使輸入信號產(chǎn)生畸變，從而帶來額外的信噪比損耗[2]。重疊加窗可以減小信噪比損耗，重疊的比例越大，加窗損耗就越小。但是重疊比例增大意味著計(jì)算量增大。實(shí)際中重疊比例的選擇取決于硬件條件和性能要求。本方案采用1/2重疊。

1.2 數(shù)據(jù)輸出

時(shí)域上窗函數(shù)從中間向兩端逐漸衰減，越靠近兩端，衰減越嚴(yán)重，引起的信號畸變越大，從而帶來的信噪比損耗越大。在兩路加窗數(shù)據(jù)合成時(shí)采用相加輸出的方法，即將兩路數(shù)據(jù)的重疊部分相加作為最終輸出，這樣相對于傳統(tǒng)的選擇輸出方法可以進(jìn)一步減小信噪比損失[3]。

1.3 門限檢測

在重疊加窗的結(jié)構(gòu)下，可以采用多種頻域干擾消除技術(shù)。目前常用的窄帶干擾抑制算法有K譜線法、N-sigma算法、最大似然法、自適應(yīng)多門限檢測算法。前3種算法在沒有窄帶干擾存在時(shí)能夠?qū)π盘柟β收_估計(jì)，可以保證不會抑制有用信號。但是當(dāng)存在窄帶干擾時(shí)，接收信號不滿足高斯分布，需要適當(dāng)調(diào)整門限。如果門限調(diào)整不合適，就會對有用信號抑制過度，或者干擾抑制不足，無法達(dá)到預(yù)期的效果。多門限檢測算法[4]可以克服上述不足，該算法步驟如下：

1）對N根譜線的模平方求取均值，作1/λ為的估計(jì)值。

2）求8/λ。

3）對N根譜線進(jìn)行統(tǒng)計(jì)，若有模平方超過8/λ的譜線則認(rèn)為其中存在著窄帶干擾可以將其置零。

4）返回1）注意是對所有未歸零的譜線求取平均值，已歸零譜線不計(jì)算在內(nèi)。重復(fù)以上過程直到?jīng)]有超過門限的譜線。

2 窄帶干擾抑制算法的FPGA實(shí)現(xiàn)

窄帶抑制算法在Xilinx公司的Virtex5芯片上實(shí)現(xiàn)。Virtex5具有資源豐富，運(yùn)算速度快的優(yōu)點(diǎn)。所用的描述語言為verilog HDL，整個(gè)窄帶抑制部分包括win1024，fft_fifio，F(xiàn)FT1024，jam_reover，IFFT1024，fft_fifo2等模塊，分為上下兩個(gè)支路。輸入信號為經(jīng)過A/D采樣和數(shù)字下變頻處理的GPS信號。程序設(shè)計(jì)如圖2所示。硬件實(shí)現(xiàn)平臺如圖3所示。

圖2 程序設(shè)計(jì)框圖Fig.2 Block diagram for programming

圖3 硬件實(shí)物圖Fig.3 Actual photo of hardware

2.1 重疊加窗的實(shí)現(xiàn)

win1024是存儲深度為1 024的rom，用于存儲加窗系數(shù)，窗函數(shù)為Blackman窗。Blackman窗函數(shù)的旁瓣抑制效果較好[5]，其旁瓣抑制可以達(dá)到-57 dB。N/2延遲通過啟動(dòng)信號start1、start2來完成，start1、start2均為高有效。start1、start2分別通過兩個(gè)計(jì)數(shù)器st1_cnt、cnt控制。計(jì)數(shù)器時(shí)鐘均為10M。計(jì)數(shù)器st1_cnt加到512時(shí)，start_1置1。 在此之后，cnt開始計(jì)數(shù)。當(dāng)cnt也計(jì)到512時(shí)，start_2置1。 start_1與start_2相差512個(gè)時(shí)鐘周期，實(shí)現(xiàn)了上下支路的重疊加窗。

2.2 FIFO的作用

fft_fifo1和fft_fifo2作用是實(shí)現(xiàn)跨時(shí)鐘域的數(shù)據(jù)轉(zhuǎn)換以及數(shù)據(jù)緩存。如圖2所示，信號流經(jīng)了若干個(gè)時(shí)鐘域（10 MHz和200 MHz）。為了避免異步時(shí)鐘域產(chǎn)生錯(cuò)誤的采樣電平，使用FIFO緩存的方法完成異步時(shí)鐘域的數(shù)據(jù)轉(zhuǎn)換。

2.3 傅里葉變換的實(shí)現(xiàn)

FFT1024和IFFT1024使用IP核XFFT （Xilinx Fast Fourier Transform）v7.1，該模塊可以計(jì)算正向及逆向復(fù)數(shù)定點(diǎn)或浮點(diǎn)型FFT，并且可以實(shí)時(shí)配置。XFFT支持4種結(jié)構(gòu)：

1）Pipelined，Streaming I/O：允許連續(xù)的數(shù)據(jù)處理，速度快，缺點(diǎn)是占用資源大；

2）Radix－4，Burst I/O，數(shù)據(jù)與計(jì)算分別進(jìn)行，占用資源較少；

3）Radix－2，Burst I/O，與Radix－4相仿，但是占用的資源更少；

4）Radix－2 Lite，Burst I/O，該結(jié)構(gòu)占用資源最少，但是速度也最慢。

設(shè)計(jì)中采用了Radix－2，Burst I/O模式。這種模式速度較快，資源占用較少。

2.4 干擾判決及抑制的實(shí)現(xiàn)

2.5 數(shù)據(jù)輸出

對于上支路，fft_fifo2的寫入的是地址為512及其以后的數(shù)據(jù)。對于下支路，fft_fifo2的fifo2_we就是IFFT1024的ifft_dv。即IFFT1024只要一有輸出，就立即寫入到fft_fifo2中。這樣就實(shí)現(xiàn)了上下兩支路的重疊部分相加作為最終輸出。

3 測試結(jié)果

采用注入式測試的方法，對窄帶干擾抑制信號處理板的抗干擾性能進(jìn)行測試。測試框圖如圖4所示。

圖4 測試框圖Fig.4 Block diagram for test

信號源注入窄帶干擾信號，經(jīng)功分器送到數(shù)字信號處理板?？刂崎_關(guān)用來控制數(shù)字信號處理板是否工作，以比較窄帶干擾抑制效果。用頻譜儀記錄窄帶干擾抑制前和窄帶干擾抑制后兩種情況下自適應(yīng)信號處理模塊輸出的干擾信號功率。測試儀器為R&S公司的FSP－40頻譜分析儀。測試結(jié)果如圖5和圖6所示。

圖5 點(diǎn)頻干擾抑制效果對比圖Fig.5 Spectrum of narrowband interference suppression effect comparison

圖6 200kHz干擾抑制效果對比圖Fig.6 Spectrum of 200kHz interference suppression effect comparison

圖5（a）為5個(gè)點(diǎn)頻窄帶干擾的頻譜，圖5（b）為窄帶干擾抑制后的頻譜。圖6（a）為200 kHz干擾的頻譜，圖6（b）為窄帶干擾抑制后的頻譜。從測試結(jié)果可以看出，多個(gè)點(diǎn)頻的窄帶干擾和部分頻帶干擾都得到了很好的抑制，輸出信號的信干比大大提高。

4 結(jié)束語

文中提出了一種應(yīng)用于GPS接收機(jī)[6－7]的窄帶干擾抑制算法的實(shí)現(xiàn)方案，該方案以FPGA為實(shí)現(xiàn)平臺。實(shí)驗(yàn)結(jié)果表明，該方案可以為GPS接收機(jī)提供至少60 dB的抗窄帶干擾能力。與傳統(tǒng)窄帶干擾方法相比，收斂速度和信噪比損耗均得到了改善。

[1]孫希延，紀(jì)元法，施滸立，等.基于GPS軟件接收機(jī)的抗窄帶干擾方法[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2009，34（1）：77－80.SUN Xi-yan，JI Yuan-fa，SHI Hu-li，et al.Algorithms for removing narrow-band interference based on software GPS receiver[J].Geomatics and Information Science of Wuhan University，2009，34（1）：77－80.

[2]李然，趙剛.FFT重疊相加抗窄帶干擾的算法研究[J].四川理工學(xué)院學(xué)報(bào)：自然科學(xué)版，2009，22（3）：87－89.LI Ran，ZHAO Gang.Study on the narrowband interference suppression algorithm using lapped FFT[J].Journal of Sichuan University of Science&Engineering：Natural Science Edition，2009，22（3）：87－89.

[3]曾祥華，李崢嶸，王飛雪.擴(kuò)頻系統(tǒng)頻域窄帶干擾抑制算法加窗損耗研究[J].電子與信息學(xué)報(bào)，2004，26（8）：1276－1281.ZENG Xiang-hua，LI Zheng-rong，WANG Fei-xue.Study on windowing degradation of frequency domain narrowband interference suppression algorithms in spread spectrum system[J].Journal of Electronics&Information Technology，2004，26（8）：1276－1281.

[4]穆立波，霍志勇.頻域抗干擾中的門限檢測和處理[J].計(jì)算機(jī)與網(wǎng)絡(luò)，2010（10）：46－48.MU Li-bo，HUO Zhi-yong.Research on threshold detection and interference processiong in frequency domain antijamming[J].Computer&Network，2010（10）：46－48.

[5]鄒寧，徐松濤，劉明園，等.一種重疊加窗頻域抑制窄帶干擾算法及研究[J].現(xiàn)代防御技術(shù)，2010，38（3）：120－125.ZOU Ning，XU Song-tao，LIU Ming-yuan，et al.An overlap window-based frequency-domain narrowband interference suppression algorithm and its analysis[J].Modern Defence Technology，2010，38（3）：120－125.

[6]石超，孫保華，魏云飛，等.GPS接收機(jī)低噪聲放大器設(shè)計(jì)[J].電子科技，2010（3）：56－58.SHI Chao，SUN Bao-hua，WEI Yun-fei，et al.Design of a GPS receiver low noise amplifier[J].Electronic Science and Technology，2010（3）：56－58.

[7]史陽春，劉佑寶.基于模糊控制的高靈敏GPS接收機(jī)設(shè)計(jì)[J].現(xiàn)代電子技術(shù)，2012，35（9）：48－51，56.SHI Yang-chun，LIU You-bao.Design of high-sensitive GPS receiver based on fuzzy control[J].Modern Electronics Technique，2012，35（9）：48－51，56.