白 燕，盧曉春，王 瑾，韓 濤

(1.中國(guó)科學(xué)院 國(guó)家授時(shí)中心，西安 710600；2.中國(guó)科學(xué)院 精密導(dǎo)航定位與定時(shí)技術(shù)重點(diǎn)實(shí)驗(yàn)室，西安 710600；

一種新的衛(wèi)星導(dǎo)航寬帶干擾抑制方法

白 燕1，2，3，盧曉春1，2，王 瑾1，2，3，韓 濤1，2，3

(1.中國(guó)科學(xué)院 國(guó)家授時(shí)中心，西安 710600；2.中國(guó)科學(xué)院 精密導(dǎo)航定位與定時(shí)技術(shù)重點(diǎn)實(shí)驗(yàn)室，西安 710600；

3.中國(guó)科學(xué)院大學(xué)，北京 100039)

寬帶干擾是衛(wèi)星導(dǎo)航系統(tǒng)中常見(jiàn)的一種壓制性干擾，也是危害較大的一種干擾方式。針對(duì)線性調(diào)頻寬帶干擾，以離散傅里葉變換為信號(hào)分析工具，并在此基礎(chǔ)上提出一種基于NEW-MDCFT的寬帶干擾抑制算法。詳細(xì)介紹了該算法的模型和實(shí)現(xiàn)方式，并針對(duì)疊加了線性調(diào)頻寬帶干擾信號(hào)的北斗系統(tǒng)B1信號(hào)，采用該算法進(jìn)行計(jì)算機(jī)仿真，給出干擾抑制前后信號(hào)的時(shí)域性能、頻域性能、捕獲性能以及跟蹤性能等特性。結(jié)果表明，基于NEW-MDCFT的抑制算法具有良好的線性調(diào)頻寬帶干擾抑制性能，可以使接收機(jī)在干擾環(huán)境下順利完成捕獲和跟蹤，大大增強(qiáng)了衛(wèi)星導(dǎo)航信號(hào)的抗干擾能力。

衛(wèi)星導(dǎo)航；線性調(diào)頻寬帶干擾；干擾抑制

1 引言

衛(wèi)星導(dǎo)航定位技術(shù)在現(xiàn)代軍事、民用領(lǐng)域中發(fā)揮越來(lái)越重要的作用，但由于衛(wèi)星導(dǎo)航信號(hào)自身信號(hào)較弱，容易受到各種干擾的影響，從而對(duì)接收信號(hào)的性能造成損失[1-2]。實(shí)際空間環(huán)境中的干擾可分為人為干擾和非人為干擾兩類。不同干擾類型有不同的特性，不同的特性對(duì)應(yīng)不同的干擾抑制方法。非人為干擾大多為來(lái)自自然界的干擾，如系統(tǒng)內(nèi)部的干擾、噪聲干擾、多徑干擾等，這些干擾是客觀存在，非敵意的。目前對(duì)于非人為干擾只能進(jìn)行削弱，而不能完全消除。對(duì)于人為干擾，是一種故意干擾，意在對(duì)敵方的通信實(shí)施干擾，達(dá)到破壞對(duì)方通信的目的，人為干擾可以利用信號(hào)處理的方法對(duì)其進(jìn)行消除或削弱。

一般來(lái)說(shuō)，對(duì)衛(wèi)星導(dǎo)航系統(tǒng)接收機(jī)威脅最大的人為干擾可能是壓制式干擾[3-4]。壓制式干擾采用一部干擾機(jī)來(lái)擾亂該地域出現(xiàn)的所有衛(wèi)星信號(hào)，這種干擾方式又可分為寬帶和窄帶干擾。寬帶干擾采用一般的時(shí)域或頻域的抑制技術(shù)，很難達(dá)到理想的抑制效果，而必須考慮采用一些其他的濾波方法[5-7]。本文主要討論衛(wèi)星導(dǎo)航、電子對(duì)抗等環(huán)境中常用的一種干擾樣式即線性調(diào)頻(line frequency modulation，LFM)干擾。LFM信號(hào)作為一種大時(shí)間-帶寬積的非平穩(wěn)信號(hào)[8-9]，其頻譜在整個(gè)帶寬內(nèi)比較平坦，而在帶外的衰減很快，能量在帶寬內(nèi)非常集中。當(dāng)干擾為L(zhǎng)FM信號(hào)且干擾帶寬與導(dǎo)航信號(hào)的帶寬相同或相差不多時(shí)，衛(wèi)星導(dǎo)航擴(kuò)頻系統(tǒng)的性能會(huì)急劇下降，采用通常的自適應(yīng)濾波算法此時(shí)很難達(dá)到有效干擾抑制效果。針對(duì)這種情況，本文提出了一種基于NEW-MDCFT的抑制方法，該方法以離散傅里葉變換為基礎(chǔ)，可以有效地對(duì)LFM寬帶干擾進(jìn)行抑制。

2 離散傅里葉變換基礎(chǔ)

離散Chirp-Fourier變換[10-12](discrete Chirp-Fourier transform，DCFT)是一種有效的線性調(diào)頻信號(hào)檢測(cè)技術(shù)，它是離散Fourier 變換(discrete Fourier transform，DFT)的推廣，特別適合于LFM信號(hào)。

假設(shè)x(n)(0≤n≤N-1)為N點(diǎn)的離散LFM信號(hào)，定義DCFT為[13]

(1)

(2)

文獻(xiàn)[13]證明，若采樣點(diǎn)N為質(zhì)數(shù)且k0、l0為整數(shù)(0≤k0,l0≤N-1)，x(n)的DCFT幅值為

(3)

(4)

考慮式(5)形式的模擬Chirp信號(hào)

xa(t)=exp(j2π(α0t+β0t2))

(5)

設(shè)總采樣點(diǎn)數(shù)為N(N為質(zhì)數(shù))，采樣時(shí)刻滿足t=nT/N1/3，n=0,1,2,…,N-1， 則總采樣時(shí)間為N2/3T， 采樣后的離散線性調(diào)頻類(Chirp-like)信號(hào)可以變?yōu)?/p>

(6)

式(6)中，k0=a0TN2/3，k0=β0TN1/3。DCFT匹配LFM信號(hào)需要滿足兩個(gè)約束條件：1)信號(hào)長(zhǎng)度必須為質(zhì)數(shù)；2)離散LFM信號(hào)的參數(shù)k0、l0必須是整數(shù)或極為接近整數(shù)。

兩個(gè)約束條件中條件1比較容易滿足，但條件2的要求極為苛刻。因?yàn)樵诮o定質(zhì)數(shù)采樣點(diǎn)N的情況下，N1/3和N2/3必定不是整數(shù)，它們分別與參數(shù)a0T和β0T相乘幾乎不可能是整數(shù)。更為嚴(yán)重的是，即使k0和l0偏離整數(shù)一個(gè)不大的數(shù)據(jù)，定義的DCFT對(duì)Chirp信號(hào)的檢測(cè)性能將急劇下降。

為了克服DCFT的缺點(diǎn)，對(duì)DCFT的定義進(jìn)行了修正，它對(duì)信號(hào)的采樣點(diǎn)數(shù)和信號(hào)參數(shù)不附帶任何約束條件，適應(yīng)對(duì)離散LFM信號(hào)進(jìn)行檢測(cè)和參數(shù)估計(jì)，即提出修正離散Chirp-Fourier變換(modifiedDCFT，MDCFT)[14]。

定義MDCFT公式為

(7)

MDCFT反變換(inverseMDCFT，IMDCFT)公式為

(8)

對(duì)模擬LFM進(jìn)行N點(diǎn)采樣(N可以為任意整數(shù))，采樣時(shí)刻為

采樣后的離散LFM信號(hào)為

(9)

式(9)，中k0=α0T,l0=β0T2， 這里k0、l0不一定為整數(shù)。

MDCFT與DCFT相比，有如下幾點(diǎn)特點(diǎn)：

1)MDCFT對(duì)信號(hào)的總采樣點(diǎn)數(shù)和信號(hào)參數(shù)沒(méi)有任何約束條件，它對(duì)信號(hào)參數(shù)的精確估計(jì)范圍也遠(yuǎn)大于DCFT，所以更有意義。

2)MDCFT對(duì)LFM信號(hào)的能量聚集性與DCFT相比略有降低。

3)DCFT和MDCFT均為線性變換，在存在多個(gè)LFM信號(hào)分量時(shí)不存在交叉項(xiàng)，這是它們的主要優(yōu)點(diǎn)。

在理想情況下，當(dāng)采樣時(shí)鐘沒(méi)有偏差時(shí)，k0、l0都是整數(shù)或接近整數(shù)，MDCFT能夠匹配k0和l0為整數(shù)的情況，干擾抑制算法能夠很好的抑制LFM干擾；但當(dāng)接收端采樣時(shí)鐘不準(zhǔn)確時(shí)，LFM干擾的k0和l0有分?jǐn)?shù)部分，這時(shí)會(huì)出現(xiàn)以下問(wèn)題：

1)MDCFT會(huì)出現(xiàn)能量泄露，從而不能準(zhǔn)確檢測(cè)實(shí)際幅度的最大值。

2)得到的干擾信號(hào)的參數(shù)估計(jì)值出現(xiàn)較大偏差，從而影響干擾抑制的性能。

針對(duì)出現(xiàn)的問(wèn)題1及問(wèn)題2，為了匹配k0和l0的小數(shù)部分，本文采用一種新穎的MDCFT變換算法，即NEW-MDCFT變換

(10)

式(10)中，δk,δl∈[-0.5,0.5]。該定義能有效地匹配由采樣時(shí)鐘偏差所引起的問(wèn)題，更好地匹配LFM干擾信號(hào)。

對(duì)于如下形式的離散Chirp信號(hào)

(11)

當(dāng)k0、l0為整數(shù)且0≤k0,l0≤N-1，x(n)的NEW-MDCFT幅值為

(12)

當(dāng)k0、l0含有小數(shù)部分時(shí)，分別改變?chǔ)膋、δl可有效地補(bǔ)償小數(shù)部分，此時(shí)NEW-MDCFT的幅值與式(11)相同。

3 基于NEW-MDCFT的干擾抑制方法

3.1 干擾抑制模型

基于NEW-MDCFT的干擾抑制模型原理如圖1所示。基本思想是對(duì)接收信號(hào)r(n)進(jìn)行NEW-MDCFT變換，形成(k,l)二維平面，在此平面上進(jìn)行二維搜索得到的坐標(biāo)位置即為(k,l)的估計(jì)值。由于NEW-MDCFT變換保留了干擾的幅度和相位特性，因此可將峰值點(diǎn)的幅度經(jīng)過(guò)變換作為干擾信號(hào)幅度的估計(jì)值，峰值點(diǎn)的相位直接作為相位的估計(jì)值。根據(jù)得到的估計(jì)參數(shù)重構(gòu)干擾j(n)， 然后從接收信號(hào)中減去重構(gòu)后的干擾信號(hào)得到干擾抑制后的信號(hào)y(n)。

圖1 基于NEW-MDCFT的抑制模型

3.2 干擾抑制的實(shí)現(xiàn)步驟

1)粗略估計(jì)干擾的參數(shù)

2)精確估計(jì)干擾的參數(shù)(a0i,φ0i,k0i,l0i),i=1,2,…,M。

(12)

3)根據(jù)步驟2)中所得到的估計(jì)的各個(gè)干擾成分的參數(shù)重構(gòu)干擾，得到干擾j(n)的估計(jì)為

(13)

4)從接收信號(hào)r(n)中減去綜合重構(gòu)的干擾信號(hào)j(n)， 即得到干擾抑制后的信號(hào)y(n)， 然后就可以對(duì)y(n)進(jìn)行解擴(kuò)和判決等其他處理。

4 干擾抑制性能仿真

為了驗(yàn)證上述干擾抑制算法對(duì)LFM干擾的抑制性能，利用國(guó)家授時(shí)中心的全球衛(wèi)星導(dǎo)航系統(tǒng)(globalnavigationsatellitesystem，GNSS)信號(hào)質(zhì)量監(jiān)測(cè)采集系統(tǒng)采集到一組我國(guó)北斗衛(wèi)星導(dǎo)航系統(tǒng)簡(jiǎn)稱北斗系統(tǒng)(BeiDounavigationsatellitesystem，BDS)B1信號(hào)，并疊加一定能量的線性調(diào)頻干擾信號(hào)，干擾的帶寬覆蓋整個(gè)導(dǎo)航信號(hào)的帶寬。采用基于NEW-MDCFT的抑制算法對(duì)混有干擾的實(shí)測(cè)導(dǎo)航信號(hào)進(jìn)行干擾抑制，給出干擾抑制前后的信號(hào)時(shí)域性能、頻域性能圖、捕獲性能、以及跟蹤性能，分析該干擾抑制算法的性能。

4.1 時(shí)域性能仿真結(jié)果

當(dāng)對(duì)BDS B1信號(hào)混疊進(jìn)入干信比為-10 dB的線性調(diào)頻干擾后，信號(hào)的時(shí)域圖如圖2所示，可以看出，BDS衛(wèi)星導(dǎo)航信號(hào)疊加了線性調(diào)頻干擾后，其波形存在明顯的干擾雜波，而經(jīng)過(guò)基于NEW-MDCFT的干擾抑制后如圖3，信號(hào)波形變得較為純凈。

圖2 干擾抑制前時(shí)域性能圖

圖3 干擾抑制后時(shí)域性能

4.2 頻域性能仿真結(jié)果

從圖4可以看出，當(dāng)疊加了LFM干擾后，干擾的能量已經(jīng)完全覆蓋了衛(wèi)星導(dǎo)航信號(hào)，而經(jīng)過(guò)干擾抑制后，如圖5，可以明顯看出北斗衛(wèi)星導(dǎo)航信號(hào)的頻譜，其中心頻點(diǎn)在37.5 kHz的中頻上，在有用信號(hào)帶寬外，存在部分其他的諧雜波分量。

圖4 干擾抑制前頻域性能圖

圖5 干擾抑制后頻域性能圖

4.3 捕獲性能仿真結(jié)果

當(dāng)將-10 dB干信比的線性調(diào)頻擾疊加在衛(wèi)星導(dǎo)航信號(hào)上，此時(shí)信號(hào)能量被干擾所覆蓋，導(dǎo)致接收機(jī)無(wú)法識(shí)別衛(wèi)星，從而無(wú)法實(shí)現(xiàn)導(dǎo)航信號(hào)的捕獲跟蹤。

而通過(guò)基于NEW-MDCFT算法的干擾抑制后，信號(hào)可以順利完成捕獲和跟蹤，如圖6所示。

圖6 干擾抑制后捕獲性能圖

4.4 跟蹤性能仿真結(jié)果

如圖7所示，經(jīng)過(guò)基于NEW-MDCFT算法的干擾抑制后，順利完成跟蹤。

圖7 干擾抑制后跟蹤性能圖

4.5 仿真分析

從上述仿真結(jié)果可以看出，當(dāng)沒(méi)有采用干擾抑制算法時(shí)，衛(wèi)星導(dǎo)航信號(hào)的時(shí)域圖和頻譜圖明顯的受到嚴(yán)重線性調(diào)頻干擾的影響，并且無(wú)法識(shí)別衛(wèi)星，無(wú)法完成導(dǎo)航信號(hào)的捕獲和跟蹤。而經(jīng)過(guò)干擾抑制算法之后，信號(hào)受到的LFM干擾基本被濾除，且順利地完成了信號(hào)的捕獲和跟蹤程序。因此，基于NEW-MDCFT的干擾抑制算法可以較好的對(duì)線性調(diào)頻干擾進(jìn)行抑制。

5 結(jié)束語(yǔ)

研究衛(wèi)星導(dǎo)航系統(tǒng)接收機(jī)抗干擾技術(shù)，對(duì)于提高衛(wèi)星導(dǎo)航系統(tǒng)在復(fù)雜電磁環(huán)境下的工作性能、增強(qiáng)其在各種環(huán)境中的適用性具有重要意義。本文以衛(wèi)星導(dǎo)航系統(tǒng)中常見(jiàn)的寬帶線形調(diào)頻干擾為研究對(duì)象，以離散Chirp-Fourier變換為信號(hào)分析工具，并在此基礎(chǔ)上提出基于NEW-MDCFT的干擾抑制算法。通過(guò)對(duì)北斗導(dǎo)航衛(wèi)星信號(hào)疊加線性調(diào)頻寬帶干擾后的混合信號(hào)進(jìn)行計(jì)算機(jī)仿真，比較該干擾抑制算法在干擾抑制前后對(duì)導(dǎo)航信號(hào)的時(shí)域、頻域、捕獲、跟蹤性能等的影響。結(jié)果表明，該干擾抑制算法具有較好的干擾抑制性能。

[1] MAY M M，LAINEZ M D，MARTIN J R，et al.GNSS，an Evolving Technology.Current Systems Evolutions and Future Perspectives[C]//Proceedings of the 24th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2011).Portland，OR：The Institute of Navigation，2011：3083-3101.

[2] 劉基余.GNSS 全球?qū)Ш叫l(wèi)星系統(tǒng)的發(fā)展與挑戰(zhàn)[EB/OL].[2014-05-23].http：//www.docin.com/p-117236353.html.

[3] KAPLAN E D，HEGARTY C J.GPS原理與應(yīng)用[M].寇艷紅，譯.2 版.北京：電子工業(yè)出版社，2007.

[4] 費(fèi)華連.導(dǎo)航戰(zhàn)中的GPS干擾與抗干擾技術(shù)[J].航空電子技術(shù)，2001，32(1)：19-22.

[5] 朱春華，穆曉敏.DSSS系統(tǒng)中線性調(diào)頻干擾抑制技術(shù)研究[J].電波科學(xué)學(xué)報(bào)，2003，18(3)：341-351.

[6] 劉艷良，孔軍輝.直接序列擴(kuò)頻通信系統(tǒng)LFM干擾抑制技術(shù)研究[J].電子信息對(duì)抗技術(shù)，2013，7(4)：56-60.

[7] LIU Yao，NING Peng，DAI Huai-yu，et al.Randomized Differential DSSS：Jamming-Resistant Wireless Broadeast Communication[C]//Proceedings of the 29th IEEE International Conference on Computer Communications (IEEE INFOCOM 2010).San Diego，CA：IEEE，2010：1-9.

[8] 朱健東，趙擁軍，唐江.線性調(diào)頻連續(xù)波信號(hào)的周期分?jǐn)?shù)階Fourier變換檢測(cè)與估計(jì)[J].電子與信息學(xué)報(bào)，2013，35(8)：1827-1833.

[9] 張賢達(dá)，保錚.非平穩(wěn)信號(hào)分析與處理[M].北京：國(guó)防工業(yè)出版社，1998.

[10]WEI Yong-mei，BI Guo-an.Broadband Interference Suppression In DS-SS System with Modified Chirp Fourier Transform[C]//Proceedings of the 9th IEEE International Conference on Communications Systems(ICCS 2004).Singapore：IEEE，2004：640-644.

[11]FAN Bo，ZHANG Kai，QIN Yu-liang，et al.Discrete Chirp-Fourier Transform-based Acquisition Algorithm for Weak Global Positioning System L5 Signals in High Dynamic Environments[J].Radar，Sonar & Navigation，IET，2013，7(7)：736-746.

[12]WU Liang，WEI Xi-zhang，YANG De-gui，et al.ISAR Imaging of Targets with Complex Motion Based on Discrete Chirp Fourier Transform for Cubic Chirps[J].IEEE Transactions on Geoscience and Remote Sensing，2012，50(10)：4201-4212.

[13]XIA Xiang-gen.Discrete Chirp-Fourier Transform and Its Application to Chirp Rate Estimation[J].IEEE Transactions on Signal Processing，2000，48(11)：3122-3133.

[14]孫泓波，郭欣，顧紅，等.修正離散Chirp Fourier變換在運(yùn)動(dòng)目標(biāo)檢測(cè)中的應(yīng)用[J].電子學(xué)報(bào)，2003，31(1)：1-4.

An Interference Suppression Method Based on NEW-MDCFT in Satellite Navigation System

BAI Yan1，2，3，LU Xiao-chun1，2，WANG Jin1，2，3，HAN Tao1，2，3

(1.National Time Service Center，Chinese Academy of Sciences，Xi’an 710600，China；2.Key Laboratory for Precise Navigation and Timing Technology，National Time Service Center，Chinese Academy of Sciences，Xi’an 710600，China；3.Universisty of Chinese Academy of Sciences，Beijing 100039，China)

Wideband interference is common oppressive which impact on the receiver performance badly in the satellite navigation system.In this paper，aiming at the linear frequency modulation (LFM) wideband interference，an interference suppression method based on NEW-MDCFT is presented，and the model and implementation method of the suppression algorithm have been introduced.By adding the LFM jamming to the COMPASS B1 signal，the interference suppression method is simulated，and time domain performance，frequency domain performance，acquisition performance and tracking performance of the signal before and after interference suppression have been given.The results indicate that，the suppression method based on NEW-MDCFT has good LFM interference suppression performance and is able to perform good interference suppression for seriously-interfered information navigation signal，so that the receiver can accomplish acquisition and tracking smoothly.

wideband interference；LFM；interference synthesis

2014-04-29

國(guó)家自然科學(xué)基金(11203027)。

白燕(1975)，女，陜西清澗縣人，副研究員，博士，主要從事衛(wèi)星導(dǎo)航技術(shù)及信號(hào)處理等方面的研究工作。

TN967

A

2095-4999(2015)-01-0061-05