張紅旗　李曉　沈欣

摘要：? ? ? 為了提高引信的最佳毀傷控制能力， 本文結(jié)合LFMCW引信特點(diǎn)對差拍信號進(jìn)行了理論分析， 得出了影響距離分辨率、 速度分辨率等指標(biāo)的因素， 并利用Simulink軟件搭建了相關(guān)的仿真模型， 對通過二維傅里葉變換獲取距離、 速度信息的方法進(jìn)行仿真驗(yàn)證。 該模型可以模擬引信回波信號、 噪聲信號的產(chǎn)生以及中頻信號處理的過程。 仿真結(jié)果驗(yàn)證了算法的可行性與模型的準(zhǔn)確性， 為引信系統(tǒng)優(yōu)化設(shè)計奠定了基礎(chǔ)。

關(guān)鍵詞：? ? ?LFMCW引信; Simulink; 測距; 測速; 二維傅里葉變換

中圖分類號：? ? ? TJ43+4.1文獻(xiàn)標(biāo)識碼：? ? A文章編號：? ? ?1673-5048（2020）02-0053-06

0引言

未來空中威脅目標(biāo)性能的提高、 作戰(zhàn)要求的變化和空戰(zhàn)環(huán)境的復(fù)雜化對空空導(dǎo)彈提出了嚴(yán)峻的挑戰(zhàn)[1]。 作為導(dǎo)彈的重要組成部分， 無線電引信要在多目標(biāo)、 全空域等復(fù)雜環(huán)境下完成對目標(biāo)的精確毀傷[2]， 最佳毀傷控制能力成為引信的發(fā)展方向[3]， 需要盡可能多地利用距離、 速度等多個維度的信息。

線性調(diào)頻連續(xù)波（Linear Frequency Modulation Continuous Wave， LFMCW）體制具有比較大的時寬帶寬積， 相比脈沖體制， 具備較高的測距精度[4]、 較低的發(fā)射功率、 高距離分辨力和無距離盲區(qū)等特點(diǎn)[5]， 在近距離、 高分辨的探測場合得到越來越廣泛的應(yīng)用[6]， 成為無線電引信的重要發(fā)展方向之一。

LFMCW體制是通過發(fā)射信號和接收信號的頻率差進(jìn)行距離測量， 運(yùn)動目標(biāo)的多普勒頻率會疊加到信號的頻率差中（即距離-速度耦合現(xiàn)象）[7]。 為提高測距精度， 文獻(xiàn)[8]采用了三角波調(diào)頻方式， 通過上升階段和下降階段頻率的不同來消除多普勒影響， 但存在多目標(biāo)時距離-速度匹配困難的問題。 本文在文獻(xiàn)[9]快速斜坡模式測距的基礎(chǔ)上， 結(jié)合引信的特點(diǎn)對差拍信號進(jìn)行理論分析， 并利用MATLAB/Simulink 建立了LFMCW體制引信的仿真模型， 所設(shè)計的系統(tǒng)能夠?qū)σ诺母鞑糠中盘栠M(jìn)行實(shí)時直觀觀察， 仿真結(jié)果能夠?yàn)樵擉w制引信工程設(shè)計提供完善的依據(jù)。

1Simulink簡介

MATLAB 是MathWorks公司推出的通用科學(xué)計算和系統(tǒng)仿真軟件， 被廣泛應(yīng)用于各類數(shù)值計算和仿真。

Simulink是基于MATLAB強(qiáng)大的可視化動態(tài)系統(tǒng)仿真工具， 其提供一個動態(tài)系統(tǒng)建模、 仿真和綜合分析的集成， 支持線性系統(tǒng)和非線性系統(tǒng)， 可以用連續(xù)采樣時間、 離散采樣時間或兩種混合的采樣時間進(jìn)行建模， 支持多速率系統(tǒng)， 也就是系統(tǒng)中的不同部分具有不同的采樣速率[10]。

3LFMCW引信的Simulink 仿真

3.1系統(tǒng)模型建立

LFMCW引信的Simulink仿真模型主要由調(diào)頻信號產(chǎn)生器、 目標(biāo)回波模擬器、 混頻濾波器、 低通濾波器、 采樣、 信號處理及波形顯示等部分組成，如圖2所示。

信號流程： 調(diào)頻信號產(chǎn)生器產(chǎn)生調(diào)頻周期T、 掃頻寬度B、 載頻f0的鋸齒波調(diào)頻信號， 作為引信發(fā)射信號; 目標(biāo)回波模擬器對引信發(fā)射信號進(jìn)行處理， 模擬距離R處相對速度為v的目標(biāo)回波; 目標(biāo)回波信號與引信發(fā)射信號進(jìn)行混頻， 并通過低通濾波器取出差拍信號， 經(jīng)零階保持器（ Zero-Order Hold） 后轉(zhuǎn)換為時間離散信號， 在信號處理器中進(jìn)行數(shù)據(jù)重排、 二維傅里葉變換， 通過波形顯示模塊觀察輸出波形。

調(diào)頻信號產(chǎn)生器模塊由Repeating Sequence模塊和Continuous-Time VCO模塊構(gòu)成， 如圖3所示。 該模塊可以對調(diào)頻周期、 掃頻寬度、 載頻、 幅度進(jìn)行調(diào)整。

目標(biāo)回波模擬器模塊主要圍繞可變延時模塊（Variable Time Delay）構(gòu)成， 圖4所示為單點(diǎn)目標(biāo)的目標(biāo)回波模擬器模塊。 通過多個模塊的并聯(lián)并設(shè)置不同模塊的速度、 距離、 幅度、 相位， 按照亮點(diǎn)法[14]進(jìn)行復(fù)雜形體目標(biāo)的回波模擬。

Random Source模塊產(chǎn)生服從高斯分布的隨機(jī)數(shù)噪聲信號[15]， 用以模擬接收通道中附加的高斯白噪聲。

信號處理模塊包括數(shù)據(jù)緩存、 數(shù)據(jù)重排、 數(shù)據(jù)加窗、 二維FFT運(yùn)算等， 見圖5。 其中數(shù)據(jù)緩存模塊（Buffer）對一幀的數(shù)據(jù)進(jìn)行緩存， 通過數(shù)據(jù)重排（Reshape）形成M×N的矩陣， 進(jìn)行加窗處理和二維FFT運(yùn)算， 從而獲得二維頻譜分布。 其中二維FFT可分解為兩次一維的FFT進(jìn)行。 在波形顯示模塊中調(diào)用MATLAB的Surfl函數(shù)進(jìn)行顯示。

3.2系統(tǒng)仿真結(jié)果

主要仿真參數(shù)如下： 調(diào)頻周期T=2 μs; 掃頻帶寬B=300 MHz; 載頻 f0=10 GHz; 彈目距離R=10 m; 相對速度v=-500 m/s; 采樣率fs =50 MHz; 第一維FFT點(diǎn)數(shù)M=100; 第二維FFT點(diǎn)數(shù)N=128。 仿真獲得的波形如圖6所示， 可見差拍信號為正弦波。

對每個掃頻周期內(nèi)的差拍信號進(jìn)行一維傅里葉變換， 可以獲得其頻譜， 見圖7， 其峰值位置在10.1 MHz處， 對應(yīng)的距離為cfm2k=cfm2B/T=3 m。

N個不同周期10 MHz處頻譜的波形見圖8。

可見波形為正弦波形， 其頻譜見圖9， 頻率為多普勒頻率， 體現(xiàn)出彈目相對速度為-531.5 m/s。

經(jīng)過對每個距離單元的信號進(jìn)行第二次傅里葉變換， 獲得差拍信號的二維分布見圖10。檢測到的目標(biāo)距離為10.1 m，? 相對速度為-531.5 m/s， 此時距離分辨率ΔR=0.5c/B=0.5 m， 速度分辨率Δv=0.5c/（f0 NT） =58.5 m/s， 可見與設(shè)定的目標(biāo)參數(shù)一致。

可見， 采用該仿真模型可以分辨出目標(biāo)1、 目標(biāo)2、 目標(biāo)3， 但由于目標(biāo)4、 目標(biāo)5的距離/速度均與目標(biāo)3接近， 因此無法區(qū)分出目標(biāo)3、 目標(biāo)4、 目標(biāo)5， 而將其認(rèn)為是一個目標(biāo)， 與理論分析結(jié)果一致， 說明模型建立的正確性。 在圖11的基礎(chǔ)上可以進(jìn)一步進(jìn)行二維的CFAR（Constant False Alarm Detection）， 進(jìn)一步提高引信的檢測性能。

4結(jié)論

本文在分析LFMCW引信回波信號的基礎(chǔ)上， 建立了LFMCW引信的Simulink動態(tài)仿真模型， 體現(xiàn)出Simulink在系統(tǒng)仿真方面的強(qiáng)大功能， 并對距離分辨率、 速度分辨率等關(guān)鍵參數(shù)進(jìn)行了分析， 通過分析引信的工作特性， 為實(shí)現(xiàn)引信系統(tǒng)結(jié)構(gòu)的進(jìn)一步完善提供了必要的理論依據(jù)， 是對引信系統(tǒng)進(jìn)行優(yōu)化設(shè)計的重要基礎(chǔ)。

參考文獻(xiàn)：

[1] 樊會濤， 張蓬蓬. 空空導(dǎo)彈面臨的挑戰(zhàn)[J]. 航空兵器， 2017（2）： 3-7.

Fan Huitao， Zhang Pengpeng. The Challenges for AirtoAir Missile[J].Aero Weaponry， 2017（2）： 3-7.（in Chinese）

[2] 樊會濤， 崔顥， 天光.空空導(dǎo)彈70年發(fā)展綜述[J] .航空兵器， 2016（2）： 3-12.

Fan Huitao， Cui Hao， Tian Guang. A Review on the 70Year Development of AirtoAir Missiles [J].Aero Weaponry， 2016（2）： 3-12.（in Chinese）

[3] 施坤林， 黃崢， 牛蘭杰， 等. 引信的三大基礎(chǔ)技術(shù)與發(fā)展要求[J].探測與控制學(xué)報， 2018， 40（1）： 1-4.

Shi Kunlin，Huang Zheng，Niu Lanjie， et al.Three Underlying Technologies and Developing Requirements of Fuze[J]. Journal of Detection and Control， 2018， 40（1）： 1-4.（in Chinese）

[4] 高野軍. 淺談國外引信系統(tǒng)研究現(xiàn)狀及發(fā)展趨勢[J]. 制導(dǎo)與引信， 2018， 39（1）： 1-5.

Gao Yejun. Analysis of the Research Status and Development Process of Fuze System in Foreign Countries[J].Guidance and Fuze， 2018， 39（1）： 1-5.（in Chinese）

[5] 劉靜， 趙惠昌， 高振儒.基于QMFBT的LFMCW引信信號識別[J].南京理工大學(xué)學(xué)報： 自然科學(xué)版， 2011（2）： 97-101.

Liu Jing， Zhao Huichang， Gao Zhenru.Recognition of LFMCW Fuze Signal Based on QMFBT[J].Journal of Nanjing University of Science and Technology： Natural Science， 2011（2）： 97-101.（in Chinese）

[6] 秦鵬， 呂輝， 李駿.LFMCW雷達(dá)信號處理及目標(biāo)模擬[J].中國軍轉(zhuǎn)民， 2018（1）： 62-64.

Qin Peng， Lü Hui， Li Jun.LFMCW Radar Signal Processing and Target Simulation[J]. Defence Industry Conversion in China， 2018（1）： 62-64.（in Chinese）

[7] 張宏宇， 王冬華， 李云飛， 等.一種提高LFMCW 體制雷達(dá)測距精度的方法[J].雷達(dá)與對抗， 2017， 37（1）： 9-12.

Zhang Hongyu， Wang Donghua， Li Yunfei， et al. A Method to Improve Ranging Accuracy of LFMCW Radar[J].Radar & ECM， 2017， 37（1）： 9-12.（in Chinese）

[8] 馬馳， 劉秋生， 彭朝琴.對稱三角波調(diào)頻定高引信設(shè)計與仿真[J].無線電工程， 2016， 46（1）： 61-64.

Ma Chi， Liu Qiusheng， Peng Zhaoqin.Design and Simulation of Symmetrical Triangle Wave Modulation Altimeter Fuze[J].Radio Engineering， 2016， 46（1）： 61-64.（in Chinese）

[9] 汪意焙.24 GHz汽車防撞雷達(dá)系統(tǒng)的研究與實(shí)現(xiàn)[D].成都： 電子科技大學(xué)， 2018： 15-20.

Wang Yibei.Research and Implementation of 24 GHz Automotive AntiCollision Radar System[D].Chengdu： University of Electronic Science and Technology of China， 2018： 15-20.（in Chinese）

[10] 薛定宇， 陳陽泉.基于MATLAB/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用[M].2版.北京：清華大學(xué)出版社， 2011： 129-130.

Xue Dingyu， Chen Yangquan.System Simulation and Application Based on MATLAB/Simulink[M].2nd Ed.Beijing： Tsinghua University Publishing House， 2011： 129-130.（in Chinese）

[11] 劉文彬， 朱名日， 鄭丹平， 等. 基于FPGA 的二維FFT 算法在LFMCW雷達(dá)信號處理中的應(yīng)用[J].電子器件， 2015， 38（3）： 597-600.

Liu Wenbin， Zhu Mingri， Zheng Danping， et al. Applications of TwoDimensional FFT Algorithm in Radar Signal Processing of LFMCW Based on FPGA[J].Chinese Journal of Electron Devices， 2015， 38（3）： 597-600.（in Chinese）

[12] 陳龍， 曹力.一種Zynq SoC片內(nèi)硬件加速的二維傅里葉變換[J].單片機(jī)與嵌入式系統(tǒng)應(yīng)用， 2018（2）： 36-40.

Chen Long， Cao Li.TwoDimensional Fourier Transform Based on Zynq SoC Hardware Acceleration [J].Microcontrollers & Embedded Systems， 2018（2）： 36-40. （in Chinese）

[13] 孫全， 高雅. 線性調(diào)頻信號脈壓加權(quán)處理對雷達(dá)測距影響分析[J].數(shù)字技術(shù)與應(yīng)用， 2017（4） ： 87-88.

Sun Quan， Gao Ya.Analysis of the Effect of Pulse Compression Weighting Processing of LFM Signal on Radar Ranging[J].Digital Technology and Application， 2017（4） ： 87-88.（in Chinese）

[14] 王俊鵬， 雷振亞， 王菲， 等. 基于亮點(diǎn)法的引信射頻系統(tǒng)仿真[J].微波學(xué)報， 2010（8） ： 270-273.

Wang Junpeng， Lei Zhenya， Wang Fei， et al. Simulation of RF System Based on Highlight Method[J].Journal of Microwaves， 2010（8） ： 270-273.（in Chinese）

[15] 張淑杰， 方宏業(yè).基于Simulink 的脈沖壓縮雷達(dá)信號處理系統(tǒng)仿真[J].湖北汽車工業(yè)學(xué)院學(xué)報， 2016， 30（1） ： 37-41.

Zhang Shujie， Fang Hongye.Analysis of SimulinkBased Signal Processing of Pulse Compression Radar[J].Journal of Hubei University of Automotive Technology， 2016， 30（1） ： 37-41.（in Chinese）

Abstract： In order to improve the optimal damage control ability of LFMCW fuze， the beat signal is analyzed theoretically combining with the characteristics of fuze in this paper. The factors affecting range and velocity resolution are obtained. The simulation model is established in Simulink， and the method of obtained distance and velocity by twodimensional Fourier transform is simulated. The model can simulate the echo signal and noise of the fuze， and IF signal processing process. The simulation results verify the feasibility of the algorithm and the accuracy of the model， which lays a foundation for the optimal design of fuze.

Key words： LFMCW fuze; Simulink; ranging; velocity measurement; twodimensional Fourier transform