王 哲, 閆 巖, 金 釗, 程思敏
(北京遙感設(shè)備研究所,北京100854)
調(diào)頻引信是一種發(fā)射信號頻率按調(diào)制信號規(guī)律變化的無線電引信,利用發(fā)射信號和回波信號頻率差值與彈目距離的對應(yīng)關(guān)系,實(shí)現(xiàn)引信定距功能[1]。相對于連續(xù)波多普勒引信,具有定距精度高、抗干擾性能好等特點(diǎn),且算法復(fù)雜度低,容易在硬件平臺上實(shí)現(xiàn)[2],在彈藥裝備中得到了廣泛應(yīng)用。
電子干擾對調(diào)頻引信的戰(zhàn)場生存能力構(gòu)成了嚴(yán)重威脅[3-4]。基于數(shù)字射頻存儲 (DRFM)的干擾機(jī)可以準(zhǔn)確地存儲、重構(gòu)和轉(zhuǎn)發(fā)調(diào)頻引信的發(fā)射信號,從而實(shí)現(xiàn)假目標(biāo)欺騙干擾。近年來,抗DRFM干擾研究的一個重要分支是發(fā)射波形的設(shè)計,包括脈沖壓縮[5]、頻率捷變[6]、脈沖多樣性[7]、脈沖重復(fù)周期跳變等方法[8]。頻率捷變系統(tǒng)可以在保留信號處理設(shè)計的基礎(chǔ)上,僅通過改變前端設(shè)計獲得良好的抗干擾效果,在各種引信體制中已得到了廣泛應(yīng)用。在工程實(shí)現(xiàn)方面,還需要克服高次諧波的影響。隨著跳頻點(diǎn)數(shù)的增加,系統(tǒng)信噪比會進(jìn)一步降低。
本文基于發(fā)射波形去周期化設(shè)計的抗干擾思路,提出一種基于雙調(diào)制頻率跳變的調(diào)頻引信抗DRFM干擾方法。這種方法可以改善現(xiàn)有調(diào)頻體制引信的距離模糊特性,從而抑制存儲轉(zhuǎn)發(fā)干擾的影響。本文以理論分析為基礎(chǔ),介紹了半實(shí)物仿真測試結(jié)果,驗(yàn)證了該方法的可行性與有效性。
如前所述,調(diào)頻引信發(fā)射信號的頻率是時間的函數(shù),發(fā)射信號通常采用周期調(diào)制,差頻信號頻譜以調(diào)制頻率整數(shù)倍的各次諧波形式分布,利用某次諧波能量的最大值與距離的對應(yīng)關(guān)系,可實(shí)現(xiàn)引信定距功能。圖1為典型單通道調(diào)頻多普勒引信原理框圖。
調(diào)制器產(chǎn)生頻率為fm的三角波調(diào)制信號對振蕩器進(jìn)行線性調(diào)頻,目標(biāo)反射信號被收發(fā)共用天線接收,經(jīng)環(huán)形器后與本振信號混頻得到差頻信號,邊帶放大器選取mfm±fD頻率成分進(jìn)行二次混頻和包絡(luò)檢波,得到幅度隨距離變化的多普勒信號包絡(luò)送至執(zhí)行級進(jìn)行起爆判決。

圖1 單通道調(diào)頻多普勒引信原理框圖
以對地連續(xù)波調(diào)頻引信為例,DRFM干擾模型示意圖如圖2所示。干擾機(jī)通過延遲轉(zhuǎn)發(fā)引信發(fā)射信號,在引信接收端形成距離假目標(biāo),導(dǎo)致引信提前啟動。

圖2 DRFM干擾模型示意圖
已知調(diào)頻引信發(fā)射信號為s(t),DRFM干擾機(jī)對偵收到的引信發(fā)射信號進(jìn)行存儲,并在附加延遲時間后轉(zhuǎn)發(fā),干擾信號sj(t)可以表示為

式中:Aj為DRFM干擾機(jī)發(fā)射干擾信號幅度;τj為干擾信號延遲時間,τj=2Rj/c+τj0,它由信號在空間中的傳播時間和干擾機(jī)附加延遲時間τj0組成;Rj為干擾機(jī)到引信的距離;c為光速。
目標(biāo)回波信號sr(t)可以表示為

式中:Aj為目標(biāo)回波信號的幅度;τd為目標(biāo)回波延遲時間,τd=2Rd/c;Rd為彈目距離,且滿足引信的起爆距離。根據(jù)DRFM干擾模型,一般情況下Rj遠(yuǎn)大于Rd,相應(yīng)的延遲時間τj也遠(yuǎn)大于τd。雖然干擾信號可以在形式上高度近似目標(biāo)回波信號,但τj在小于調(diào)制周期Tm內(nèi),無法導(dǎo)致引信接收端輸出滿足起爆距離的假目標(biāo)信號。因此DRFM干擾機(jī)成功干擾調(diào)頻引信的必要條件是τj至少要大于一個Tm,即滿足以下條件:

式中:N*為正整數(shù)。DRFM干擾機(jī)在獲取Rj的前提下,通過合理設(shè)置τj0,使得式(3)成立,從而利用s(t)的周期調(diào)制帶來的距離模糊實(shí)現(xiàn)欺騙干擾,使引信在距離目標(biāo)很遠(yuǎn)的位置探測到滿足起爆距離的假目標(biāo)而提前啟動。
根據(jù)上述失效機(jī)理,可以將發(fā)射波形去周期化設(shè)計作為調(diào)頻引信抗DRFM干擾研究的基本思路。本文提出一種雙調(diào)制頻率隨機(jī)跳變的三角波線性調(diào)頻發(fā)射波形,采用m序列作為控制fm跳變的偽隨機(jī)碼,實(shí)現(xiàn)發(fā)射波形去周期化,并設(shè)計基于FFT的諧波時序檢測方法完成目標(biāo)判別。圖3為抗DRFM干擾調(diào)頻引信系統(tǒng)原理框圖。

圖3 抗DRFM干擾調(diào)頻引信系統(tǒng)原理框圖
圖3 中調(diào)頻引信系統(tǒng)與傳統(tǒng)方案主要有兩點(diǎn)差異:其一,采用雙調(diào)制率調(diào)制器產(chǎn)生復(fù)合調(diào)制波形;其二,采用基于FFT的諧波時序檢測方法處理差頻信號。
已知傳統(tǒng)三角波線性調(diào)頻引信的差頻信號頻率fi與彈目距離R之間的關(guān)系由調(diào)頻測距公式確定,即

為了提高調(diào)頻引信抗DRFM干擾性能,可以改變調(diào)制參數(shù)實(shí)現(xiàn)發(fā)射波形的去周期化,可變參數(shù)包括載頻f0、調(diào)制頻率fm和半調(diào)制頻偏ΔFm。
通過改變f0可以在保持原有信號處理方法不變的同時,取得比較理想的抗干擾效果[9]。但發(fā)射信號頻率隨時間非連續(xù)變化,且存在大幅度跳變,對VCO的響應(yīng)性能要求很高。同時,由于差頻信號在時頻域存在較大“毛刺”,進(jìn)而會帶來高次諧波,增加了信號處理的難度。
采用變fm實(shí)現(xiàn)發(fā)射波形去周期化,由式(6)可知,單一改變fm或ΔFm會導(dǎo)致fi與R不再滿足線性對應(yīng)關(guān)系,給引信定距帶來不便。本文設(shè)計ΔFm在fm跳變的同時也發(fā)生跳變,但調(diào)制斜率β=4ΔFmfm始終不變,從而保持fi與R的線性對應(yīng)關(guān)系。
圖4為本文設(shè)計雙調(diào)制頻率三角波線性調(diào)頻發(fā)射信號頻率隨時間變化曲線,調(diào)制周期Tm=1/fm在兩種取值之間隨機(jī)跳變,且Tm1=2Tm2;ΔFm也存在兩種取值,且滿足ΔFm1=2ΔFm2。調(diào)制參數(shù)的跳變由二進(jìn)制m序列生成的偽隨機(jī)碼控制。為便于工程實(shí)現(xiàn),取碼元寬度Tc=Tm1,偽隨機(jī)碼周期Tr=PTc,P為偽隨機(jī)序列長度。假設(shè)調(diào)頻信號初始相位φ0=0,載頻為f0,信號幅度為At,則發(fā)射信號在一個調(diào)制周期[-Tm1/2,+Tm1/2]內(nèi)可以表示為


圖4 雙調(diào)制頻率三角波線性調(diào)頻發(fā)射信號時頻圖
三角波線性調(diào)頻信號u1(t)可以表示為

同理,可得到Tm2周期內(nèi)發(fā)射信號及其對應(yīng)的調(diào)頻信號表達(dá)式分別為v2(t)、u2(t)。
m序列產(chǎn)生器一般的數(shù)學(xué)表達(dá)式為

式中:Ci={1,0}為二進(jìn)制m序列,通過N位線性移位寄存器產(chǎn)生,偽隨機(jī)序列長度P=2N-1。假設(shè)Ci,n表示對m序列Ci進(jìn)行邏輯非運(yùn)算的結(jié)果,代入可得到其反碼表達(dá)式為pn(t)。
綜上所述,可得雙調(diào)制頻率三角波線性調(diào)頻發(fā)射信號為

由式(8)可知,發(fā)射信號由調(diào)制參數(shù)不同的兩種調(diào)頻信號組成,分別在m序列及其反碼的控制下隨機(jī)跳變。式中兩項(xiàng)調(diào)頻信號存在半個周期的平移,目的是使疊加的發(fā)射信號瞬時頻率隨時間連續(xù)變化。
雙調(diào)制頻率調(diào)頻引信差頻信號的頻率隨時間變化關(guān)系如圖5,圖中上半部分的實(shí)線表示發(fā)射信號,虛線表示目標(biāo)回波信號;τ=2R/c表示目標(biāo)回波延遲時間。下半部分為差頻信號頻率fi隨時間變化的曲線。

圖5 目標(biāo)回波作用下差頻信號時頻圖
由圖5可見,由于β=4ΔFmfm不變,fi的頻率成分集中且與R線性對應(yīng)。不同fm差頻信號的不規(guī)則區(qū)具有不同的分布特點(diǎn),且fi變化的周期等于Tm,即fi變化規(guī)律也在m序列的控制下,在兩種狀態(tài)之間隨機(jī)跳變。首先對圖5中碼元Tc1對應(yīng)的Tm1內(nèi)差頻信號作傅里葉變換,得到頻率間隔為fm1的離散譜表達(dá)式[10]為

式中:m 為諧波次數(shù);α1(m,τ)為傅里葉系數(shù)表達(dá)式,具體形式如式(10)所示。差頻信號中m次諧波傅里葉系數(shù)α1(m,τ)包含兩個sinc函數(shù)項(xiàng),其中第二項(xiàng)的自變量取值范圍遠(yuǎn)離包絡(luò)主瓣而可以忽略。差頻信號的m次諧波幅度隨延遲時間τ變化具有類似sinc函數(shù)包絡(luò)的形式,主瓣寬度 Δτ1=1/ΔFm1,且在 m=8ΔFm1R/c時出現(xiàn)峰值。

對于圖5,Tc2內(nèi)的差頻信號只要將式(9)與式(10)中的Tm1、ΔFm1替換成Tm2、ΔFm2即可得到相應(yīng)表達(dá)式,且差頻信號的m次諧波幅度隨延遲時間τ的變化規(guī)律相同,差異在于選定的m次諧波頻率mfm2=2mfm1,且sinc函數(shù)包絡(luò)主瓣寬度Δτ2=2Δτ1。
下面討論彈目存在相對運(yùn)動時,諧波幅度隨時間的變化關(guān)系。假設(shè)引信相對于目標(biāo)作勻速直線運(yùn)動且相對速度為vR,彈目初始距離為R0,則目標(biāo)回波延遲時間為τ=2(R0-vRt)/c,代入式(10)中得到Tc1、Tc2內(nèi)m次諧波幅度隨時間變化的表達(dá)式分別為α1,D(m,t)、α2,D(m,t)。相當(dāng)于對原式進(jìn)行自變量的平移和尺度變換,因此包絡(luò)仍具備sinc函數(shù)的形式。另外由于存在相對運(yùn)動,各次諧波幅度會受到多普勒頻率fD的幅度調(diào)制[11]。結(jié)合m序列偽隨機(jī)碼及其反碼的表達(dá)式p(t)、pn(t),可以得到彈目接近過程中,差頻信號中m次諧波幅度隨時間變化的表達(dá)式為
圖6為差頻信號m次諧波幅度α(m,t)隨時間變化的關(guān)系,近似sinc函數(shù)包絡(luò),不規(guī)則部分是由兩種主瓣寬度不同的sinc包絡(luò)合成導(dǎo)致的。隨著彈目距離不斷減小,從高次到低次諧波信號會依次出現(xiàn)峰值,可以利用各次諧波的時序性作為目標(biāo)判別依據(jù)。

圖6 m次諧波幅度隨時間的變化
圖7 為DRFM干擾作用下差頻信號頻率隨時間變化曲線,圖中表示干擾信號延遲一個Tm1的干擾效果。可以看出由于發(fā)射波形采用了去周期化設(shè)計,相對于目標(biāo)回波,DRFM干擾作用下的差頻信號頻率成分更加豐富,變化范圍增大,不再具有明顯的規(guī)律性。差頻信號諧波能量分散,無法輸出諧波幅度隨距離變化的完整包絡(luò),干擾效果會受到抑制。

圖7 DRFM干擾下差頻信號時頻圖
圖8 為基于FFT的雙通道諧波時序檢測方法原理圖,首先通過ADC以采樣頻率fs對差頻信號進(jìn)行N點(diǎn)采樣。然后,通過N點(diǎn)的FFT運(yùn)算獲得差頻信號離散頻譜,提取頻譜中各次諧波頻率對應(yīng)的譜線幅值作為各通道的諧波幅值。以N點(diǎn)為周期重復(fù)上述過程,即可得到各次諧波幅值序列,經(jīng)檢波后輸出各次諧波幅度隨距離變化的包絡(luò)進(jìn)行啟動判決。

圖8 雙通道諧波時序檢測方法原理框圖
該方法中FFT點(diǎn)數(shù)的選取應(yīng)考慮兩方面因素。首先,在差頻信號分析中,通過解析單一Tm內(nèi)頻譜得到各次諧波幅度隨距離變化關(guān)系,進(jìn)行FFT運(yùn)算的采樣時長應(yīng)至少大于一個Tm。其次,考慮彈目交會過程中產(chǎn)生的多普勒效應(yīng)作用于差頻信號,相當(dāng)于抑制載波的幅度調(diào)制,原有各次諧波譜線會被頻率為mfm±fD兩根譜線取代。為了減小多普勒效應(yīng)的影響,設(shè)計FFT的頻率分辨率fs/N 應(yīng)遠(yuǎn)大于多普勒頻率fD,從而利用FFT的泄漏使分散到mfm±fD的兩根譜線上的能量集中到mfm諧波的譜線上。綜上所述,F(xiàn)FT點(diǎn)數(shù)N的選取應(yīng)滿足條件:

為了驗(yàn)證本文設(shè)計的調(diào)頻引信信號處理系統(tǒng)的定距與抗干擾性能,采用Simulink仿真模型生成的差頻信號作為激勵信號,通過Modelsim仿真環(huán)境觀測系統(tǒng)各級輸出信號。仿真參數(shù):f0=3 GHz,fm1=100 k Hz,fm2=200 k Hz,ΔFm1=100 MHz,ΔFm2=50 MHz,偽隨機(jī)序列長度P=4 095,碼元寬度Tc=10μs。采樣頻率5 MHz,F(xiàn)FT點(diǎn)數(shù)64,選取400 k Hz和800 k Hz諧波通道,對應(yīng)定距距離為3 m和6 m。彈目距離變化范圍15 m~0 m,彈體飛行速度500 m/s,仿真時間30 ms,信噪比-10 d B。

圖9 調(diào)頻引信系統(tǒng)Modelsim仿真結(jié)果
圖9 (a)為目標(biāo)回波作用下系統(tǒng)仿真結(jié)果,雙通道諧波包絡(luò)幅值大于幅度門限后輸出方波信號,包絡(luò)峰值對應(yīng)距離分別為5.5 m和2.5 m。時序邏輯模塊判斷雙通道輸出方波信號的出現(xiàn)次序?yàn)楦叽沃C波先于低次諧波,且間隔時間滿足彈目接近速度變化范圍,在彈目距離2 m附近輸出啟動信號。
圖9(b)為DRFM干擾作用下的系統(tǒng)仿真結(jié)果,雙通道諧波幅值明顯小于目標(biāo)回波,并且隨時間分布規(guī)律相對雜散。雙通道諧波包絡(luò)信號峰值均小于幅度門限,無幅度判決方波信號輸出,不滿足啟動判決條件從而無啟動信號輸出。
利用現(xiàn)場可編程邏輯門陣列(FPGA)作為主芯片,設(shè)計制作引信信號處理電路,對該方案進(jìn)行硬件驗(yàn)證。采用半實(shí)物仿真的方法,利用函數(shù)發(fā)生器分別產(chǎn)生目標(biāo)回波與DRFM干擾作用下差頻信號,作為信號處理電路的中頻輸入,通過示波器觀測DAC輸出的雙通道包絡(luò)檢波信號。系統(tǒng)半實(shí)物仿真測試結(jié)果如圖10所示。

圖10 半實(shí)物仿真測試檢波包絡(luò)
可以看出,測試結(jié)果與仿真結(jié)果基本一致。目標(biāo)回波作用下的雙通道諧波包絡(luò)完整,信號幅度峰值大于100 mV,且滿足高次諧波峰值先于低次出現(xiàn)的時序關(guān)系,雙通道諧波峰值時間間隔約為6 ms,對應(yīng)的距離變化為3 m,經(jīng)幅度與時序邏輯判決后可以輸出啟動信號。DRFM干擾作用下諧波包絡(luò)與目標(biāo)作用下差異明顯,且信號幅度均小于50 mV,無法滿足幅度判決條件。
設(shè)置DRFM干擾信號的初始延遲時間在0.975 Tm1~1.025 Tm1之間的敏感范圍內(nèi),以10 ns為延時步進(jìn),進(jìn)行50組抗干擾半實(shí)物仿真測試,雙調(diào)制頻率調(diào)頻引信被干擾6次,抗干擾成功率為88%;傳統(tǒng)調(diào)頻引信被干擾27次,抗干擾成功率為46%。測試結(jié)果表明,雙調(diào)制頻率調(diào)頻引信系統(tǒng)可以在滿足定距性能的前提下,提高調(diào)頻引信抗DRFM干擾的能力。
本文在分析傳統(tǒng)調(diào)頻引信在DRFM干擾下失效機(jī)理的基礎(chǔ)上,設(shè)計了以雙調(diào)制頻率隨機(jī)跳變的三角波線性調(diào)頻信號作為發(fā)射波形的調(diào)頻引信抗DRFM干擾方法,提出了基于FFT的諧波時序檢測方法判別目標(biāo),并通過Modelsim仿真分析與硬件測試驗(yàn)證了該抗干擾方法可以在滿足系統(tǒng)定距性能的前提下,提高調(diào)頻引信的抗DRFM干擾性能。